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INTRODUCTION 

Recent developneutB in the structural synthesis area point to the 
developing use of advanced optimisation techniques, combined with sophis- 
ticated structural analysis, to provide a system which can be used for 
production application. One such system has been described recently by 
Sobieski and Bhat (1). Their system makes use of the program, CONMIN, 
as the optimizer, which is based on mathematical nonlinear programming 
techniques of feasible-usable directions in combination with a structural 
analyzer called SPAR, (2,3) which is a program of high modularity and — 
computer efficiency. These two programs have been condiined through Opti- 
mizer-to-Analyzer and Analyzer-to-Optimizer Processors and the use of a 
standard computer operating system, in this case CDC NOS. This implemen- 
tation provides a system which runs in a large computing environment (CDC 
Cyber). While this implementation of the system is quite workable, it 
has all of the problems of making use of a large computer utility, in 
that turnaround time becomes of concern, and the ability to Interact 
with the system during execution is generally lost. 

Recently, some attractive advances have been made in the computer 
area, as it relates to general structural analysis capability. The intro- 
duction of large scale mini-computers available at very reasonable 
prices, now makes it possible to perform extensive structural analysis 
on these machines. A number of computer manufacturers are now making 
what can be called maxi -mini-computers , with large central memories and 
virtual operating systems, whose speeds are in the medium range of all 
computing equipment. This type of computing hardware sells for, on the 
order of, one tenth the cost of some of the large systems, such as IBM 
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f 
I 

and CDC. Consequently, even though the processing speeds may not be 
as fast as the large scale machines, the price performance certainly 
is far better than anything we have seen previously. This availability 
has led to decentralization of computing and allows a structural 
analysis-design group to have its own computing capability under their 
direct control. Even though the processing speed may not be as fast 
as the larger systems, the turnaround time and the Interactiveness of 
this type of system is Improved tremendously. It is just this kind of 
system which NASA Langley has Installed in their structural analysis- 
design group and has Implemented the program SPAR, a large general pur- 

i 

pose linear structural analysis package, on this system. | 

The question then becomes whether or not it is viable to move the 1 

concept, as demonstrated by Sobleski and Bhat to the PRIME or similar j 

computer as a total system. The problem generated is that the overall 
optimization problem generally Involves a considerable number of re- 
analysis steps as one approaches an optimal solution in the design 
sense. Consequently, of concern is the operating speed of the analyzer 
program SPAR on the PRIME P400 hardware at Langley. It certainly appears 
that the analysis speeds must be Increased by an order of 5 to 10, in 
order for the total system concept to be implemented on the PRIME and 
to provide reasonable turnaround and Interactiveness. 

One potential for providing this Increased speed on the PRIME or 
similar hardware is that of parallel or array processing Implemented 
1 in hardware which has become available recently at very reasonable 

I 

prices. There has been considerable Interest in array processing as 
illustrated by the early machines, such as ILLIAC IV and the conmier- 
clally available CDC STAR. However, these two pieces of hardware are 
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very expensive in their design and implementation. On the other hand, 
there are now manufacturers, such as Floating Point Systems, who make 
array processors with execution speeds on the order of 7 to 10 million 
instructions per second, and which are available for a price of around 
SO to 60 thousand dollars and which can be interfaced directly to some 
of the mini-computers presently available. 

It is the potential of the connection of a Floating Points Systems 
AP120 processor to a mini-computer which needs exploration. If, indeed, 
the connection of an array processor can provide significant increase 
in execution speed for the SPAR analysis program, then the potential 
for Implementing the Sobieski-Bhat concept on a maxi-mini-computer 
would be greatly enhanced. 

This research project proposes to perform research in the area of 
simulating the effects of the use of array processor techniques within 
the program SPAR, in order to evaluate the potential speedups which may 
result. This research is possible because of the availability of an 
AP120 simulator package which now executes on a PRIME computer and the 
availability of a PDP-11/40 with an AP120B attached at RPI. Consequently 
the potential exists for simulating portions of the SPAR analysis pro- 
gram both with the PRIME simulator and the PDP-11/40-AP! 20B h rdware. 

PREVIOUS WORK 

I 

During the 1978-79 academic year, RPI has worked with the Floating 
Point Systems array processor simulator on our PRIME computer to obtain 
some general information concerning potential speedups of small pieces of 
computer code as lifted from some SPAR Processors (4) . This preliminary 
work shows results only of the floating point processor execution times. 
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but does point towards considerable potential for significant speedup 
of some of the tight loop areas of the SPAR code. Several sections of 
code from SPAR have been recoded in Floating Point System's assembly 
language and run on the AP120 simulator on the PRIME. The code is that 
taken from a section of the matrix Inversion processo'r of SPAR. A 
comparison of computation times between the PRIME P500 and the AP120 
for completing identical routines for a given joint was made, varying 
several important parameters. The results are shown in Table 1. 


! Trial 

CONRNG 

NOF 

NZERO 

CPU(PRIME) 


> 

1 

2 

3 

3 

3.00 ms 

0.07 ms 

2 

6 

3 

3 

24.24 ms 

0.80 ms 

3 

10 

3 

3 

84.80 ms 

2.26 ms 

; ^ 

15 

3 

3 

187.88 ms 

3.37 ms 

5 

h 

6 

6 

3 

90.90 ms 

2.39 ms 

r 

NOF 

- Number of 

degrees of 

freedom per Joint 


CCNRNG - Number of non-sero submatrices in a particular 
row of the stiffness matrix 

NZERO - Number of non-zero degrees of freedom for the 
given calculation (Constraint conditions) . 

TABLE 1 

These results, even though they do not include any I/O time re- 
quired to pass the data arrays between the PRIME and the AP120, pointed 
to potential for significant speedup with the use of the array processor, 

SCOPE OF PRESENT WORK 

The present work has concerned itself with a systematic approach 
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to simulating and measuring the AP120 performance In relationship to 
a nundier of SPAR processors. The past simulation data Is reevaluated 
and some additional data generated to better understand the potential 
and feasibility of using an array processor such as the AP120 to speed- 
up the analysis process for large structural systems.' 

The past simulator study measured only the execution time of the 
array processor instruction set for some selected portions of trans- 
lated SPAR code. There was not direc.. measuring of the input/output 
and data transfer times between the PRIME and the simulated AP120 con- 
nection. In this study, a limited amount of actual measurements of 
input/output and data transfer times were made using a PDP-11/40 with 
an AP120B attachment. Using these data, estimates are made as to the 
relative speedups that can be executed In a more complete Implementation 
on an array processor - maxi-mlnl computer system. 

SPAR CODE USED 

Since the previous study showed that a significant amount of man 
power was required to rewrite existing FORTRAN routines from the SPAR 
Processors in the assembler language of the AP120, it was decided to 
test those sections of SPAR code which had already been translated for 
the simulation study. This involved the actual Implementation on the 
AP120B of this code. Significant differences between the actual 
AP120B and the simulator codes had to be resolved by debugging. 

Seven SPAR subroutines or portions of SPAR subroutines have been 


implemented in the AP120 assembler language and timing Infonnation 
gathered for FORTRAN execution and AP120 assesdiler execution. The 
sections of code were t^en from five different SPAR Processors; EIG, 
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INV, SSOL, K ani H. Each piece of code will be described individually. 
EIG 


The EIG Processor is used to solve vibration and blfuncation 
buckling eigenproblems. Three short subroutines were' selected from 
EIG and portions of each were coded in AP120 assembler and implentented 
on the PDP-11/40-AP120B system. The routines converted do not repre- 
sent a significant contribution to the running time of EIG, but were 
selected for their relative ease of program conversion and the fact 
that they did use tightly coupled loops of vector floating point 
multiplications and additions. The experience gained in using these 
routines helped greatly in the conversion of more complex jperatlons. 

The first routine tested was subroutine EIGLD, which is u.ied by 
SPAR to set problem dependent parameters. The particular code con- 
verted was seven lines of FORTRAN which generates a vector with random 
valued components. Figure 1 shows these seven lines of FORTRAN code. 
In order to perform timing testing, three stand alone programs were 
developed. The first represents the FORTRAN code to be converted. 
Dimension, timing and input/output have been built around the seven 
statements. This program is called FORI. FOR and allows us to obtain 
timing Information for the FORTRAN execution. The second program de- 
veloped replaces the FORTRAN statements with FORTRAN calls to Array 
Processor utilities and to the AP assend>ler replacement code. This 
program also contains timing generating elements and is used to ob- 
tain timing information for the AP execution. This program is desig- 
nated APFORl.FOR. The third program is the AP assembler program which 
performs the operations on the AP120B. It is developed as a FORTRAN 


callable routine and is called from APFORl. Ttie asaei^ler routine is 
named APEGLD. These three progrras are presented in Appendix A. 

The second routine tested was subroutine EXPNDl, vdiich is used by 
SPAR to perform the addition of one scaled vector to another and sub- 
stitute the result into the original vector. The particular code con- 
verted was 13 lines of FORTRAN which is shown in Figure 2. The timing 
tests were performed using three stand alone programs. Progrra 
F0R2.F0R represents the 13 lines of code with dimensioning, timing 
calls and input/output added. This program provides timing informa- 
tion for the FORTRAN execution. The second program, APF0R2.F0R, re- 
places the FORTRAN statements with FORTRAN calls to AP utilities and 
the AP assembler code. This program also contains timing generating 
elements and is used to time the AP execution. The third program is 
the AP assembler routine which performs the operations on the AP120B 
and is named APXPDl.APM. These three programs are presented in 
Appendix B. 

The third routine from EIG which was tested is NEWX, which per- 
forms a number of operations on vectors, including zeroing, normal- 
ization, and multiplication by a constant. Three separate portions 
of the subroutine NEWX have been tested. The first one is five lines 
of FORTRAN used to zero a vector. This code is shown as Figure 3, 

The three programs for this portion are F0R3.F0R, APF0R3.F0R and 
APNWXl.APM. These three programs are presented in Appendix C. 

The second portion of NEWX is 11 lines of FORTRAN which performs 
a normalization operation on two vectors. This code is shown as 
Figure 4. The three programs for this portion are F0R4.F0R, 

APF0R4.F0R and APNWXl.APM. These three programs are presented in 
Appendix D. 


The third portion of NEWX Is 6 lines of FORTRMl which perfoms 
a Bulclpllcatlon of a vector by a constant. This code Is shown as 
Figure 5. The three programs for this portion are FOR5.FUR, 
APFORS.FOR and APHWX3.APM. These three programs are presented In 
Appendix E . 

INV 


INV Is the stiffness matrix inversion or reduction processor 
for SPAR. Two routines which represent over 50 percent of the 
total CPU time used in INV have been implemented. The first Is the 
major portion of the Subroutine RED which performs the core of the 
reduction process. Elghty-two lines of FORTRAN code have been coded 
in AP assembler, representing the major test implementation for this 
prcject. Figure 6 lists this code. Again three programs have been 
developed for testing purposes. They are FOR6.FOR, AFFOR6.FOR and 
APRED.AFM. These three programs are presented in Appendix F. 

The second routine from INV which has been implemented Is AFEX. 
Seventeen lines of FORTRAN have been converted to AP asseiid>ler and 
this FORTRAN is shown as Figure 7. The three prograias developed 
for testing this code are FOR8.FOR, APFORS.FOR and APAPEX.APH. These 
three programs are presented in ApperJix G. 

SSOL 


SSOL is the static solution processor for SPAR. It performs 
the matrix multiplication operations required to obtain the displace- 
ment vectors for each loading condition using the reduced matrix 
from INV. The basic routine used is MULTEX, a routine found in the 


SPMl library. Ttw •lav«n lines of FORTRAN shown in Figure 8 have 
been isipleaented by the usual three programs designated as F0R9.F0R. 
APF0R9.F0R and APMLTX.APM. These three progruss are presented in 
Appendix H. 

K & M 

The K and M processors asseid)le the system stiffness and mass 
matrices, respectively. Both the R and M processors call the routine 
TRAN6 to perform the transformation 

TK^ * GKL * TL 

for each node. Fifty-two lines of FORTRAN code representing the 
major portion of TRANS have been Implemented. This code is shown 
in Figure 9. The three programs used to test timing have been de- 
signated FORC.FOR, APFORO.FOR and APTRN6.APM. Listings of these 
three programs are presented in Appendix I. 

TEST PROGRAM 

Six of the nine program sets described above were tested ex- 
tensively on the PDP-11/40 - AP120B system. The most important im- 
plementations are F0R6, which relates to the stiffness matrix reduc- 
tion, F0R9» which relates to the matrix multiplication process and 
/ORO, which performs rotational and translational transformations 
of matrices. 

In order to obtain some comparison of the application to real 
problems, two problems were solved using SPAR. These are FUSEL and 
LUT, listings of ffhich are provided in Appendix J. The CPU timing 


result! of these runs are sho%m in Tables 2 and 3. 


It is iniN>rtant to note that the IMV proceesor uses the nost 
CPU resources, representing approximately 40% of the total CPU time 
in LUT and approximate!'' 75% of the total in FUSEL. 

TEST RESULTS 

The results of running comparisons between FORTOAN exacMtlon and 
AP execution of the sane code are susaarised in Table 4. These ratios 
include both tia execution times and the data transfer tines in all 
cases. It is this fact that accounts for the wide range of specdups 
obtained. 

In the case of the INV reduction process is^lenentation, the re- 
sult is quite favorable. Since the execution of routine RED repre- 
sents more than one half of the CPU time for the processor IHV, the 
result indicates that one could expect on the order of doubling of 
the speed in the processor using only this one small section of code 
loplesantation. 

CONCI’JSIONS 


The results of this liralted testing program add to the evidence 
presented in Reference 4 as to the appropriateness of doing a full 
blown inplementation. The evidence is not conclusive that sufficient 
speed can be obtained for an analysis-design speedup of sufficient 
magnitude to warrant full scale testing. However, the results 
generated here are for ^ very small part of the SPAR system being 
Implemented on the AP. One must keep in mind that this limited im- 
plementation has taken a considerable amount of manpower. The 


FUSEL run with SPAR on PRIME 65C 


Processor 


CPU Time (Seconds) 


TAB 

ELD 

E 

EKS 

TOFO 

K 

INV 

AUS 

DCU 

SSOL 

GSF 

PSF 


13.6 

6.9 
15.5 
79.4 
28.0 
67.1 

768.2 

2.7 

0.5 

32.3 

7.3 

7.9 


T&ble 2 


LUX run with SPAR on PRIME 750 


Processor 


CPU Time (Seconds) 


TAB 

ELD 

TOPO 

E 

EKS 

R 

INV 

M 

AUS 

SSOL 

GSF 


28.9 

43.3 

12.3 

26.9 

32.7 

76.9 

318.6 

164.6 

25.9 
33.1 

30.7 


Table 3 


Run Time Comparisons 


Program Set 

P0R2 - XPDl 
F0R4 - NWX2 
F0R6 - RED 
F0R8 - AFEX 
F0R9 > MLTX 
FORD - TRN6 


iit?RTRAN TIME 
AP TIME 

26 

19 

51 

6 

40 

22 


Table 4 


problem is that the level of sophistication of the computer programmer 
must be high enough to be able to work with the very complex assembler 
language of the AP. Even with an expert programmer, the programming 
of the AP must be described as tedious and tricky, au« to the 
parallelism Inherent in the system. 

The final conclusion is that the purchase of an array processor 
for attachment to a maxl-mlmi computer cannot yet be justified solely 
on the evidence to date in the structural analysis area. Other appli- 
cations would have to be included to justify the purchase. 
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FIGURES 


ZRB-0 

DO 1070 J-l.JT 
DO 1065 I-l.JDF 
K-INEX( I ) 

IRBPI-IRB+I 

lOiiS A(IRBPI)-RAN(0,0)*RSCALE(K) 
lOVO IRB-IRB+JDP 


Figure 1 


DO 260 J-1,NV2 

JJZ-J+J2 

DO 250 K-l.NVl 

KKZ-K+KZ 

QKJ-Q{KKZ,JJZ) 

IP(KTEST( JJ2).NE.O) GO TO 150 
IF( JJZ.NE.KKZ) GO TO 250 
QKJ-1 . O 
150 CONTINUE 

DO 200 1-1,LVEC 

200 V2(I, J)-V2(I, J)+V1(I, K)*QKJ 
250 CONTINUE 
260 CONTINUE 


Figure 2 


DO 70 N-1,NL0ADS 
IF(LSDO(N) .EQ. 0)GO TO 70 
DO 60 1-1,I.VEC 
60 V2{I,N)-0, 

70 CONTINUE 


Figure 3 


DO 1300 N-1,NL0ADS 
IF(LSDO(N).EQ.O) GO TO 1300 
SUM-.O 

DO 1100 I-1,LVEC 
1100 SUM-SUM-«-V2(I,N)*Vl(I,N) 
SUM-1 ./SQRT( ABS( SUM ) ) 
Z(N)-SUM 
DO 1200 I-1,LVEC 
V1(1,N)-V1(I,N)*SUM 
1200 V2(I,N)-V2(I,N)*SUM 
1300 CONTINUE 


Figure 4 


DO 1850 N*1,NL0ADS 

IF (LSDO(N) .EQ, 0) GO TO 1850 

SDM-2(N) 

DO 1800 I«1,LVEC 
1800 V1(I,N)-V1(I,N)*SUM 
1850 CONTINUE 


Figure 5 
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DO 1000 K"1,NZER0 
M- XABS(MAP(K)) 

DO 1000 J-l,CONRNG 
L- SUBHAP(J) 

DO 1000 1-1. NDF 
1000 S(M,I,L) 

C 

c*»«* preliminary b modification. 

DO 1400 K-l.NZERO 
M- MAP(K) 

IF(M.LT.O) GO TO 1400 
RA-BB(M,K) 

IF(RA.GT. ZEROD) GO TO 1025 

C NEX-INEX(M) 

IF(RA.LT. -ZEROD) GO TO 1015 
C KSING-KSZNG41 

C NSING-KSING 

C WRITE( lOUT, 1010 ) JOINT.NEX 

ClOlO FORMAT(49H •*» VJARNING. SYSTEM K SINGULAR. JOINT/COMPONENT-15,12) 
RA-.O 
GO TO 1030 
C1015 KNEG-XNEG41 
C NNEG-KNEG 

C IF(IPRT.LT.2) GO TO 1025 

C WR1TE( lOOT, 1020 ) JOINT, NEX 

C1020 FORMAT( 38HONEGATIVE DIAG TERM. JOINT/COMPONENT- 15,12) 

C1025 RA-l./RA 
1015 CONTINUE 
1025 CONTINUE 
1030 BB(H,K)-RA 

IF(K.EQ.NZERO) GO TO 1200 
LA-K+1 

DO 1100 L-LA,NZER0 
IA-IABS(MAP(D) 

RAB- RA*BB(IA,K) 

DO 1100 I-IA,NDFCON 
1100 BB(I,L)- BB(I,L) -RAB*BB(I,K) 

1200 IF(H.EQ.NDP) GO TO 1400 
INEXT-M41 

DO 1300 I-INEXT,NDF 
1300 BB(Z,K)- BB(Z,K)*RA 
1400 CONTINUE 
C 

ZF(CONRNG.EQ.l) GO TO 2500 


Figure 6 


DO 2100 K-l.NZERO 
M- MAP(K) 

ZP(M.LT.O) GO TO 2100 
RA- BB(M,K) 

NSUB- CONRNG 
DO 2000 I-2,C0NRNG 
NSUB- NSUB-fl 
I.S- SUBMAP(NSUB) 

C 

C**»» MODIFY Eli 

DO 1600 ICOL-l.NDF 
RAB- RA*B( ICOL,I,K) 

DO 1600 IR0W-1,IC0L 

1600 S(IR0W,IC0L,LS)- S( IROW, ICOL, LS ) -RAB»B( IROW, I,K) 

COS CALL CALCS(CONRNG,B, SUBMAP, NDF,S) 

CDC 12 CDS OMITTED 
C 

C »•« THE COMPASS ROUTINE CALCS REPLACES THE FOLLOWING ON CDC 
C 

DO 1700 IR0W-1,NDF 
1700 B(IR0W,1,K)-RA*B(ZR0W,I,K) 

IF(I.EQ.CONRNC)GO TO 2000 
C 

C***« MODIFY EIJ S 
JA-I+1 

DO 1900 J-OA,CONRNG 

NSUB-NSUB+l 

LS-SUBMAP(NSUB) 

DO 1800 IC0L-1,NDF 
DO 1800 ZR0W-1,NDF 

1800 S( IROW, ICOL, LS )-S( IROW, ICOL, LS )-B( IROW, I , K )*B( ICOL, J, K ) 
1900 CONTINUE 
2000 CONTINUE 
2100 CONTINUE 
2500 DO 2600 L-1,C0NRNG 
K- SUBMAP(L) 

DO 2600 J-1,NDF 
DO 2600 I-1,NDF 
2600 S(I,J,K)-.0 


Figure 6 (Continued) 


DO 1000 ISUB-1,NSUBS 
Xr(ISUB.EQ.l) CX> TO 600 
J-XPZX(AX(LK)) 

400 XP( J.EQ,K4(ZX:0N}) GO TO 600 
LCON-LCON41 

XF(LCON.LT.liCONX) GO TO 400 
WRITE(6,500) 

500 FORMAT( 2930*** MPILE/KMAP XNCONSISTENCY ) 
STOP 

600 LSOB-LCON4CONRNG 
K-K4( I*S0B ) 

DO 700 J-1,NDF 
DO 700 I-1,NDF 
S(I,J,K)-S(I,J,K)+AK(LKSUB) 

700 LKSUB»UCSUB-fl 
LCON-LCON+1 
1000 LX-LX-t-1 


Figure 7 


I-JLIST( 1 ) 

DO 100 ]>1,N 
DO 100 M-1,N 

100 V0UT{ L, X )-VOOT( L, X >+A( L, M, 1 ) ‘VINC M, X ) 
IF( NSUBS.lt. 2) GO TO 300 
DO 200 K-2,NSUBS 
J-JLIST(K) 

DO 200 L-1,N 
DO 200 M-1,N 

VDOT( L, X )-VDDT( L, I )+A( L, M, K )*VIN( M, J > 
200 VOUT(L,J)-VODT(L,J)+A(M,L,K)*VIN(M,I> 


Figure 8 


DO 2000 ]>1,NN00ES 
XX*ITRANS(L) 

DO 2000 K-1,L 
KX-ZTRANS(K) 

C 

DO 100 J-1,6 
DO 100 1-1,6 
GKLTL(Z,J)- .0 
100 HXL(I.J)- .0 

C FORM HKL- TK( TRANSPOSE ) *GKL *TL 

C FIRST, GKL«TL. 

C 

DO 1100 J-1,3 
DO 1100 1-1,3 
DO 1100 N-1,3 
TLMJ-T(M,J,LL) 

GaCLTL( 1 , J)- GKLTL( I, J) +S( I, M,N) *TLMJ 

GXLTL( I,J+3)- GKLTL< I,J+3) +S( I,M+3,N) *TU*J 

GXLTL(I-»-3, J)- GXLTL(I+3, J) +S(I+3, M,N) *TLMJ 
1100 GKLTL(I-<-3,J-»-3)- CaCLTL( 7+3, J+3 ) +S( I+3,M+3,N ) ‘TIiHJ 
C 

C TK( TRANSPOSE ) *( CaCL'TI. ) 

DO 1200 J-1,3 
DO 1200 1-1,3 
DO 1200 M-1,3 
TKMI-T(M,I,KK) 

HKL( I, J)- HKL( I, J) +TKMI *GJCLTI.( M, J) 
HKL( I,J+3)- Hia.( I,J+3) +TKMI*GXLTL( M,J+3) 
HKL(I+3, J)- HKL(I+3, J) +TKMI »GKLTL( M+3 , J) 
1200 HXL(Z+3,J+3)- BXL(I+3,J+3) +TKMI *GKLTL( M+3 , J+3 ) 

C 

C TRANSPOSE, IF REQ. 

IF(MAP(N).GT.O) GO TO 1400 

DO 1300 J-2,6 

JM-J-1 

DO 1300 I-1,JM 
EIJ- BRIi(Z,J) 

HKL(I,J)-Bia.(J,I) 

1300 HKL( J, I }-EZ J 
1400 LOC-IABS(NAP(N)) 

IF(NDF.LT.6) 00 TO 1600 
DO 1500 J-1,6 
DO 1500 1-1,6 

1500 B(I,J,IOC)- B(I, J,LOC)+HKL(I,J) 

GO TO 2000 
1600 DO 1700 I-1,NDF 
NROW-INEX( X ) 

DO 1700 J-1,NDF 
NCOL-INEX( J) 

1700 B( I, J,LOC)-H( I, J,LOC)+HKI,(NROW,NCOL) 

2000 CONTINUE 


Figure 9 
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FORTRAN XV V02.S-2 Tu« 04-Nov-tO 22t05t39 


PAGE 001 


CXm 05-14-00 RPX0 6066S «6660 66665 77777 11111 

C 
c 
c 

C This is Prograa FORI. FOR Which rsprsssnts s portion 
C of Subroutins EXGLD in Proesssor EXG. Used to obtain 

C timing information for FORTRAN axscution. 

c 

C Corresponding Programs arst 

C APFORl.FOR - FORTRAN of FORI with AP Calls. 

C APEGli).APN - FP120 Assssblsr Program rsplacamsnt 

C of FOOTRAN portions. 

C APEGLD.ABJ - Object Code of APEGLD.APN 

C APEGLO.SAV - Linked version of APEGLD.Am. 

c 

c •**•••»,»*********••*******••*•*•••******•*•**•* 

c 

C SUBRCKirZIOE EXGLD 

c 


0001 


DINENSXm A(6000),XNEX(6),RSCAX£(6) 

0002 


DIMENSION XTIK1( 2 ) , XTXN2( 2 ) 

0003 


DO 11 XX-1,6 

0004 


XNEX( XX )-XX 

0005 

11 

RSCALE( XX )«FL0AT( XX ) 

0006 


JT-1000 

0007 



0008 


WRITE(6,3001) 

0009 

3001 

F0RNAT( ' XNPOT THE NUMBER OF TXNE8 TO BE EXECUTED* ) 

0010 


READ(5,«) NTEST 

0011 


CPU1-SECND8( 0 . ) 

0012 


CALL GTXM(XTIMl) 

0013 


DO 270 L-1, NTEST 

0014 


xn^ 

0015 


DO 1070 J-1,JT 

0016 


DO 1065 X-1,JDF 

0017 


R-XNEX( X ) 

0018 


XRBPX-XRB4-X 

0019 

1065 

A(XRBPX )-RAN(0,0)«RSCAI£(K) 

0020 

1070 

. RB-IRB-fOnr 

0021 

270 

o-yriNOE 

0022 


CALL OTXM(XTXN2) 

0023 


CPU2-SECNDS(0. ) 

0024 


WRITE(6,50) CPUl 

0025 


WRXTE(6,50) CPU2 

0026 


CALL CVTTXM(XTXN1,XRR.XNI,X8E,XTI) 

0027 


NRXTE(6,70) XSR, XNX, XSE, XTX 

0028 


CALL CVn>XM(XTXN2,XHR.Xhl,XSE,XTX) 

0029 


HRXTE(6,70) XRR, XNX . XSE, ITX 


99990 



25 


FORTRAN XV V02.&-2 Tu« 04-NOV-8C 22t05i39 PAGE 002 

0030 50 PORNATC NUMBER OP SSC.^5 PAST NXmiGHT' • MTIS.?) 

0031 70 fORMAT( • TIME - *,I2, *1 M2, • I M2, • I M2) 

0032 CPU-CPU2-CPU1 

0033 WR1TE(6,22)CPU 

0034 22 P0RMAT(5X, *EUpS*d TiM - * ,P10.S, ‘Saconds* ) 

0035 STOP 

0036 END 
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FORTRAN IV V02.5-2 Tue 04-Nov-BO 10 s 17 j 36 


PAGE 001 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 


COH 05-14-80 RPlt 66666 66666 66666 77777 11111 

C 

C *••*************••*••****•«•••*«**»*•«*••• 

c 

C Thi» is Program APFORl.FOR which contains the AP Calls 

C as replacements for FORTRAN code in FORI. FOR. 

C Represents a portion of Subroutine EIGLD in Processor 
C £IG. Obtains timing information for AP execution. 

C 

C Corresponding Programs are: 

C FORI. FOR - Portion of Subroutine EIGLD in 

C Processor EIG. 

C APEGLD.APM - 

C APEGLD.ABJ - 

C APEGLD.SAV - 

C 

C ************************************************* 

c 

DIMENSION A( 6000 ) , INEX( 6 ) , R5CALE( 6 ) 

DIMENSION ITIM1( 2 ) , ITIM2( 2 ) , ITIM3( 2 ) , ITIM4( 2 ) 

DO 11 11-1,6 
INEX( II )-II 

11 RSCAI£( II )-FLOAT( II ) 

JT-IOOO 

JDF-6 

SEED-O. 

t»RITE(6,3001) 

3001 FORMAT( ' INPUT THE NUMBER OF TIMES TO BE EXECUTED’ ) 
READ(5,*) NTEST 
CPU1-SECNDS(0. ) 

CALL GTIM(ITIMl) 

DO 270 L-1, NTEST 
C 

C What follows are the AP calls which replace the 
C comnented code below. 

C 


0015 

CALL 

APCLR 

0016 

CALL 

APPUT( SEED, 12, 1,2) 

0017 

CALL 

APPUT( INEX(1),0, JDF,1) 

0018 

CALL 

APPUT( RSCALE( 1 ) , 6 , JDF , 2 ) 

0019 

CALL 

APWD 

D 

CALL 

GTIM( ITIM2 ) 

0020 

CALL 

APEGLD(JT,JDF) 

0021 

CALL 

APWR 

D 

CALL 

GTIM( ITIM3 ) 

0022 

CALL 

APGET( A( 1 ) , 13 , 6000 , 2 ) 

0023 

CALL 

APWD 


C 


99999 
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PAGE 002 


0024 

0025 

0026 

0027 

0028 
0029 


0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 


C End Of AP Calls, 

C 

270 CONTINUE 

CALL GTIN( ITIN4 ) 
in>U2-SECNDS( 0 . ) 

CALL CVTT1M( ITIMl , IHR, IMI , ISE, ITI ) 

VR1TE(6.70) 

WRITE(6,75) IHR, IMI, ISE, ITI 
D CALL CVTTIM( ITIM2 , IHR, IMI , ISE, ITI ) 

D VmiTE(6,80) 

D VmiTE(6,75) IHR, IMI, ISE, ITI 

D CALLCVTTIM(ITIM3,IHR,IMI,ISE,ITI) 

D WRITE(6,90) 

D WR1TE(6,75) IHR, IMI , ISE, ITI 

CALL CVTTIM( ITIM4, IHR, IMI, ISE, ITI ) 

WRITE(6,100) 

WRITE(6,75) IHR, IMI, ISE, ITI 
WRITE(6,110) CPUl 
»VRITE(6,110) CPU2 
CP0-CPU2-CPU1 
VnUTE( 6 , 120 ) CPU 

70 P0RMAT( • 5TIME AT START OF DATA INPUT - ’ ) 

80 FORMAT{ ‘STIME AT COMPLETION OF DATA INPUT AND EXECUTION’, 
1 ' START - • ) 

90 FORMAT( ' $TIMF AT END OF EXECUTION AND START OF DATA' 

1 • OUTPUT - ' I 

100 FORMAT( • STIME AT COMPLETION OF DATA OUTPUT - ’ ) 

75 FORMAT( , 12, ' i ' , 12, ' » ’ , 12, ' j ’ , 12 ) 

110 FORNAT( ' NUMBER OF SECONDS PAST NIIXdGHT - *,F15.5) 

120 FORMAT(5X, 'ELAPSED TIME - ',F10.5,' SECONDS') 

C 

C FORTRAN Code replaced - Begin. 

C 

C lRB-0 

C DO 1070 J-1,JT 

C DO 1065 1-1, JDF 

C K-INEX( I ) 

C IRBPI-IRB+I 

C 1065 A( IRBPI )-RAN(0,0)*RSCALE(K) 

C 1070 IRB-IRB-KIDF 
C 

C FORTRAN Code replaced - End. 

C 


0044 

0045 


STOP 

END 


$TZTUC APEGLD 
SENTRY APEGLD, 2 


THIS ROUTINE REPLACES A SECTION Or CODE LOCATED IN THL EIGLD 
SUBROUTINE IN THE EIG PROCESSOR 


Corresponding Programs are: 
FORI. FOR 
APFORl.FOR 

AUTHOR: K. PERSON 
DATE: FEBRUARY 1979 
Revised: L.J. Peeser and K. Nat is 
Date: Nay 1980 

USAGE 

FORTRAN: CALL APEGLD( JT, JDF) 
SPAGE 


MAIN DATA MEMORY MAP 


« 

• A ARRAY 
« 


* SEED 

* RSCALE ARRAY 

* 


« 

• INEX ARRAY 

« 


(STARTING ADDRESS) 


13 


12 


6 


0 


ARRAY DESCRIPTIONS 


A( JDF* JT ) 


REAL ARRAY. RANDOM VECTOR OUTPUT 
THAT IS RETURNED ONLY ONCE. 


SEED 


RANDOM NUMBER STORED ZN PROGRAM 
MQ«)RY. 


** RSCAIiE(6) REAL ARRAY. PASSED ONCE. 

•• 

” INEX(6) INTEGER ARRAY. PASSED ONCE. 


" S-PAD PARAMETERS 


JT 

SEQU 

0 

"TOTAL NUMBER OF 

JOINTS 

JDF 

SEQU 

1 

"NUMBER OF DEGREES OF FREEDOM 

K 

5EQU 

2 

"ADDRESS POINTER 

FOR RSCALE 

RADDR 

SEQU 

3 

"BASE ADDRESS OF 

RSCALE ARRAY 

CADDR 

SEQU 

4 

"BASE ADDRESS OF 

OUTPUT ARRAY A 

INEX 

SEQU 

5 

"BASE ADDRESS OF 

INEX ARRAY 

SADDR 

SEQU 

6 

"ADDRESS OF SEED 


INCNT 

SEQU 

7 

"INNER lOOP COUNTER 

OUTCNT 

SEQU 

10 

"OUTER LOOP COUNTER 


** THE RANDOM NUMBER GENERATOR USED IN THIS ROUTINE IS IDENTICAL 
” TO FLOATING POINT SYSTEM'S RANDOM NUMBER ROUTINE 

M 

" FORTRAN I IRB-0 

" DO 1070 

" DO 1065 1-1. JDF 

" K-INEX( I ) 

- IRBPI-IRB+I 

" 1065 A(IRBPI)-RND(0)*RSCALE(K) 

" 1070 IRB-IRB+JDF 


APEGLD: 

LDSPI SADDR; DB-12 . 

"LOAD SEED AuDH 


NOV SADDR. SADDR; SETMA 

"GET SEED 


RPSF B; DPX(0)<DB 

"GET MULTIPLIER U 


RPSF FMASK; DPX(2)<DB 

"GET FRACTION MASK 


DPX(1)<DB; DB-40000; WRTMAN 



DPX(1)<DB; DB-1015; NRTEX 

"DPX( 1 )-l . 


PMUL DPX(0).MD 

"MPY SEED*B 


MOV JT. OUTCNT 

"LOAD OUTER COUNT 


LDSPI RADDR; DB-5. 

"LOAD RSCALE ADDR -1 


LDSPI CADDR; DB-12. 
CLR INEX 
INEX 

"LOAD A ADDRESS -1 

OUTLOPt 

MOV JDF. INCNT 

"LOAD INNER COUNT 

INLOP 1 

FMUL; 

"POSH 


GTlt 


CONTI i 


CONT2S 


INC INEX; SETHA 

FNUL 

FSUB FH,DPX(l)t DPY(0)<FN 
FADD; 

LDSPl K; DB-ND 
FAND DPX(2),DPY(0)} 

ADD* X.RADDR; SETMA 
FADD ZERO.DPY(O); 

BFGT GTl 

FADD 

DPY(2)<FA 

FMUL DPY(2),MD 

FNUL 

FNUL) 

DEC INCNT 

M1<FN) INC CADDR) SETHA; 
BEQ CONTI 

FMUL DPX(0),DPY(2) 

JNP INLOP 

DEC OUTCNT 
BEQ CONT2) 

SUB JDF.INEX 
FMUL DPX(0),DPY(2) 

IMP OUTLOP 

NOP 


"GET X 
-PUSH 

"A*B ~X. SAVE A*B 
"PUSH 
"STORE K 
"ASSUME A«B>1 

-FRACTIOM EXTRACTED WITH MASK 
"GET RSCALE(K) 

"GET FRACTION DIRECTLY 
"IF B*A<1 

"GET FRACTION IMMEDIATELY 
"IF B*A>1 

"GET NEXT FA RESULT 
"SAVE FRACTION 

"MULT RSCALE*RND( 0 ) 


"DEC COUNT 

"SAVE RESULT IN A 
"GET NEXT NUMBER 


"RESET INEX 


Bs 

FMASKi 


SFP 27.0 
SFP .9999999925 
SEND 
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TYPE P0R2 

FORTRAN IV V02.5-2 Tue 04-NOV-BO 10 i 18 i 01 PAGE 001 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 


0021 

0022 

0023 

0024 

0025 

0026 
0027 


CON 05-14-80 RPI# 66666 66666 66666 77777 22222 

C 

c •***••**************•*•.**••*•***•*****••*•*••• 

c 

C This is Program POR2.POR which represents a portion of 

C Subroutine EXPNDl in Processor EIG. Used to obtain 

C timing information for FORTRAN execution, 

c 

C Corresponding Programs are; 

C APF0R2.F0R - FORTRAN of F0R2.F0R with AP Calls. 

C APXPDl.APM - FP120 Assembler Program replacement 

C of FORTRAN portions. 

C APXPOl.ABJ - Object code of APXPDl.APM. 

C APXPDl.SAV - Linked version of APXPDl.APM. 

C 

C •••***«******.********.*******«r**il,*****»**«*«*«******* 

c 

DIMENSION Q( 3 , 3 ) , V2 ( 1000 , 3 ) , Vl( 1000 , 3 ) , KTEST( 3 ) 

DIMENSION ITIM1( 2 ) , ITIN2( 2 ) 

JZ-0 

KZ-0 

NVl-3 

NV2-3 

LVEC-1000 

DO 1 JJ-1,NV1 

DO 1 II-1,NV2 

1 Q(II,JJ)-1.0 
DO 2 II-1,KV2 

2 KTEST(II)-1 
KK-LVEC 
VIRITE(6,3001) 

3001 FORMAT( * INPUT THE NUMBER OF TIMES TO BE EXECUTED' ) 
READ(5,*) NTEST 
DO 3 JJ-1,NV1 
DO 3 II-1,KK 
V1(II, JJ)-2.0 

3 V2(II,JJ)-3.0 
C 

D NR1TE(6,2001) 

D WRITE(6,2000) ( ( V2( I, J), J-1, 3 ), I-l, 3 ) 

C 

CPU1-SECNDS(0. ) 

CALL GTIN(ITINl) 

DO 270 L-1, NTEST 
DO 260 J-1,NV2 
JJZ-JtJZ 
DO 250 K-1,NV1 
KKZ-K4KZ 


11111 


FORTRAN IV 


V02 . 5-2 


Tue 04-NOV-80 10 t IB t 01 


0028 


QKJ-Q(KXZ.JJZ) 

0029 


IF(KTEST( JJZ).NE.O) GO TO 150 

0031 


IF( JJZ.NE.KKZ) GO TO 250 

0033 


QKJ-1.0 

0034 

150 

CONTINUE 

0035 


DO 200 I-l.LVEC 

0036 

200 V2(I, J)-V2(I, J)+V1(I, K)*QKJ 

0037 

250 

CONTINUE 

0038 

260 

CONTINUE 

0039 

270 

CONTINUE 

0040 


CALL GTIN(ITIM2) 

0041 


CPU2-SECNDS( 0 , ) 


C 

D WRXTE(6,2002) 

D WRITE(6,2000) ( (V2( I. J), J-1,3 ;,I-1, 3 ) 

C 


0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 


WRITE(6,50) CPUl 
WRITE(6,50) CPU2 

CALL CVTTIM(ITIM1,IHR,IMI,ISE,ITI> 
WRITE(6,70) IHR.IMI.ISE.ITI 
CALL CVTTIM(ITIM2,IHR,1MI,ISE,ITI) 
WRITE(6,70) IHR,IMI,1SE,ITI 
50 PORMAT( • NUHBER OF SECONDS PAST MIDNIGaTT 
70 FORMAT( • TIME - ’,12, ' t • , 12, ' : M2. • : ’ , 12 
CPU-CPU2-CP01 
WRITE(6,22)CPU 
22 FORMAT(5X,F16.9) 

2000 FORMAT(3(lX,3F15.6,/>) 

2001 FORMAT( ' V2 MATRIX BEF<JRE CALLS • ) 

2002 FORMAT( ' V2 MATRIX AFTER C ALLS ' ) 

STOP 

END 
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TYPE APP0R2 

FORTRAN IV V02.5-2 Tue 04-NOV-80 10 ties 27 PA(X 001 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 


0021 

0022 


COM 05-14-80 RPI« 66666 66666 66666 77777 22222 11111 

C 

C •••••*««•••«•«*•••«»•»••«••*«•«•••»«•«*«*•«•••* 

c 

C This is Program APFOR2.POR which contains the AP Calls 

C as replacements for FORTRAN code in F0R2.F0R. Represents 

C a portion of Subroutine EXPNDl in Processor EIG. Obtains 
C timing information for AP execution. 

C 

C Corresponding Programs ares 

C P0R2.F0R - 

C APXPDl.APM - 

C APXPDl.ABJ - 

C APXPDl.SAV - 

C 

C a*********************************************** 

C 

DIMENSION Q( 3 , 3 ) , V2( 1000 , 3 ) , VI ( 1000 , 3 ) , KTEST( 3 ) 

DIMENSION ITIN1( 2 ) , ITIM2( 2 ) , ITIM3( 2 ) , ITIM4( 2 ) 

JZ-0 

K2-0 

NVl-3 

NV2-3 

LVEC-1000 

DO 1 JJ-1,NV1 

DO 1 II-1,MV2 

1 Q(II,JJ)-1.0 
DO 2 II-1,NV2 

2 KTEST(II)-1 
KK-LVEC 
WRITE(6,3001) 

3001 PORMAT( ' INPUT THE NUMBER OF TIMES TO BE EXECUTED’ ) 

READ(5,*) NTEST 
DO 3 JJ-l.NVl 
DO 3 II-1,KK 
Vl(II,JJ)-2.0 

3 V2(II,JJ)-3.0 
C 

D WRITE(6,2001) 

D WRITE(6,2000) ( ( V2( I, J), J-1, 3 ), I-l, 3 ) 

C 

CALL GTIM(ITIMl) 

DO 270 L-1, NTEST 
C 

C What follows are the AP calls which replace the 

C commented code below, 

c 


0023 


K-NV1*MV2 


FORTRAN tv 


V02.5-2 


Tu« 04-Hov-eO 10 1 18 1 27 


0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 


0035 

0036 

0037 

0038 

0039 

0040 


0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 


KBASE-0 
CALL APCLR 

CALL APPUT(1CTEST(1),0,NV2,1) 

CALL APPUT(Q(1,1),100,X,2) 

CALL APPUT(V1(1,1),109,3000,2) 

CALL APPUT( V2( 1 , 1 ), 3109 , 3000 , 2 ) 

CALL APWD 

D CALL GTIN( ITIM2 ) 

CALL APXPD1(KZ,JZ.NV2,NV1,LVEC, 100, 0,109, 3109) 
CALL APWR 

D CALL GTIM( ITIN3 ) 

CALL APGET(V2( 1,1 >,3109,3000, 2) 

CALL APND 
C 

c End of AP call*, 

c 

270 aWTINOE 

CALL GTIH(ITIN4) 

C 

WRITE( 6,2002) 

WRITE{6,2000) ( < V2( 1 , J >, J-1, 3 ). I-l, 3 ) 

C 

CALL CVTTIM(ITIM1,IHR,IM1,ISE,ITI) 

WRITE(6,70) IRR,ZMI,ISE,ITZ 
D CALL CVTTIN( 1TIN2 , IBR, INI , ISE, ITI ) 

D WRITE(6,75) IHR, IMI , ISE, ITI 

D CALL CVTTIN( ITIM3 .IHR, INI , ISE , ITI ) 

D WRITE(6,80) IHR, INI , ISE, ITI 

CALL CVTTIN( ITIM4, IBR, INI , ISE, ITI ) 

WRITE(6,85) IBR, INI, ISE, ITI 
70 F0RMAT( * TINE AT START OF DATA INPUT - • ,12, • J 
l,I2,'i M2,-iM2) 

75 FORNAT( ' TINE AT CCBIPLETION OF DATA INPUT AND 
1,’EXECUTICW START - ' , 12, * i * , 12, * i ’ , 12, * i • , 12 ) 
80 FORNAT( ' TIKE AT END OF EXECUTION AND START ' , 
I’OF DATA OUTPUT - • , 12, ’ j * , 12, ’ « * , 12, ' t ’ , 12 ) 

85 FORMAT( ' TINE AT CONPLETKRi OF DATA OUTPUT - ' 
1I2,’«M2,*|’,I2,’«*,I2) 

2000 F0RNAT(3(1X,3F15.6,/)) 

2001 FOIO<AT( ' V2 MATRIX BEFORE CALLS ‘ ) 

2002 FOIOIAT( ' V2 MATRIX AFTER CALLS ’ ) 

STOP 

END 
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4TZTLE M>)0^D1 
GENTRY MPXPDl, 9 

H 

" THIS K>UTINE PERFORMS A SECTION OP CXMX IN THE EXPNOl SUBROUTINE 
" LOCATED IN THE EIG PROCESSOR. 


** AUTHOR! K. PERSON 
» DATE) MARCH 1979 

” Revised! L. J. Peeser and K. Matis 
” Date! May 1980 j 


" USAGE 

" PORTRAN! CALL APXPD1(KZ, JZ,NV2,NV1,LVEC,QBASE,ICBASE,V1BASE,V2BASE) 

M 

- all parameters are INTE<XRS] 

w 

" SPA(Z 

" ^MAIN DATA MEMORY MAP 


V2 ARRAY 


VI ARRAY 


Q ARRAY 


KTEST ARRAY 
K R Y 


(STARTING A1»RESS) 


100-Hl*N-t-LRl 


1004N*N 


100 


" S-PAD PARAMETERS 


KZ 

SEQU 0 

“ NTEGER ADDRESS POINTER 


JZ 

SEQU 1 

"INTEGER AIH>RESS POINTER 


NV2 

SEQU 2 

"OUTER LOOP COUNT 


NVl 

SEQU 3 

"INNER LOOP CTUT 


LVEC 

SEQU 4 

"VECTOR LENGTH 

1 

QBASE 

SEQU S 

"BASE ACORESS OP Q ARRAY 

< 

V 


1 



QADDR 

SEQU 

5 

XBASE 

SEQU 

6 

XAIR}R 

SEQU 

6 

CRT 

SEQU 

7 

VIBASE 

SEQU 

7 

V2BASE 

SEQU 

10 

VIADDR 

SEQU 

11 

V2AIOR 

SEQU 

12 

V2FADR 

StQO 

13 

J 

SEQU 

14 

K 

SEQU 

15 

VALUE 

SEQU 

16 

VITEHP 

SEQU 

17 


"ADDRESS POINTER FOR Q 
"BASE ADMUESS OP KTE8T ARRAY 
"ADDRESS POINTER OF KTEST 
"INNER HOST COUNTER 
"BASE ADDRESS OF VI ARRAY 
"BASE ADDRESS OF V2 ARRAY 
"ATORESS POINTER OF VI 
"ADDRESS POINTER OF V2 
"AIX>RESS POINTER FOR V2 TARGET 
"OUTER LOOP COUNT 
"INNER LOOP COUNT 
"TENPORARY STORAGE 
"BASE POINTER FOR VI 


" FORTRAN! DO 260 J-1,NV2 
" JJZ-J-fJZ 

M 

APXPDli CLR J 

DEC QBASE "LOOP SET UP 

DPY(0)<SPFN| HOV VIBASE , VIBASE "STORE VI BASE IN DATA PAD 

m 

" OUTER LOOP 

m 

LOOPlt CLR K 
INC JZ 

LDSPI VIBASE) DB-DPY(O) 

NOV VIBASE, VITEHP) 

DPY(0)<SPFN 


FIRST INNER LOOP 

FORTRAN) DO 250 K-1,NV1 
KKZ-K-^KZ 

IF (KTEST(JJZ) ,NE. 0)00 TO ISO 
IF (JJZ .NE. nCZ) GO TO 250 
QKJ-1.0 
150 amriNUE 


L00P2I NOV KADDR,KADDR) SETNA 
INC KZ 

INC QADDR; SETNA 
LDSPI VALUE) DB^(D 
MOV VALUE, VALUE 
BEQ aWTl 


"GET KTEST(JJZ) 

"GET KZ+K 

"GET g<ICZ-fZ,JZ-fJ) 

"STORE KTEST(JJZ) 

"TEST FOR KTE8T(JJZ)*0 


"GET JZ+J 
"RESTORE VI BASE 
"LOAD BASE POINTER FOR VI 
"SAVE VI BASE IN DATA PAD 
"THIS FREES SPAD 7 


Di>X(0)<MD 
JMP 8TMrF3 
CONTI t SUBt JZ.KZ 
WB CONT2 
RPSF ONE 

DPX(0)<DB ”8ET QJK< 

M 

” DO THE 200 wot CJOiCUlATZONS 

M 

” FORTRAN I DO 200 I-1,LVEC 

" 200 V2(l,J)-V7(l.J)-fVl(I,K)*t^J 


STXR.J.- WJ'f V1TENP,V1A1»R; SETMR 
NOV LVEC.CNT 
NOV V2BASE,V2AIX>R; SETHX 
nOIL DPX(0),ia>} 

NOV Vru *>R,V2PMm 

PMOL; 

raC V2FADR 

LOOPS: FMUL 

INC VIADDR) SEIMR) 

FJUX) FN,MD 
FADD 

INC V2A0DR; SETNA 
FNUL DPX(0),KD} 

DEC CHT 

ENDS: INC V2FADR; SEIMA; MI<FA} 

FMIL; 

BNE LOOPS 


** FORTRAN: 250 CONTINUE 


a»IT2 : INC K 

SOUt K,NV1 
BEQ CONTSi 

AI8> LVEC.VITENP 

JMP LOOP2 


” FORTRAN: 260 CONTINUE 

M 

ONITS: INC J 

INC KADim 
SUBt J.NV2 
BEQ CMIT4I 

AZE> LVEC.V2BASE 

JMP LOOT! 

ONE: SFP 1.0 


"STORE QJK-Q(KKZ.JJZ) 
"GOTO INNER MOST LOOP 
"TEST JJZ-XXZ 


1.0 


"CET VI 

"LOAD INNER COtMT 
"GET V2 
"DO QJK«V1 

"LOAD TARGET ADDRESS 
"PUSH 

"LOOP SETUP 
"PUSH 

"GET NEXT VI 
"DO V24V1*QJK 
‘PUSH 

GET NEXT V2 
"DO V1*0JK 
"TEST CNT-O 
"STORE RESULT 
"PUSH 


"TEST INNER LOOP 
"lOAJUST VI MX)RESS 


"INC OUTER LOOP COUNT 
"INCREMENT KTEST AIM>RESS 
"TEST OUTER LOOP 

"REAJUST V2 ADDRESS 


39 


C0RT4I RETUIW 
SEND 
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APPENDIX C 


Listings of: 


F0R3.F0R 

APF0R3.F0R 


APNWXl.APM 


A1 


.TYPE F0R3 

FORTRAN IV V02.5-2 Tue 04-NOV-80 10 j 18 i 59 PAGE 001 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 


COH 05-14-80 RPIt 66666 66666 66666 77777 22222 

C 

C 

c 

C This is Program F0R3.F0R which represents a portion 
C of Subroutine NEWX in Processor EIG. Used to obtain 
C timing information for FORTRAN execution. 

C 

C Corresponding Programs aret 

C APF0R3.F0R - FORTRAN of F0R3.F0R with AP Calls, 

c APNWXl.APN - FP120 Assembler Program replacement 

C of FORTRAN portions. 

C APNWXl.ABJ - Object code of APNWXl.APN. 

C APNWXl.SAV - Linked version of APNWXl.ABJ. 

c 

c *•*****•»**••***»****•**•*****»•**•*********•»********* 
c 

DIMENSION LSD0(3),V2(2000,3) 

DIMENSION ITIN1(2).ITIN2(2) 

NLOADS-3 
LVEC-20p0 
DO 10 li-l,NLOADS 
10 LSD0(II)-II 
KNTLS-0 
WRITE(6,3001) 

3001 F0RMAT( * INPUT THE NUMBER OP TIMES TO BE EXECUTED* ) 
READ(5,*) NTEST 
CPU1*SECNDS( 0 . ) 

CALL GTIM( ITIMl ) 

DO 270 L-1, NTEST 
DO 70 N-l,NLOADS 
IP(LSDO(N) .EQ. 0)G0 TO 70 
DO 60 I-l.LVEC 
60 V2(I,N)-0. 

70 CONTINUE 
270 CONTINUE 

CALL GTIM( ITIM2 ) 

CPU2-SECNDS( 0 . ) 

WRITE(6,50) CPUl 
WRITE(6,50) CPU2 

CALL CVTTIM( ITIMl, IHR,IMI.ISE,ITI) 

WRITE(6,80) IHR,IMI,ISE,1TI 
CALL CVTTIM( ITIM2 , IHR, IMI , ISE , ITI ) 

WRITC(6.80) IHR, IMI, ISE, ITI 

50 F0RMAT( ' NUMBER OP SECONDS PAST MIDNIGHT - ',P15.7) 

80 F0RMAT( • TIME - ' , 12, • i ’ , 12 , ' j ' , 12, • i * , 12 ) 

CPU-CPU2-CPU1 
WRITE(6,22 )CPU 
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0033 22 F0RMAT(5X, ‘TIME-* ,ri6.9) 

0034 STOP 

0035 END 


FORTRAN IV 


V02.5-2 


Tue 04-Nov-eO 10 1 19 1 23 


PAGE 001 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 


CON 05-14-80 RPIt 66666 66666 66666 77777 22222 44444 

C 

C *»****••*«******..*«****••**••*.«••*••**•»••*•*••••• 

c 

C This is Program APF0R3.F0R which contains the AP Calls 

C as replacements for FORTRAN code in FOR3.FQR. Represents 

C a portion of Subroutine NEWX in Processor EIG. 

C 

C Corresponding Programs arei 
C FOR3.FOR - 

C APNWXl.APN - 

C APNWXl.ABJ - 

C APNNXl.SAV - 

C 

C ••*•****«******••*•****••*••*•*****•*••*»***••**•***• 

C 

DIMENSION ^ ),V2( 2000. 3 ) 

DIMENSION IT1N1( 2 ) , ITIM2( 2 ) , ITIN3( i . ITIM4( 2 ) 

NLOADS-3 
LVEC-2000 
DO 10 II-l.NLOADS 
10 I.SDO(II)-II 
KNTLS-0 
VnUTE(6,3001) 

3001 FORMAT( ' INPUT THE NUMBER OF TIMES TO BE EXECUTED* ) 

READ(5,*) NTEST 
CPU1-SECNDS( 0 . ) 

CALL GTIM(lTIMl) 

DO 270 L-1, NTEST 
CALL APCLR 

CALL APPUT(LSDO(1),0,NLOADS,1) 

CALL APND 

D CALL GTIN(ITIM2) 

CALL APNWXKNLCADS.KNTLS.LVEC.NLOADS) 

CALL APNR 

D CALL GTIM(ITIM3) 

CALL APGET(V2(l,l).NLOADS,6O0O,2) 

CALL APWD 
270 CONTINUE 

CALL GTIN(ITIM4) 

CPU2-SECNDS(0. ) 

CALL CVTTIM(ITIM1,IHR,IMI,ISE,ITI) 

WRITE(6,70> 

WRITE(6,75) IHR.INI.ISE.ITI 
D CALL CVTTIM(1TIM2, IHR.INI.ISE.ITI) 

D WRITE(6,80) 

D WRITE(6,75> IHR.INI.ISE.ITI 

D CALL CVTTIM(ITIM3.IHR.IMI.1S1.ITI) 
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PAGE 002 


0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 


D WR1TE(6,90) 

D WRITE(6,75> IHR, IMI , ISE, ITI 

CALL CVTTIM(ITIM4,IHR,IMI,ISE,ITI) 

VfRlTE(6,100) 

WR1TE(6,75) IHR, IMI, 13E, ITI 
WRITE(6,110) CPUl 
WRITE(6,110) CPU2 
CPU-CPU2-CP01 
VIRITE(6,120) CPU 

70 F0RMAT( * STIME AT START OF DATA INPUT - * ) 

80 FORKAT( 'STIME AT COMPLETION OP DATA INPUT AND EXECUTION', 
1 ' START - ' ) 

90 FORMAT( 'STIME AT END OF EXECUTION AND START OP DATA' 

1 ' OUTPUT - • ) 

100 FORNAT( ' STIME AT COMPLETION OF DATA OUTPUT - ' ) 

75 PORMAT( ' + ',12, • j ',12, ' I ' ,12, ' : ' ,12) 

110 PORMAT( ' NUMBER OF SECONDS PAST MIDNIGHT - ',P15.5) 

120 FORMAT(5X, 'ELAPSED TIME - ',F10.5,' SECONDS') 

C 

C FORTRAN Code replaced - Begin 

c 

C DO 70 N-1,NL0ADS 

C IF(LSDO(N) .EQ. 0)GO TO 70 

C DO 60 I-1,LVEC 

C 60 V2(I,N)-0. 

C 70 CONTINUE 

C 

C FORTRAN code replaced - End 
C 


0041 

0042 


STOP 

END 


$T1TI£ APNWXl 
SENTRY APNWXl, 4 


THIS ROUTINE PERFORMS A SECTION OF CODE IN THE NEWX SUBROUTINE 
LOCATED IN THE EIG PROCESSOR 

AUTHOR; K. PERSON 
DATE: MARCH 1979 

Revised: L. J. Feeser and K. Matis 
Date: May 1980 

USAGE 

FORTRAN: CALL APNWXl (NL0ADS,KNTLS,LV£C,V2ADDR) 

SPAGE 


—MAIN DATA MEMORY MAP 


*«*■»*****■»*«***••****»***** (STARTING ADDRESS) 

* * 

* 2 ARRAY * 

* * 

************************** NL0ADS+LRM+2*(LVEC*NL0ADS ) 
« « 

* VI ARRAY * 

♦ 

***»***«**♦♦***«♦***«**«« NLOADS+IiRM+LVEC*N1jOADS 

« 

V2 ARRAY * 

* 

*****#**»********«*****iM( NIjOADS+LRM 

* 

B ARRAY * 

««**«**«*«*«««*«*«*«*««*« NLOADS 

* 

LSDO ARRAY * 

it 

itHUHiiltitititititititltlcitltltitltititititii Q 


—ARRAY DESCRIPTIONS- 


Z( NLOADS ) 


REAL ARRAY. NEVER USED 


•• Vl(LVEC«NIiOADS) . , . REAL ARRAY. PREVIOUSLY PASSED. 

M 

" V2( LVEC*MLOADS ) . . . REAL ARRAY. RETURNED ONCE. 


H 

" B(LRM) REAL ARRAY. NEVER USED. 

” LSDO(NLOAOS) INTEGER ARRAY, PREVIOUSY PASSED. 


" S-PAD PARAMETERS 


NLOADS 

$EQU 

0 

"NUMBER OF LOADS 

KNTLS 

$EQU 

1 

"INTEGER PASSED 

LVEC 

$EQU 

2 

"TOTAL JOINTS ‘DEGREES OF FREEDOM 

V2ADDR 

SEQU 

3 

"ADDRESS POINTER FOR V2 

LSDO 

$EQU 

4 

"INTEGER MAP 

VALUE 

$EQU 

5 

• xEMPORARY STORAGE 

OUTCNT 

SEQU 

6 

"OUTER LOOP COUNT 

INCNT 

SEQU 

7 

"INNER LOOP COUNT 


" FORTRANs DO 70 N-l.NLOADS 
" IF(LSDO(N) .EQ. 0 )G0 TO 70 

" KNTLS-KNTLS+1 

” DO 60 1-l.LVEC 

" 60 V2(I,N)-0. 

" 70 CONTINUE 


DEC V2ADDR 
MOV NLOADS.OUTCNT 
CLR LSDO 
DEC LSDO 

OUTLOP: INC LSDO) SETMA 

MOV LVEC,INCNT 
NOP 

LDSPI VALUE) DB-ND 
MOV VALUE, VALUE 
BNE CONTI 
ADD LVEC,V2ADDR 
JMP CONT2 


"SET UP V2ADDR FOR LOOP 
"LOAD OUTER COUNT 
' i LEAR ADDRESS POINTER 

"GET LSDO(N) 

"LOAD INNER COUNTER 

"VALUE-LSDO( N ) 

"TEST VALUE-O 

"AJUST V2ADDR 


CONTI I 


INC XNTLS 


LOOPS 

DEC INCNT 



INC V2ADDR;SETMAsMI<ZERO;BNE LOOP 

"CLR V2 

CONT2S 

DEC OUTCNT 

"TEST OUTER LOOP 


BNE OUTLOP 



SEND 



APPENDIX D 


Listings of: 


FORA. FOR 
APFORA.FOR 


APNWX2..\PM 
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•m>E F0R4 

FOKTRMI IV V02.5-2 Tue 04-Nov-BO 10 i 20 >24 


001 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0CI2 

0013 

0014 


0015 

0016 

0017 

0018 
0019 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 
0029 


COM 05-14-80 RPI8 66666 66666 66666 77777 22222 

C 

C **«..**.*.••*.»**.*«***..*«**«•***•*•*•••*••**•*• 

c 

C This is Frograun F0R4.F0R which represents a portion of 
C Subroutine NEVfX in Processor EIG. Used to obtain 
C timing information for FORTRAN execution. 

C 

C Corresponding Programs are> 

C APF0R4.F0R - FORTRAN of F0R4.F0R With AP Calls. 

C APNWX2.APM - FP120 Assembler Program replacement 

C of FORTRAN portion. 

C APNWX2.ABJ - Object code of APNNX2.APM. 

C APNWX2.SAV - Linked version of APNWX2.ABJ. 

C 

C **«**•**♦«***•**•***•••«*******«*••**•****«********«*.««* 

c 

DIMENSION V2(1000,4),V1(1000,4),LSOO(4),Z(4) 

DIMENSION ITIM1(2),ITIM2(2) 

NLOADS-4 

LVEC-1000 

DO 11 II-1,NL0ADS 

11 LSD0(II)-II 
DO 12 JJ-1,4 
DO 12 11-1,1000 
Vl(II,JJ)-2. 

V2(II,JJ)-1. 

12 CONTINUE 
VnUTE(6,3001) 

3001 F0RMAT( ’ INPUT THE NUMBER OP TIMES TO BE EXECUTED* ) 
READ(5,*) NTEST 
C 

D WRITE(6,90) 

D WRITE(6,95) ( ( Vl( I , J), J-1, 4), I-l, 3 ) 

CFU1-SECNDS( 0 . ) 

CALL GTIM(ITIMl) 

DO 270 L-1, NTEST 
DO 1300 N-1,NL0ADS 
IF(LSDO(N).EQ.O) GO TO 1300 
SUM-.O 

DO 1100 I-l.LVEC 
1100 SUM-SUM+V2(I,N)* *1(I,N) 

SUM-1 . /S(3RT( ABS( SUM ) ) 

Z(N)-SUN 

DO 1200 I-1,LVEC 
V1(I,N)-V1(I,N)*SUM 
1200 V2(I,N)-V2{I,N)*SUM 
1300 CONTINUE 
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50 


FORTRAN IV V02.5-2 Tu« 04-Nov>80 10 t 20 i 24 PA(^ 002 

0030 270 CONTINUE 

0031 CAU. GT1M( IT1H2 ) 

0032 CPU2-SECNDS(0. ) 

C 

D WRITE(6,100) 

D WRITE(6,95) ( ( Vl( I, J ), J-1, 4 ) , 1*1, 3 ) 

D NRITE(6,110) 

D WRITE(6,95) ( ( V2( I , J), J-1, 4 ), I-l, 3 ) 

U033 WRITX(6,50) CPUl 

0034 WRITE(e,50) CPU2 

0035 CALL CVTTIH(1TIM1,IHR,IMI,ISE,ITI) 

0036 WR1TE(6,70) IHR, IMI , ISE, ITI 

0037 CALL CVTTIN(ITIM2,IHR,IMI,1SE,ITI) 

0038 WR1TE(6,70) IHR, IMI , ISE, ITI 

0039 50 PORMAT( ' NUMBER OF SECC»n)S PAST MIIBIKfflT - *,F15.7) 

0040 70 FORMAT( • TIME- * ,12, ' j ’ ,12, ’ i ’ ,12, * i ’12) 

0041 90 FORMAT( ' VI MATRIX BEFORE CALL* ) 

0042 95 FORMAT(3(lX,4F15.7,/)) 

0043 100 FORMAT( ' VI MATRIX AFTER RETURN* ) 

0044 110 FORMAT( * V2 MATRIX AFTER RETURN* ) 

0045 CPU-CPU2-CPU1 

0046 VmiTE(6,23)CPU 

0047 23 FORMAT( 5X, * TIME- * , F16 . 9 ) 

0048 STOP 

0049 END 
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TYPE APF0R4 

FORTRIUS IV V02.S-2 Tu« 04-Nov-BO 10 1 20 1 48 PAGE 001 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 


0015 

0016 

0017 

0018 

0019 

0020 
0021 

0022 

0023 


COM 05-14-80 RPI8 66666 66666 66666 77777 22222 

C 

c *••«•••••.«•***••*•••*****•••.***•*••••••**••**•*•*•• 

c 

C This is Program APP0R4.P0R which contains ths AP Calls 

C as rsplacsments for FORTRAN cods in P0R4.P0R. Rsprssents 

C a portion of Subroutino NEWX in Procsssor EZG. Obtains 
C timing information for AP sxscution. 

C 

C Corrssponding Programs arst 

C P0R4.P0R 

C APNNX2.APM 

C APNNX2.ABJ 

C APMNX2.SAV 

C 

c *•**•*••.«.**«***•••«.*.**•****«******••*•••***•*«*••• 

c 

DIMENSION V2( 1000, 4 ) . ^a( 1000, 4 ) , LSD0( 4 ), Z( 4 ) 

DIMENSION ITIM1( 2 ) , a TIM2( 2 ) , ITIM3( 2 ) , ITIM4( 2 ) 

NZiOADS-4 

LVEC-1000 

DO 11 II«1,NL0ADS 

11 LSD0(IZ)-II 
DO 12 JJ-1,4 
DO 12 ii-i,io:,o 
Vl(ll,JJ)-2. 

V2(II,JJ)-1. 

12 CONTINUE 
WRITE(6,3001) 

3001 PORMAT( ' INPUT THE NUMBER OP TIMES TO BE EXECUTED* ) 
READ(5,*> NTEST 
D WRITE(6,100) 

D WRITE(6,90) ( ( Vl( I , J ), J-1, 4 ), I-l, 3 ) 

C 

C THE FOLLOWING CALLS REPLACE THE FORTRAN CCOOfENTED BELOW 
C 

CALL Gl'ZM(ITZMl) 

DO 270 L-1, NTEST 
CALL APCLR 

C\LL APPOT(LSDO(1),0,NLOADS,1) 

CALL APPUT(V1(1,1),100,4000,2) 

CALL APPUT(V2( 1,1 >,4101,4000, 2) 

CALL APND 

D CALL GTIN(ITIN2) 

CALL APNMX2( 100, 4101, LVEC,NLOADS, 8101) 

CALL APWR 

D CALL GTIN(I7IN3) 

CALL AP<aET(Vl(l,l),100,4000,2) 
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0024 
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PAGE 002 


0025 

0026 

0027 

0028 


CALL APGET(V2(1,1),410;.4000,2) 
CALL APWD 
270 CONTINUE 

CALL GTIM(ITIH4) 

C 

C FORTRAN cod« replaced - Begin 

C 

C DO 1300 N-l.NLOADS 

C IF(LSDO(N).EQ.O) GO TO 1300 

C SUM-.O 

C DO 1100 I-1,LVEC 

C 1100 SUM-SUM+V2(I,N)*V1(1,N) 

C SUM-1 . /SQRT( ABS( SUM > ) 

C Z( N )-SUN 

C DO 1200 I*1,LVEC 

C V1(1,N)-V1(I,N)*SUM 

C 1200 V2(I,N)-V2(I,N)*SUM 
C 1300 CONTINUE 
C 

C FORTRAN Code replaced - End 

C 


0029 CALL CVTTIM(ITIK1,IHR,IMI.ISE,IT1) 

0030 WR1TE(6,70) IHR, IMI. ISE« ITI 

D CALL CVTTIM( ITIM2,IHR,IMI,ISE,ITI) 

D WR1TE(6,75) IHR, IMl , ISE, ITl 

D CALL CVTTTN(ITIN3,IHR,IMI,ISE,ITI) 

D WRITE(6,80) IHR, IMI , ISE, ITl 

0031 CALL CVTTIN(ITIM4,IHR,IMI,ISE,ITI) 

0032 WRITE(6,85) IHR, INI , ISE, ITI 

0033 70 FORMAT( ' TIME AT STAOT OF DATA TRANSFER - ‘ , 

112, ' j M2, ' I M2, • I *12) 

0034 75 FOR»5AT( ' TIKE AT COMPLETION OF DATA INPUT AND ' 

1, *E:3'_UTI0N START - M2, ' j M2, ' t M2, * i M2 ) 

0035 80 FORMAT( ' TIME AT END OF EXECUTION AND START ' , 

l*OF DATA OUTPUT - M2, ' t M2, • t M2, ’ J M2 ) 

0036 85 FORMAT( ' TIME AT COMPLETION OF DATA OUTPUT - ' , 

112, • : M2, • j M2. • : M2) 



D 

WRITE(6,95) 



D 

HRITE(6,90) 

((V1(I,J),J-1,4),I-1. 


D 

WRITE(6,97) 



D 

WRITE(6,90) 

((V2(I.J).J-1,4),I-1, 

0037 

90 

rontnT( 3( IX, 

,4F15.7./)) 

0038 

95 

FORMAT( • VI 

MATRIX AFTER RETURN' ) 

0039 

97 

FORMAT( * V2 

MATRIX AFTER RETURN' ) 

0040 

100 

FORNAT( • VI 

MATRIX BEFORE CALL' 

0041 


STOP 


0042 


END 



(TITLE APNWX2 
(ENTRY APNWX2, 5 
(EXT DIV 
(EXT SQRT 


THIS ROUTINE PERFORMS A SECTION OP CODE IN THE NEWX SUBROUTiNE 
LOCATED IN THE FIG PROCESSOR 


AUTHORi K. PERSON 

DATE: FEBRUARY 1979 

Revised: L. J. Peeser and K. Hatis 

Date: May I960 

USAGE 

FORTRAN: CALL APNWX2( VlBASE,V2BASE,LVEC,NLOAD,ZBASE) 
ALL PARAMETERS ARE INTEGERS] 

(PAGE 


MAIN DATA MEMORY MAP 


Z ARRAY 


VI ARRAY 


V2 ARRAY 


B ARRAY • 

* 

4 * 

LSDO ARRAY • 


(STARTING ADDRESS) 


NLOADS+LRM+2 •( LVEC*NLOAOS ) 


NL0A0S4LRM-fLVEC*NL0ADS 


NLOADS-t-LRM 


NLOAOS 


0 


•• ARRAY DESCRIPTIONS 

M 

” Z(NLOADS) REAL ARRAY. NEVER PASSED. 

M 

" Vl(LVEC*NLOADS). . . REAL ARRAY. PREVIOUSLY PASSED. 
'• RETURNED ONLY ONCE. 

fl 

" V2(LVEC*NLOADS). . . REAL ARRAY. RETURNED ONCE. 


B(LRM) REAL ARRAY. NEVER USED. 

ff 

■■ LSDO(NIiOADS) INTEGER ARRAY. PREVIOUSY PASSED. 

ff 


” S-FAD PARAMETERS 


VIBASE $EQU 0 
V2BASE $EQU 1 
LVEC SEQU 2 
NLOAD SEQU 3 
ZBASE SEQU 4 
VIADDR SEQU 5 
V2ADDR SEQU 6 
CNT SEQU 7 
LSADDR SEQU 10 
OUTCNT SEQU 10 
VALUE SEQU 11 
VIDEST SEQU 12 
V2DEST SEQU 13 


" FORTRAN: DO 1300 N-1,NL0ADS 
" IF(LSDO(N) .EQ. 0) GO TO 

" SUM-0 . 


"BASE OF VI ARRAY 
"BASE OF V2 ARRAY 

"TOTAL JOINTS * DEGREES OF FREDDOM 
"NUMBER OF LOADS 
"ADDDRESS OF Z ARRAY 
"ADDRESS OF VI 
"ADDRESS OF V2 
"INNER COUNTER 
"ADDRESS OF LSDO ARRAY 
"OUTER LOOP COUNTER 
"TEMPORARY STORAGE 
"VI ADDRESS FOR WRITE 
"V2 ADDRESS FOR WRITE 


1300 


APNWX2; CLR OUTCNT 
LOOPl: MOV LSADDR, LSADDR; SETMA 

DPX(0)<DB; DB-ZERO 
NOP 

LDSPI VALUE; DB-MD 
MOV VALUE, VALUE 
BNE CONTI 
JMP C0NT2 


"CLEAR 1300 LOOP COUNTER 
"GET LSDO(N) 

"CLEAR DPX(0)-SUH 

"STORE LSDO(N) 

"TEST LSDO(N) 




1100 LOOP CALCULATIONS 


55 


" FORTRAN I DO 1100 I-l.LVEC 

" 1100 SUM-SUM4V2(1,N)*V1(I,N) 


-DPX(O) SHOULD BE ZERO 


CONTI: 


L00P2: 


MOV VIBASE.VIADDR; SETMA 

MOV LVEC.CNT 

MOV V2BASE.V2ADDR; SETMA 

DPX(1)<MD 

INC VIADDR; SETMA 

FMUL DPX(1),MD 

INC V2ADDR; SETMA; 

FMUL 

DPX(1)<MD; 

FMUL 

INC VIADDR; SETMA; 

FADD FM,DPX(0) 
FMUL DPX(1),MD; 

FADD; 

DEC CUT 


INC V2ADDR; SETMA; 

FMUL;DPX(0)<FA;BNE L00P2 


"GET V(1,N) 

"LOAD INNER COUNTER 
"GET V2(1,N) 

"SAVE VI 
"GET V1(2,N) 

"DO V1*V2 
"GET V2(2,N) 

"PUSH 

"SAVE V1(LVEC,N) 
"POSH 

"GET V2(LVEC,N) 
"V2*V1+SUM 
"V1*V2 
"PUSH 
"TEST LOOP 
"GET V2(LVEC,N) 
"SAVE SUM 


" FORTRAN: SUM-1 . 0/SQRT( ABS( SUM) ) 
" 2( N )-SUM 


JSR SQRT 
UVTMA;DB-IONE 
NOP 

DPY(0)<TM 
JSR DIV 


"SQRT(ABS(SUM)> 


"SET DPY-1 
"DPY(l) DIV DPX(O) 


OR THE INVERSE OF SQRT( ABS( SUM) ) 

ADD* OUTCNT,ZBASE; SETMA; MI<DPX(0) "STORE Z(N)-SUM 


THE RESULT IS IN DPX( 0 )-SUM 
DO THE 1200 LOOP CALCULATIONS 


FORTRAN: DO 1200 I-1,LVEC 

Vl( I,N)-V1(1,N)*S0M 
1200 V2( I,N)-V2( I,N)*SUM 
1300 CONTINUE 


5 #. 


L00P3t 


CX1NT2J 


ONE I 
END I 


NOP 

NOP 

HOV VlBASE.VlADDRi SETNA 
HOV VIADDR.VIDEST 
NOV V2BASE.V2ADDR) SETNR 
FNUl. DPX(0),ND> 

WiV LVEC.CNT 
MOV V2ADDR.V2DEST 
FMUIi DPX(0),ND; 

DEC V2DEST 
PMUL; 

DEC VIDEST 

FNUL; DPX(1)<FM 
INC VIADDR) SETNR) 

DPY(1)<FM 

NOP 

INC V2RODR) SETNR 

FNUL DPX(0).MD 

INC VIDEST) NI<DPX(1>) SET»a 

FNUli DPX(0),ND 

DEC CNT 

INC V2DEST) SETNR) MI<DPY(1)) 
BNE LOOP3;FMUI, 

ROD LVEC,V1BRSE 
INC OUTCNT 
SUBt NliORD, OUTCNT 
BEg END) RDD LVEC,V2BRSE 
JMP LOOPl 
SFP 1.0 
RETURN 
SEND 


"WRIT TO ACCESS HD 
"WRIT SOKE MORE 
"GET V1(1,N) 

"LORD TARGET ADDRESS 
"GET V2<1,N) 

"SUN*V1 

"LORD INNER COUNT 
"LORD T.RRGET RDDMTSS 
"SUN*V2 
"LOOP SETUP 
"PUSH 

"LOOP SETUP 
"SAVE V1»SUN 
"GET NEXT VI 
"SAVE V2«SUN 

"GET NEXT V2 
"V1*SIW 

"STORE V1*SUM 

"V2*SUM 

"TEST LOOP COUNT 
"STORE V2*SUM 

"INCREMENT BASE 
"INC COUNTER AND LSDO ADDRESS 
"TEST COUNTER 




58 


FORTRAN IV V02.5-2 Tue 04-NOV-80 10 ; 21 t 09 PAGE 001 

COM 05-14-80 RPIi 66666 66666 66666 77777 22222 66666 

C 

C *••*•********•******»•**»***««»•*•**•••****»********•**»•• 

c 

C This is Program F0R5.F0R which represents a portion of 

C Subroutine NEWX in Processor EIG. Used to obtain timing 

C infojrmation for FORTRAN execution. 

C 

C Corresponding Programs are: 

C APF0R5.F0R - FORTRAN of F0R5. FOR with AP Calls. 

C APNWX3.APH - FP120 Assembler Progr 2 un replacement 

C of FORTRAN portions. 

C APNWX3.ABJ - Object Code of APNWX3.APM. 

C APNWX3.SAV - Linked version of APNWX3.ABJ. 

C 

Q *«**«*•«•*•*«***«*•»*«•«••**«••«****•»****•***••*•**«•••«••• 

C 

0001 DIMENSION tSDO( 3 ),'Z( 3 ),V1( 1000, 3 ) 

0002 DIMENSION ITIM1( 2 ), ITIM2( 2 ) 

0003 NLOADS-3 

0004 LVEC-1000 

0005 DO 1 11-1,3 

0006 LSDOC^I )-II 

0007 1 Z(II)-FL0AT(II) 

0008 DO 2 JJ-1,3 

0009 DO 2 11-1,1000 

0010 2 Vl(II,JJ)-2.0 

0011 WRITE( 6 , 3001 ) 

0012 3001 FORMAT( * INPUT THE NUMBER OP TIMES TO BE EXECUTED* ) 

0013 READ(5,*) NTEST 

0014 CPU1-SECNDS(0. ) 

0015 CALL GTIM(lTIMl) 

0016 DO 270 L-1, NTEST 

0017 DO 1850 N-l,NLOADS 

0018 IF (LSDO(N) .EQ. 0 ) GO TO 1850 

0020 SUM-Z( N ) 

0021 DO 1800 I-1,LVEC 

0022 1800 Vl( I,N)-V1( I,N)*SUM 

0023 1850 CONTINUE 

0024 270 CONTINUE 

0025 CAIiL GTIM(ITIM2) 

0026 CPU2-SECNDS( 0 . ) 

0027 WRITE( 6,50) CPUl 

0028 WRITE(6,50) CPU2 

0029 CALL CVTTIM(ITIM1,IHR,IMI,ISE,ITI) 

OOiO WRITE(6,70) IHR, IMI , ISE, ITI 

0031 CALL CVTTIM(ITIM2,IHR,IMI,ISE,1TI) 

0032 WR1TE(6,70) IHR, IMI , ISE, ITI 
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0034 

0035 

0036 

0037 

0038 

0039 
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50 F0RMAT( • NUMBER OF SECONDS PAST MIDNIGHT - *,F15.7) 

70 FORMAT( • TIME - M2, • J M2, ' » M2, ' : • , 12 ) 

CPU-CPU2-CPU1 
VmiTE(6,22 )CPU 
22 FORMAT( 5X, ' TIME- • , F16 . 9 ) 

STOP 

END 
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COM 05-14-80 RPIt 66666 66666 66666 77777 22222 66666 

C 

C *************************************************** 

c 

C This is Program APF0R5.F0R which contains the AP Calls 

C as replacement for FORTRAN code in FOR5.FOR. Represents 
C a portion of Subroutine NEWX in Processor EIG. Obtains 
C timing information for AP execution, 
c 

C Corresponding Programs are; 

C F0R5.P0R 

C APNWX3.APH 

C APNWX3.ABJ 

C APNWX3.SAV 

C 

^ *«******«***iMkir*******»*«******trir****ft** *********** 

c 

0001 DIMENSION LSD0( 3 }, Z( 3 ) , Vl( 1000, 3 ) 

0002 DIMENSION ITIM1( 2 ), ITIM2( 2 ), ITIM3( 2 ) , ITIM4( 2 ) 

0003 NLOADS-3 

0004 LVEC-1000 

0005 DO 1 11-1,3 

0006 LSDO( II )-II 

0007 1 Z(II)-FLOAT(II) 

0008 DO 2 JJ-1,3 

0009 DO 2 11-1,1000 

0010 2 Vl(II,JJ)-2.0 

0011 VnaTE(6,3001) 

0012 3001 F0RMAT( ’ INPUT THE NUMBER OP TIMES TO BE EXECUTED’ ) 

0013 READ(5,*) NTEST 

0014 CPU1-SECNDS( 0 . ) 

0015 CALL GTIM(ITIMl) 

0016 DO 270 L-1, NTEST 

0017 CALL APCLR 

0018 CALL AP?OT(LSDO(1),0,NLOADS,1) 

0019 CALL APPUT(V1(1,1),10,3000,2) 

0020 CALL APPUT(Z(1),6000,NLOADS,2) 

0021 CALL APWD 

D CALL GTIM(ITIM2) 

0022 CALL APNWX3(10,6000,LVEC,NL0ADS) 

0023 CALL AFWR 

D CALL GTIN(ITIM3) 

0024 CALL APGET(V1(1,1),10,3000,2) 

0025 CALL APWD 

0026 270 CONTINUE 

0027 CALL GTIN(ITIM4) 

0028 CPU2-SECNDS( 0 . ) 

0029 CALL CVTTIM(ITIM1,IHR,IMI,ISE,ITI) 
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0030 

0031 


0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 


WRITE(6,70) 

WRITE(6,75) IHR,IMI,ISE,ITI 
D CALL CVTTIM( ITIM2, IHR, IMI , ISE, ITI ) 

D VmiTE(6,80) 

D VfRITE(6,75) IHR, IMI , ISE, ITI 

D CALLCVTTIM(ITIM3,IHR,IMI,ISE,ITI) 

D WRITE(6,90) 

D WRITE(6,75) IHR, IMI , ISE. ITI 

CALL CVTTIM(ITIM4,IHR,IMI,ISE,ITI) 

WRITE(6,100) 

WRITE(6,75) IHR, IMI, ISE, ITI 
WRITE(6,110) CPUl 
WRITE(6,110) CPU2 
CPU-CPU2-CPU1 
VnaTE(6,120) CPU 

70 FORMAT( • STIME AT START OF DATA INPUT - * ) 

80 FORMAT( * $TIH£ AT COMPLETION OF DATA INPUT AND EXECUTION', 
1 • START - • ) 

90 FORMAT( 'STIME AT END OF EXECUTION AND START OF DATA' 

1 ' OUTPUT - • ) 

100 FORMAT( ' STIME AT COMPLETION OF DATA OUTPUT - ' ) 

75 FORMAT( ' + ',12, ' : ' , 12, ' : ' , 12, ' : ' , 12 ) 

110 FORMAT( ' NUMBER OF SECONDS PAST MIIMJIGHT » ',F15.5) 

120 FORMAT(5X, 'ELAPSED TIME - ',F10.5,' SECONDS') 

C 

C FORTRAN code replaced - Begin 
C 

C DO 1850 N-1,NL0ADS 

C IF (LSDO(N) .EQ. 0) GO TO 1850 

C SUM-Z(N) 

C DO 1800 I-1,LVEC 

C 1800 V1(I,N)-V1(I,N)*SUM 
C 1850 CONTINUE 
C 

C FORTRAN code replaced - End 
C 


0046 

0047 


STOP 

END 


3T1TLE APNWX3 
SENTRY APNWX3,4 


THIS ROUTINE PERFORMS A SECTION OF CODE IN THE NEWX SUBROUTINE 
LOCATED IN THE EIG PROCESSOR 


AUTHOR! K. PERSON 
DATE! MARCH 1979 

Revised: L. J. Feeser and K. Metis 
Date: May 1980 

USAGE 

FORTRAN: CALL APNWX3( V1BASE,ZBASE,LVEC,NL0AD) 
ALL PARAMETERS ARE INTEGERS] 

SPACE 


MAIN DATA MEMORY MAP 


*«*«*«««**««»«««**««**«««« 
« * 

* 2 ARRAY * 

* « 


* * 

* VI ARRAY * 

« * 

n * 

* V2 ARRAY * 

* « 

*«******«***««*«***««*«*«* 


* B ARRAY 
« 


« 

« 

« 




« 


* 

LSDO ARRAY * 

* 


(STARTING ADDRESS) 


NL0ADS+LRN42 *( LV£C*NLOADS ) 


NLOADS4>LRM4-LVEC UNLOADS 


NLOADS*«-LRN 


NLOADS 


0 


ARRAY DESCRIPTIONS 


Z( NLOADS ) 


REAL ARRAY. NEVER PASSED. 


" VI ( LVEC* NLOADS ).. . REAL ARRAY. PASSED ONCE. 
" RETURNED ONCE. 

" V2(LVEC*NLOADS). . . REAL ARRAY. 


" B(LRM> REAL ARRAY. NEVER USED. 

M 

" LSDO( NLOADS) INTEGER ARRAY. PREVIOUSY PASSED. 

•» 


" S-PAD PARAMETERS 


VIBASE 

SEQU 0 

"BASE OF VI ARRAY 


ZBASE 

$EQU 1 

"BASE OF Z ARRAY 


LVEC 

SEQU 2 

"TOTAL NUMBER OF JOINTS * 

DEGREES OF FREEDOM 

NLOAD 

SEQU 3 

"TOTAL NUMBER OF LOADS 


LSDO 

SEQU 4 

"ADDRESS POINTER FOR LSDO 

ARRAY 

OUTCNT 

SEQU 4 

"OUTER LOOP COUNTER 


CNT 

SEQU 5 

"INNER LOOP COUNTER 


VALUE 

SEQU 6 

"TEMPORARY STORAGE 


VIDEST 

SEQU 7 

"VI TARGET ADDRESS 


VIADDR 

SEQU 10 

"VI ADDRESS POINTER 



** FORTRAN: DO 1850 N-1, NLOADS 

" IF(LSDO(N) .EQ. 0) GO TO 1850 

" SUN-Z( N ) 


CLR OUTCNT 

LOOPl: MOV LSDO,LSDO; SETKA 

MOV LVEC.CNT 
ADD# OUTCNT, ZBASE; SETMA 
LDSPI VALUE; DB-f4D 
NOV VALUE, VALUE 
BEQ CONT2 


"GET LSDO(N) 

"LOAD INNER COUNTER 
"GET Z(N) 

"SAVE LSDO(N) 

"TEST LS0b(N)-0 


" DO LOOP CALCULATIONS 

m 

" FORTRAN: DO 1800 I-1,LVEC 

" 1800 V1(I,N)-V1(I,N)*SUM 

" 1850 CONTINUE 


CONTI j 


LOOP2I 


H 

CONT2J 


MOV VIBASE.VIADDR; SETMR) DPX(0)<KD 
MOV VIBASE.VIDEST 
DEC VIDEST 

FMUL DPX(0),MD 
INC VIADDR; SETMX; 

FMUL 

PMULj 

DEC CNT 

INC VIDEST; SETMR; MI<FM; 

BNE LOOP2 


"SAVE SUM IN DPX(O) 
"LOAD ADDRESS REGISTERS 

"V1*SUM 
"GET NEXT VI 
"PUSH 
"PUSH 

"TEST COUNTER 
"SAVE RESULT 


INC OUTCNT 
SUB* NLOAD, OUTCNT 
BEQ END; 

ADD LVEC,V1BASE 

OMP LOOPl 


"INC LSDO ADDRESS 
"TEST OUTER LOOP 

"REAJUST VI BASE 
"GO TO 1300 LOOP 


ENDt 


SEND 


APPENDIX F 


Listings of: 


FOR6.FOR 

APF0R6.F0R 


APRED.APM 
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0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 


COM 05-14*80 RPX« 66666 66666 66666 77777 22222 77777 

C 

c 

C Thia is Progran P0R6.P0R Which raprcaanta a portion of 
C Subroutina RED in Procaasor IMV. uaad to obtain 
C timing information for FORTRAN axacution. 

C 

C Corraaponding Programa arat 

C APP0R6.F0R - FORTRAN Of F0R6.F0R With AP Calla. 

C APRED.APH - FP120 Aaaamblar Progran replacamant 

C of FORTRAN portion. 

C APRED.ABJ - Objact coda of APRED.APM. 

C APRED.SAV - Linked veraion of APRED.ABJ. 

C 

c *••**•**••••*•*••••*•**•*•****•••***•*••••*****•*•*•** 

c 

DIMENSION BB(90,50),MAP(6),5UBMAP(250),S(6,6,50),B(6,15,50) 
DIMENSION ITIM1(2), ITIM2(2) 

INTEGER SUBMAP, CONRNG 
EQUIVALENCE ( B( 1, 1, 1 ) , BB( 1 , 1 ) ) 

NDFCC»I«90 
CONKNG-15 
ZEROD-. 00005 
NDF-6 
NZERO-6 
DO 5 11-1,6 

5 MAP(IX)-II 

DO 6 IX-1,250,10 
SUBMAP(II )-l 
SUBMAP( II41 )-2 
5UBMAP( 1142 )-3 
SUBMAP( 1X43 )-4 
SUBMAP( XI44 }-5 
SUBMAP( 1145 )-6 
SUBMAP{ 1146 )-7 
SUBMAPC 1147 )-e 
SUBMAP( 1X48 )-9 
SUBMAP(X149 )-10 

6 CONTXNUE 
XSTAGE-1 

NDFCON-NDF*CONRNG 
DO 10 XX-1,6 
DO 10 JJ-1,6 
DO 10 nC-1,50 
10 S(XX,JJ,KX)-1. 

WRITE(6,3001) 

3001 F0RMAT( ' INPUT THE NUMBER OF TIMES TO BE EXECUTED' ) 
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0032 READ(5,«) NTEST 

C XF(H2ERO.EQ.O) GO TO 2S00 

C**<* FILL B WITH UNMODIFIED DATA FROM SU3MATR1X LINE 1. 

0033 CPU1-SECNDS( 0 . ) 

0034 CALL GTIM(ITIMI) 

0035 DO 270 LL-1, NTEST 

0036 DO 1000 K-1,NZER0 

0037 M- IABS(MAP(K)) 

0038 DO 1000 J*1,C0NRNC 

003S L- SUBMAP(J) 

0040 DO 1000 1-1. NDF 

0041 1000 B(I,J,K)- S(M,I,L) 

C 

PRELIMINARY B MODIFICATION. 

0042 DO 1400 K-l.NZERO 

0043 M- MAP(K) 

0044 IF(M.LT.O) GO TO 1400 

0046 RA-BB(M.K) 

0047 IF(RA.GT. ZEROD) GO TO 1025 

C NEX-INEX(M) 

0049 ZF(RA.LT. -ZEROD) GO TO 1015 

C KSING-KSING-fl 

C NSING-KSING 

C WRITE(IOUT.IOIO) JOINT, NEX 

ClOlO F0RMAT(49H ••• WARNING. SYSTEM X SINGULAR. JOINT/COMPONENT-15,12) 

0051 RA-.O 

0052 GO TO 1030 

C1015 KNEG-KNEG41 
C NNEG-KNEG 

C IF( IPRT.LT.2 ) GO TO 1025 

C WRITE( IOUT,1020) JOINT. NEX 

C1020 FORMAT( 38HONEGATIVE DIAG TERM. JO I NT/ COMPONENT- 15,12) 

C1025 RA-l./RA 


0053 

1015 

CONTINUE 

0054 

1025 

CONTINUE 

0055 

1030 

BB(H,K)-RA 

0056 


IF(K.EQ.NZERO) GO TO 1200 

0058 


LA-K+1 

0059 


DO 1100 L-LA,NZERO 

0060 


ZA-ZABS(MAP(L)) 

0061 


RAB- RA*BB(ZA,K) 

0062 


DO 1100 I-IA,NDFCCN« 

0063 

1100 

BB(Z.L)- BB(Z.L) -RAB*BB(1,K) 

0064 

1200 

IF(M.EQ.NDF) GO TO 1400 

0066 


INEXT-M+1 

0067 


DO 1300 I-INEXT.NDF 

0068 

1300 

BB(I,K)- BB(I,K)*RA 

0069 

1400 

CONTINUE 
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C 

0070 IF(CONRNG.EQ.l) GO TO 2500 

0072 DO 2100 K-1,NZER0 

0073 N- MAP(K) 

0074 IF(M.LT.O) GO TO 2100 

0076 RA- BB(M,K) 

0077 NSUB- CONRNG 

0078 DO 2000 1-2, CONRNG 

0079 NSUB- NSUB+1 

0080 LS- SUBMAP( NSUB ) 

C 

modify Eli 

0081 DO 1600 IC0L-1,NDF 

0082 RAB- RA*B(ICOL,I,K) 

0083 DO 1600 IROW-l,ICOL 

0084 1600 S(IROW,ICOL,LS)- S( IRON, ICOL,LS ) -RAB*B( IRON, I,K) 

COS CAIJ, CALCS( CONRNG, B, SUBMAP, NDF,S) 

CDC 12 CDS OMITTED 
C 

C *** THE COMPASS ROUTINE CALCS REPLACES THE FOLLOWING ON CDC 
C 

0085 DO 1700 IROW-l,NDF 

0086 1700 B(IROW,I,K)-RA*B(IROW,I,K) 

0087 IF( I. EQ. CONRNG )G0 TO 2000 
C 

c*»»* modify EIJ S 

0089 JA-I+1 

0090 DO 1900 J-JA, CONRNG 

0091 NSUB-NSUB+1 

0092 LS-SUBMAP( NSUB ) 

0093 DO 1800 1C0L-1,NDF 

0094 DO 1800 IROW-l.NDF 

0095 1800 S(IROW,ICOL,LS)-S(IROW,ICOL,LS)-B(IROW,I,K)*B(ICOL,J,K) 

0096 1900 CONTINUE 

0097 2000 CONTINUE 

0098 2100 CONTINUE 

0099 2500 DO 2600 L-1, CONRNG 

0100 K- SUBMAP(L) 

0101 DO 2600 J-1,NDF 

0102 DO 2600 I-1,NDF 

0103 2600 S(I,J,K)-.0 

0104 270 CONTINUE 

0105 CALL GTIM( ITIM2 ) 

0106 CPU2-SECNDS( 0 . ) 

0107 WRITE(6,50) CPUl 

0108 WRITE( 6,50) CPU2 

0109 CALL CVTTIM(IT1M1,1HR,IMI,ISE,ITI) 

0110 WRITE(6,70) IHR, IMI , ISE . ITI 
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0111 

0112 

0113 

0114 

0115 

0116 
0117 

one 

0119 

0120 
0121 


CALL CVTTIM(ITIM2,IHR,IMI,ISE,ITI) 

WRITE(6,70) IHR,IMI.1SE,ITI 

50 PORMAT( • NUMBER OF SECONDS PAST MIDNIGHT - '.F15.7) 
70 FORMAT( • TIME - • , 12, ' : ' , 12 , • : M2, • : • , 12 ) 
CP0-CPU2-CP01 
WRITE(6,33 )CPU 
33 F0RMAT(5X, ‘TIME-* ,F16.8) 

WRITE( 6 , 11 )( BB( II , 1 ) , II-l , 27 ) 

11 FORMAT( F16 . 7 ) 

STOP 

END 
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0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 


COM 05-14-80 RPIt 66666 66666 66666 77777 22222 77777 

C 

c 

C This is Program APFOR6.FOR which contains the AP Calls 

C as replacements for FORTRAN code in FOR6.FOR. Represents 

C a portion of Subroutine RED in Processor INV. Obtains 
C timing information for AP execution. 

C 

C Corresponding Programs are« 

C F0R6.F0R 

C APREO.APM 

C APRED.ABJ 

C APRED.SAV 

C 

c ***************************************************** 

c 

DIMENSION BB( 90,50), MAP( 6 ) , SUBMAP( 250 ) , S( 6 , 6 , 50 ) , B( 6 , 15 , 50 ) 
DIMENSION ITIM1(2), ITIM2( 2 ), IT1M3( 2 ) , ITIM4( 2 ) 

INTEGER SUBMAP,CONRNG 
EQUIVALENCE ( B( 1, 1, 1 ) ,BB{ 1, 1 ) ) 

NDFCON-90 
CONRNG-15 
ZER0D-. 00005 
NDF-6 
NZERO-e 
DO 5 11-1,6 

5 MAP(II)-II 

DO 6 11-1,240,10 
SUBMAP(II)-1 
SUBMAP(II+1 )-2 
SUBMAP(II42 )-3 
SUBMAP(II+3 )-4 
SUBMAP( II+4 )-5 
SUBMAP(II+5 )-6 
SUBMAP(I146 )-7 
SUBMAP( II47 )-8 
SUBHAP( 114-8 )-9 
SUBMAP( 1149 )-10 

6 CONTINUE 
ISTAGE-1 

NDFCON-NDF *CONRNG 
DO 10 11-1,6 
DO 10 JJ-1,6 
DO 10 KK-1,50 
10 S(II,JJ,KK)-1. 

WRITE(6,3001) 

3001 FORMAT( * INPUT THE NUMBER OF TIMES TO BE EXECUTED* ) 
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0032 


0033 

0034 

0035 

0036 

0037 

0038 
0030 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 


0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 

0061 

0062 

0063 

0064 


READ(5,*) NTEST 

C IF(NZERO.EQ.O) GO TO 2500 

C**«* fill B with unmodified DATA FROM SUBMATRIX LINE 1. 
CPU1-SECNDS( 0 . ) 

CALL GTIM(ITINl) 

DO 270 L-1, NTEST 
CALL APCLR 

CALL AFFUT(MAP(1),0,NZERO,1) 

CALL APPUT(SUBMAP( 1), 12, 250,1) 

CALL APPUt,;S( 1,1,1), 262, 1800, 2) 

CALL AFWD 

D CALL GTIM(ITIM2) 

CALL APKED( 2062 ,262,12, NZERO, NDF , CONRNG , ISTAGE , NDFCON ) 
CALL APWR 

D CALL GTIM(ITIM3) 

CALL APGET(B(1,1,1),2062,4500,2) 

CALL APWD 
270 CONTINUE 

CALL GTIM(ITIM4) 

CPU2-SECNDS( 0 . ) 

CALL CVTTIM(ITIM1,IHR,IMI,ISE,ITI) 

WRITE(6,70) 

WRITE(6,75) IHR,IMI,ISE,ITI 
D CALL CVrriM(ITIM2,IHR,IMI,ISE,ITI) 

D WRITE(6,80) 

D WR1TE(6,75) IHR, IMI , ISE, ITI 

D CALLCVTT1M(ITIM3,IHR, IMI,ISE,ITI) 

D WRITE(6,90) 

D WRITE(6,75) IHR, IMI, ISE, ITI 

CALL CVTTIM( ITIM4 , IHR, IMI , ISE , ITI ) 

WRITE(6,100) 

WRITE(6,75) IHR, IMI, ISE, ITI 
WR1TE(6,110) CPUl 
WRITE(6,110) CPU2 
CPU-CPU2-CPU1 
WRITE(6,120) CPU 

70 FORMATC * 5TIME AT START OF DATA INPUT “ • ) 

80 FORMAT( ’ STIME AT COMPLETION OF DATA INPUT AND EXECUTION ’ , 
1 ' START - • ) 

90 FORMAT( 'STIME AT END OF EXECUTION AND START OF DATA' 

1 • OUTPUT - • ) 

100 FORMAT( 'STIME AT COMPLETION OF DATA OUTPUT - ' ) 

75 FORMAT( ' + ' , 12, ' ; ' , 12 , ' < ' , 12 , ' ; ’ , 12 ) 

110 FORMAT( ' NUMBEF. OF SECONDS PAST MIDNIGHT - ',F15.5) 

120 F0RMAT(5X, 'ELAPSED TIME - ',^10.5,' SECONDS') 

C 

C FORTRAN code replaced - Begin 
C 
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C DO 1000 K-1,NZER0 

C M- 1ABS(MAP(K)) 

C DO 1000 J-1,C0NRNG 

C L- SUBMAP(J) 

C DO 1000 I-1,NDF 

C 1000 S(M,I,L) 

CC 

CC**»* PRELIMINARY B MODIFICATION. 

C DO 1400 K«1,NZER0 

C M- MAP(K) 

C IF(M.LT.O) GO TO 1400 

C RA-BB(M,K) 

C IF(RA.GT. ZEROD) GO TO 1025 

CC NEX-INEX(M) 

C IF(RA.LT. -ZEROD) GO TO 1015 

CC KSING-^KSlNG+l 

CC NSING-K.3ING 

CC WRITE(IOOT.IOIO) JOINT, NEX 

CClOlO P0RMAT(49B **• WARNING. SYSTEM K SINGULAR. JOINT/COMPONENT-15,12) 

C RA-.O 

C GO TO 1030 

CC1015 KNEG-KNEG+1 

CC NNEG-KNEG 

CC IF(IPRT.LT.2) GO TO 1025 

CC WR1TE(IOUT,1020) JOINT, NEX 

CC1020 PORMATC 38HONEGATIVE DIAG TERM. JOINT/COMPONENT- 15,12) 

CC1025 RA-l./RA 
C1015 CONTINUE 
C1025 CONTINUE 
C 1030 BB(M,K)-RA 
C IP(K.EQ.NZERO) GO TO 1200 

C LA-K+1 

C DO 1100 L-LA,NZERO 

C IA-IABS(MAP(L)) 

C RAB- RA*BB(IA,K) 

C DO 1100 I-IA,NDFCON 

C 1100 BB(I,L)- BB(I,L) -RAB*BB(I,K) 

C 1200 IP(M.EQ.NDF) GO TO 1400 
C INEXT-M+1 

C DO 1300 I-1NEXT,NDP 

C 1300 iB(I,K)- BB(I,K)*RA 
C 1400 CONTINUE 
CC 

C IP(CONRNG.Ee.l) GO TO 2500 

C DO 2100 K-1,NZER0 

C M- MAP(K) 

C IP(M.LT.O) GO TO 2100 

C RA- BB(M,K) 


73 


FORTRAN IV 


V02.5-2 Tue 04-NOV-80 10 i 23 i 43 


PAGE 004 


0065 

0066 

0067 

0068 


C NSUB- CONRNG 

C DO 2000 1-2, CONRNG 

C NSUB- NSUB41 

C liS- SUBMAP( NSUB ) 

CC 

CC«*»* MODIFY Eli 
C DO 1600 ICOL-l,NDF 

C RAB- RA*B(iCOL,I,K) 

C DO 1600 IROW-l,ICOL 

C 1600 S(IROW,ICOL,LS)- S( IRON, ICOL,LS ) -RAB*B( IRON, I ,K ) 

CCOS CAI.L CA1,CS( CONRNG, B, SUBMAP, NDF,S) 

CCDC 12 CDS OMITTED 
CC 

CC »«» THE COMPASS ROUTINE CALCS REPLACES THE FOLLOWING ON CDC 
CC 

C DO 1700 1R0W-1,NDF 

C 1700 B(IRDW,I,K)-RA*B(IROW,I,K) 

C IF( I. EQ. CONRNG )GO TO 2000 

CC 

CC**** MODIFY EIJ S 
C JA-I+1 

C DO 1900 J-JA, CONRNG 

C NSUB-NSUB+1 

C LS-SUBMAP( NSUB ) 

C DO 1800 1C0L-1,NDF 

C DO 1800 IROW-l,NDF 

C 1800 S(IROW,ICOL,LS)-S(IROW,ICOL,LS)-B(IROW,I,K)*B(ICOL,J,K) 
C 1900 CONTINUE 
C 2000 CONTINUE 
C 2100 CONTINUE 
C 2500 DO 2600 L-1, CONRNG 
C K- SUBMAP(L) 

C DO 2600 J-1,NDF 

C DO 2600 I-1,NDF 

C 2600 S(I,J,K)-.0 
C 

C FORTRAN code replaced - End 

C 

WR1TE( 6 . 200 )( BB( II , 1 ) , II-l , 27 ) 

200 FORMAT(F16.7) 

STOP 

END 


$TITLE APRED 
$ENTRY APRED,8 
«EXT DIV 

THIS ROUTINE PERFORMS THE "RED" SUBROUTINE LOCATED IN THE INV 
MATRIX INVERSION PROCESSOR 

AUTHOR) K. PERSON 

DATE: DECEMBER 1976 

Revised: L. J. Feeser and K. Matis 

Date; May 1980 

USAGE 

FORTRAN: CALL APRED( BBASE, SBASE,MPBASE,NZERO,NDF, CONRNG, ISTAGE.NDFCON ) 

ALL PARAMETERS MUST BE INTEGER] 

DATA PAD STORAGE 

DPX( 2 ) ... CONTAINS NDF 

DPX( 3 ) ... CONTAINS NDF*( K-1 )*CONRNG 

DPY(l) ... CONTAINS NDF»CONRNG 

DPY( 2 ) ... CONTAINS RA 

DPY( 3 ) ... CONTAINS NDF»NDF 


MAIN DATA MEMORY MAP 

«•••««***««*«*••«**•«*«**«••* (STARTING ADDRESS) 
* * 

* » 

* B OR BB ARRAY * 


***«*«*»*****««*«************ 12+LR4+ISIZE*NDP*NDP+LRK 


* 

* 


AK ARRAY 


* 


****«•«****»*«**«*«**»*****•• 124LR4'H'1SIZE*NDP*NDP 
* * 


* S ARRAY 

* 


* 

♦ 


****•******«»»************«**■ 12+LR4 

* <r 




* 


* 




SUBMAP ARRAY 
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m 


M 

n 




« * 

♦ * 

* INEX ARKAY * 

* « 


12 


•* ***»************iMrttK>*«*«***tt« 0 

•• * ft 


ft 

* MAP ARRAY 

ft 


ft 

ft 

ft 


ftftftftftftftftftftftftftftftftftftftftftftftftftftftftft Q 


- ^ARRAY EXPLANATIONS 


n 

** MAP(6) INTEGER ARRAY. TRANSFERED ONCE FOR EACH CAI^L 

- TO THE *APRED* ROUTINE. 

M 

” INEX{6) INTEGER ARRAY. TRANSFERED ONCE FOR 

” EACH CALL TO THE 'APRED* ROUTINE. 


” SUBMAP(LR4) INTERGER ARRAY. TRANSFERED ONLY AFTER 

" BEING READ FROM THE DATA BASE. 

It 

" S(ISIZE*NDF*NDF) ... REAL ARRAY. NEVER TRANSMITTED. 


" AK(LRK) REAL ARRAY. TRANSFERED ONLY AFTER BEING 

” READ FROM THE DATA BASE. 

•t 

” B( NZERO*NDFCON ) REAL ARRAY. RETURNED ONCE FOR EVERY 

" CALL TO THE ROUTINE. 


" S-PAD PARAMETERS 


BBASE 

SEQO 

0 

SBASE 

$EQU 

1 

MPBASE 

$EQU 

2 

NZERO 

SEQU 

3 

NDF 

SEQU 

4 

CONRNG 

SEQU 

5 

NSUB 

SEQU 

6 

SUBMAP 

SEQU 

6 

ISTAGE 

$EQU 

6 

NDFCON 

$EQU 

7 

ICNT 

$EQU 

7 

K 

$EQU 

10 

N27 

SEQU 

11 

BBADDR 

SEQU 

11 


••base ADDRESS OF BB ARRAY 

"BASE ADDRESS OF 5 ARRAY 

"BASE ADDRESS FOR SUBMAP ARRAY 

"NUMBER OF NONZERO SUBMATRICES 

"NUMBER OF DEGREES OF FREEDOM PER JOINT 

"NUMBER OF NONZERO SUBMATRICES IN THIS ROW 

"NUMBER OF SUBMATICES 

"ADDRESS POINTER FOR SUBMAP ARRAY 

"CURRENT JOINT NUMBER 

"NDF*CONRNG 

"LOOP COUNTER 

"LOOP COUNTER 

"CONSTANT-27 . 

"ADDRESS POINTER OP DB ARRAY 


76 


RADDR 

$EQU 

12 

"ADDRESS POINTER 

FOR 

RA-BB(M,K) 

NDFS 

3EQU 

12 

"SCRATCH REGISTER 


J 

$EQU 

12 

"LOOP COUNTER 



BBADRl 

$EQU 

12 

"ADDRESS POINTER 

FOR 

BB ARRAY 

BADDR 

«EQU 

12 

"ADDRESS POINTER 

FOR 

BB ARAY 

BBADRF 

$EQU 

12 

"TARGET ADDRESS FOR BB ARRAY 

SADDR 

SEQU 

13 

"ADDRESS POINTER 

FOR 

S ARRAY 

BBADR2 

5EQU 

13 

"ADDRESS POINTER 

FOR 

BB ARRAY 

BIADDR 

$EQU 

14 

"ADDRESS POINTER 

FOR 

BB ARRAY 

LA 

5EQU 

14 

"LOOP COUNTER 



B2ADDR 

SEQU 

15 

"ADDRESS POINTER 

OF BB ARRAY 

ZERO 

SEQU 

15 

"CONSTANT-O. 



II 

SEQU 

15 

"LOOP COUNTER 



JCNT 

SEQU 

15 

"LOOP COUNTER 



ICOL 

SEQU 

16 

"LOOP COUNTER 



N 

SEQU 

16 

"LOOP COUNTER 



CNT 

SEQU 

16 

"LOOP COUNTER 



lA 

SEQU 

16 

"LOOP COUNTER 



IRON 

SEQU 

17 

"LOOP COUNTER 



COUNT 

SEQU 

17 

"LOOP COUNTER 



M 

SEQU 

17 

"INTEGER-MAP(K) 



INEXT 

SEQU 

17 

"COUNTER 



I 

SEQU 

17 

"LOOP COUNTER 




" FIND NDF*NDF AND STORE THE RESULT IN DPY( 3 ) 

ft 

" FIND NDF»CONRNG AND STORE THE RESULT IN DPY( 1 ) 

H 

" THESE DATA PADS MAY NOT BE USED FOR OTHER PURPOSES 


APRED: LDSPI N27; DB-27. 

MOV CONRNG,CONRNGj DPX(1)<SPFN 
MOV NDF,NDF; DPX(0)<SPFN 
FADD ZER0,MDPX(0); MOV N27,N27 

FADD ZERO,MDPX(l); MOV N27.N27 

DPX(2)<FA; 

FADD 

FMUL DPX(2),FA 

FMUL DPX(2),DPX(2) 
FMUL 

DPY(1)<FM; 

FMUL 

DPY(3)<FM 


"FLOAT CONRNG 
"FLOAT NDF 

"STORE NDF 

"NDF*CONRNG 

"NDF*NDF 

"NDF*CONRNG-DPY( 1 ) 
"SAVE NDF*NDF-DPY{ 3 ) 


" PERFORM THE 1000 LOOP CALCULATIONS 

M 

" FORTRAN j DO 1000 X-1,NZER0 
" M-IABS(MAP(K)) 
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" DO 1000 J-1,C0NRNG 

" L-SUBMAP( J) 

" DO 1000 I-l.NDP 

'• 1000 B(I,J,K)-S(M,I,L) 

•• 

MOV BBASE,BADDR 
DEC BADDR 
CLR K 

LOOPl: MOV K,K; SETMA 

MOV CONRNG.JCNT 
MOV MPBASE, SUBMAP 
DEC SUBMAP 
LDSPI M; DB-MD 
DEC M 


"LOAD ADDRESS OF B( 1,1,1) 
"SET UP FOR LOOPS 

"GET MAP(K) 

"LOAD J LOOP COUNT 
"LOAD SUBMAP BASE 
"GET BASE-1 
"SAVE M-MAP(K) 

"GET M-1 


" PERFORM THE J LOOP BY FINDING THE 
" ADDRESS OF S(M,I,L) 

M 

L00P2i INC SUBMAP > SETMA 
RPSF ONE> DPY(0)<DB 
LDSPI N27; DB-27. 

DPX(0)<MD 

FSUBR DPY(0),MDPX(0); MOV 1 
FADD 

FMUL DPY(3),FA 

FMUL 

FMUL 

DPX(0)<FM 
FIX DPX(O) 

FADD; 

MOV NDF,CNT 

DPX(0)<FA 

LDSPI SADDR; DB-DPX(O) 

ADD SBASE, SADDR 
ADD M, SADDR 
SUB NDF, SADDR 


"GET SUBMAP(J) 

"DPY( 0 )-l . 

"LOAD CONSTANT-27. 
"SAVE L-SUBMAP(J) 

,N27 "FLOAT(L-l) 

"PUSH 

"(L-1)*NDF*NDF 

"PUSH 

"PUSH 

"SAVE RESULT 
"INT( ( L-1 )*NDF*NDF ) 
"PUSH 

"LOAD INNER COUNT 
"SAVE RESULT 

"LOAD ADDRESS POINTER 
"ADD BASE TO POINTER 
"ADD (M-1) TO POINTER 
"LOOP SET-UP 


" PERFORM THE INNER I LOOP 


LOOP 3; NOP 

ADD NDF, SADDR; SETMA 
NOP 

DEC CNT 

ENDS: INC BADDR; SETMA; MI<MD; 

BNE LOOP3 


"GET S(M,I,L) 

"TEST LOOP 
"STORE RESULT 


DEC JCNT 
BEQ CONTI 
JKP LOOP2 

H 

CONTI t INC K 

SUB# NZERO.K 
BEQ CONT2 
JMP DOOPl 


"TEST J LOOP 

"ZP DONE, CONTINUE 

"IP NOT, BRANCH BACK 

"INCREMENT COUNT 
"TEST K LOOP 
"IP DONE, CONTINUE 
"IP NOT, BRANCH BACK 


BEGIN THE 1400 LOOP 


FORTRAN) DO 1400 K-1,NZER0 
M-MAP(K) 

IF (M .LT. 0) GO TO 1400 
RA-BB(M,K) 

IP(RA .CT. ZEROD) GO TO 1025 
NEX-INEX(N) 

IF (RA .LT. -ZEROD) GO TO 1015 


CONT2) MOV BBASE,BBADDR 
CLR K 

L00P4) l«}V K,K; SETHA 
INC K 
NOP 

LDSPI Ml DB-MD 
MOV M,M 
BNE C0NT3; 

DEC M 

JMP CONTI 3 

C0NT3: ADD M,BBADDR; SETMA 

RPSF ZEROD; DPY(0)<DB 
NOP 

DPX(0)<MD 
FSUB DPX(0),DPY(0) 
FADD DPX(0),DPY(0) 
PADD 

BFGT CONT6 
BFGE C0NT4 
JMP CONT5 


"LOAD BB BASE ADDRESS 
"GET M-MAP(K) 


"SAVE M 
"TEST M-0. 

"IP-0, CONTINUE 
"GET M-1 

"IF NOT, GO TO 1400 
"GET BB(M,K) 

"DPY( 0 )-ZEROD 

"RA-DPX( 0 ) 

"RA-ZEROD 

"RA+ZEROD 

"PUSH 

"IF RA>ZEROD, GOTO 1025 
"IP RA>-ZEROD, CONTINUE 
"RA< -ZEROD 


" -ZEROD <RA< ZEROD 

H 

" FOR THIS CASE THE SYSTEM IS SINGULAR 

M 

" FORTRAN: KSING-KSING+1 
" NSING-KSING 


WRITE( lOUT.lOlO) JOINT, HEX 

FORMAT(49H •••WMWING. SYSTEM K SINGUIAR. JOINT/COKPONENT-X5, 12 
RA-0. 

GO TO 1030 


CONT4I DPX(0)<DB; DB-ZERO 
JMP CONT7 


"RA-0. 

"GO TO 1030 


" RA<-ZER0D 

M 

" FOR THIS CASE THE DIAGONAL IS NEGATIVE 

m 

- FORTRAN: XNEG-KNEG4-1 
" NNEG-KNSG 

" 1F(IPRT .LT. 2) GO TO 1025 

" MRITE( IOOT,1020) JOINT, NEX 

** 1020 FORMAT(38HO NONNEGATIVE DIAG TERM. JOINT/COMP-X5, 12 ) 


CONT5 : NOP 


" RA>ZEROD 


" FOR THIS CASE THE DIAGONAL IS O.K. 


" FORTRAN: RA-l.O/RA 
" BB(M,K)-RA 

" IF(K .EQ. NZERO) GO TO 1200 

C0NT6: RPSF ONE; DPY(0)<DB 

NOP 

JSR DIV 

M 

** DIV USES S-PADS 13,14, AND 15 
** THE ANSWER IS RETURNED IN DPX(O) 

M 

CONT7: IK>V BBADDR,BBADDR; SETMA; M1<DPX(0) 

SUBi X, NZERO 
BNE CONT** 

JMP CONTIO 


"0PY(0)-1.0 

"DO DPY(0)/DPX(0) OR 1.0/RA 


"BB(M,K)-RA 
"IS K-0 77 
"IP NOT, CONTINUE 
"IF YES, BRANCH TO 1200 


REAJU5T BB(M,K) ADDRESS TO BB(1,K) 
BY SUBTRACTING M 
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" PERFORH THE 17.00 LOOP CALCULATIONS 
•• 

" FORTRAN ( LA-K-t-X 

" DO 1100 L-LA,N^ltO 

" ZA-IABS(MAP(L)) 

" RAB-RA*BB(IA,K) 

" DO 1100 I-IA,NDPCON 

" 1100 BB(I,L)->BB(I,L)-RAB«BB(I,K) 


CONTSi 


LOOPS: 


M 


SUB H,BBADDR 

MOV BBADDR.BBADRl 

ADD NDFC0N,BBADR1 

M3V K,LA 

MOV LA,LA| SETMA 

NOP 

CLR ZElfO 
LDSPI lA; DB-MD 
MOV IA,IA 
BGE CONT9 
SUB IA,ZERO 
MOV ZERO.IA 


"REAJUST BB(N,K) AIH7RESS 
-LOAD BB(Z,L)-BB(1,K) 
-GET BB(1,L) ADDRESS 
-LA-K 

"GET MAP(L) 


-IA-MAP(L) 

"TEST FOR NEGATIVE lA 
"IF >0, CONTINUE 
"OTHERWISE ZA— lA 


CONT9 j DEC lA 

ADD IA,BBADDR; SETMA 
NOP 

MOV IA,ZI 
FMUL DPX(0),MD; 

MOV BBADDR, BBADDR) SETMA 

FHUL 

FMUL; 

ADD lA.BBADRl; SETMA 

DPX(1)<FM; 

FMUL FM,MD; 

MOV BBADR1.BBADR2 
FMUL; 

DEC BBADR2 


"IA-IA-1 

"GET BB(IA,K) 

"LOAD INNER COUNT 
"RA«BB(ZA,K) 

"GET BB(I,K) 

"PUSH 

"PUSH 

"GET BB(Z,L) 

"SAVE DPX(1)-RAB 
"RAB*BB(Z,K) 

"LOAD BB(I,L) TARGET ADDRESS 
"PUSH 

"LOOP SET-UP 


FMUL 

INC BBADDR; SETMA; 


LOOPGi 


PUSH 

"GET NEXT BB(I,K) 


RI 


FSUBR PM,MD 
FADD; 

INC II 

INC BBAORli SETKA 
FHUL DPX(1),MD) 

SUB# NDFCON.II 
INC BBADR2; SETMA; HI<FA; 

FMUL; 

BNE LOOP6 


"BB(I,L)-RAB*BB(I,K) 

"PUSH 

"INC COUNTER 
"GET NEXT BB(I,E) 
"RAB*BB(7,K) 

"TEST UJOP 
"STORE RESULT 
"PUSH 


SUB NDFCON.BBADDR 
INC LA 

SUB» LA,KZERO 
BEQ CONTIO 
JMP LOOPS 


" FORTRAN* 1200 IF(M .EO- NDF) GOTO 
" INEXT-M+1 

" DO 1300 I-INEXT,NDF 

" 1300 BB(I,K)-BB(I,K)*RA 


CONTIO j INC M 

SUB# M,NDF 
BNE CONTll 
JMP C0NT12 


"REAJUST BB(I,K) ADDRESS 
"INC LA 

"TEST OUTER LOOP 
"IF DONE. CONTINUE 
"IF NOT, BRANCH BACK 


1400 


"TEST M-NDF 

"IF NOT EQUAL, CONTINUE 
"OTHERWISE, BRANCH 


CONTll* ADO INEXT.BBADDR "INEXT-N-t-1 ,G£T ADDRECS 

MOV BBADDR.BBADRF; SETMA "GET BB(INEXT,K) 

DEC BBADRF "LOOP SET-UP 

MOV INEXT.CNT "LOAD COUNTER 


LOOP?* FMUL DPX(0),HD; 

INC CNT 

INC BBADOR; SETMA) 

FMUL 

FMUL) 

SUB# CNT, NDF 

INC BBADRF) SETMA) M1<FN) 

BNE LOOP? 


"RA*BB( I,K) 

"INC COUNTER 
"GET NEXT BB(I,K) 
"PUSH 
"PUSH 

"■"^ST LOOP 
"STORE RESULT 
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" TEST THr. 1400 LOOP 

H 

" FORTRAN: 1400 CONTINUE 

II 

CONT12: SUB NDF.BBADDR 

ADD NDFCON,BBADDR 
SUB# NZERO,K 
BEQ CONTIS 
JMP LOOP4 


"REAJUST BB ADDRESS 
"GET BB(1,K+1) ADDRESS 
"K-NZERO- ? 

"IF SO, CONTINUE 
"IF NOT, BRANCH BACK 


II 


" PERFORM THE 2100 LOOP CALCUALATIONS 


FORTRAN: DO 2100 K=l,NZERO 
M=MAP(K) 

IF(M .LT. O) GO TO 2100 
RA>BB(M,K) 

NSUB=CONRNG 


CONTIS: LDSPI N27; DB-27. 

CLR K 


"CLEAR OUTER COUNT 


" CALCUALATE THE ADDRESS OF BB(M,K) BY FINDING 
" M+(K-l)*NDFCON 

«t 

LOOPS: MOV K,K; SETMA; "GET M«MAP(K) 

DPX( 1 ) <SPFN "STORE ( K-1 ) 

FADD ZERO,MDPX(l);MOV N27,N27 "FLOAT K-1 

FADD 

DPX(0)<MD; "FLOAT M 

•I 

" UPY(l) HAS THE VALUE NDF*CONRNG FROM BEFORE 


C0NT14: 


FMUL DPy(l),FA 

FADD 2ER0,MDPX(0); MOV N27,N27; 
FMUL 


FADD; 

FMUL 

DPY(0)<FA; 

FADD FM,FA 
BFGE CONT14; 

FADD;DPX(S)<FM 

JMP CONT’17 
DPX(0)<FA 
FIX DPX(O) 

FADD; 

MOV CONRNG,NSUB 


"K-1*NDF*C0NRNG 


"M-t-CON' 

"TESV 

"STORK . ONRNG 

"IF <0, BRANCH 

"STORE M+CON*NDF*K-l 


M 


DPX(0)<FA 


LOAD NSUB 
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LDSPI RADDR)DB-DPX(O) 
DEC RADDR 

ADD DBASE, RADDR; SETMA 

N 

" SET UP FOR THE 2000 LOOP 


"STORE NDF«CONRNG*(K-l)+(M-l 
"LOOP SET-UP 

"GET BB(M,K)-RA 


" FORTRAN: DO 2000 I-2.C0NRNG 
" NSUB-NSUB+1 

" LS-SUBMAP( NSUB ) 


CLR ICHT "RESET ICNT 

INC irNT 

DPY(2)<MD "STORE RA IN DPY( 2 ) 

LOOPS: MOV ICNT, ICNT; DPX(0)<SPFN "FLOAT ICNT-1 

FADD ZERO,MDPX(0); MOV N27,N27 

ADD# MPBASE,NSUB; SETMA "GET LS°SUBMAP(NSUB) 

DPY(-1)<DB; DB-40000; WRTMAN; "SET DPY(-1)-1.0 

INC NSUB 

DPY(-1)<DB; DB-1015; WRTEX "SET DPY(-1)-1.0 

DPX(-1)<MD; 

CLR ICOL "CLR COUNTER 


" START ADDRESS CALCULATIONS FOR THE LOOPS 
" TO FOLLOW. FIND (LS-l )*NDF*NDF FOR THE S ARRAY 
" ADDRESS OFFSET. FIND (K-1 )-»CCNRNG»NDF^NDF*( I-l ) 
" FOR THE ADDRESS OFFSET OF THE B ARRAY. 

M 


FSUBR DPY(-1),MDPX(-1); MOV N27,N27 
FADD; 

FMDL DPX(2),FA; 

CLR IROW 

FMUL DPY(3),FA; 

CLR NDFS 

FMUL 

FMUL; 

FADD FM,DPX(3) 

DPX(-1)<FM; 

FADD 

FIX DPX(-l); 

DPY(-3)<FA 
FIX DPY(-3) 

DPX( -1)<FA; 

FADD 

LDSPI SADDR; DB-DPX(-l); 

DPY(-3)<FA 

LDSPI BIADDR; DB-DPY(-3) 

n 

" SET UP ADDRESSES FOR THE 1600 LOOPS 
" AND PERFORM THE 1600 LOOP CALCULATIONS 


"LS -1 
"(I-1)*NDP 
"(LS-1)*NDF*NDF 


"( I-l )*NDF+NDF*CONRNG»( K-1 ) 
"STORE LS-l *NDF»NDF 

"GET INTEGER OF ( LS-l )*NDF*NDF 
"STORE I-l*NDF+NDF»CONRN 

"STORE RESULT 

"STORE OFFSET ADDRESS FOR S 
"STORE RESULT 
"STORE OFFSET FOR B ARRAY 
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" FORTRAN I DO 1600 ICOL-l.NDF 
•• RAB-RA*B(ICOL,I,K) 

•• DO 1600 IR0W-1,IC0L 

" 1600 S( IRON, 1C0L.LS)-S( IIK>W,ICOL,LS )-RAB*B( IRON, 1,K) 

M 

" CALCULATE THE ADDRESS OFFSETS FOR THE S AND B ARRAYS 


ADD SBASE.SADDR "ADD BASE AND OFFSET FOR S 

ADD BBASE,B1ADDR "ADD BASE AND OFFSET FOR B1 

MOV BIADDR, B2ADDR "LOAD B2 ADDRESS 

DEC BIADDR 

LOOPIO: INC BIADDR; SETMA "GET B1(IC0L,I,K) 



INC ICOL 

"INC LOOP COUNT 


SUB IROW,NDFS 

” ge:t ndps-irow 

FMUL DPY(2),MD; 


"MPY RA*B1(IC0L,I,K) 


ADD NDFS, SADDR 

"CHANGE BASE FOR S 

FMUL; 


MOV B2ADDR,B2ADDR; 

SETMA "GET B2(IR0W,I,K) 

FMUL; 


MOV NDF,NDFS 

"LOAD NDFS 

DPY(0)<FM; 


"STORE RAB 

MOV SADDR, SADDR; SETMA 

"GET S( IROW, ICOL, LS) 

FMUL 

DPY(0),MD 

"MPY RAB*B2 

FMUL; 


CLR IROW 

"BESET IROW 


" PERFORM THE INNER 1600 LOOP 


LOOPll! 


ENDlls 


ENDlOi 


FMUL; 


INC IRON 


"INC INNER LOOP 


INC B2ADDR; SETMA; 

FSUBR FM,MD 
FADD 

INC SADDR; SETMA; 

DPX(0)<FA 
FMUL DPY(0),MD; 

SUB# IRON, I COL 


FMUL; 

DECMA; MI<DPX(0); 

BNE LOOPll 


"S-RAB*B2 


"SAVE RESULT IN DATA PAD 
MPY RAB*B2(+1) 

"TEST INNER LOOP 

"STORE RESULTS IN S 


SUB# ICOL.NDF "TEST OUTER LOOP 

SUB ICOL,B2ADDR; BEQ GOT "RESET B2 ADDRESS 
JMP LOOPIO 


" GET READY FOR THE 1700 LOOP 

" B2ADDR HAS THE CORRECT ADDRESS FOR B(iROW.I,K) 

I* 

" FORTRAN! DO 1700 IROW-I.NDF 
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" 1700 B( IR0W,I,K)-RA*B( IR0W,I,K) 

n 

GOT: MOV B2ADDR,B1ADDR; SETMA 

" IROW ALREADY CONTAINS NDF 
” IROW-COUNT (SAME SPAD) 

n 

MOV ICNT,J 
INC J 


"GET FIRST B(IROW,I,K) 


• "ICNT-1 TO J 
"J HAS VTOiUE OF I 


LOOP12: 


FMUL DPY(2),MD 
FMUL; 

INC BIADDR; SETMA 


FMUL; 


DEC COUNT 


"MPY RA*B(IROW,I,K) 
GET NEXT B1 ELEMENT 


END12: MI<FM; DECMA; 

BNE LOOP12 
" STORE RESULT IN B(IROW,I,K) 

•• GET READY FOR THE 1900 LOOP 

•* BASE OF BIADDR IS EQUAL TO THE BASE OF B2(ICOL,J,K) 

" IN THE 1900 LOOP, ALSO B2( ICOL, J,K)-NDF WILL GIVE THE 

•* BASE OF Bl(ICOL,I,K) 

•• 


FORTRAN: IF( I .EQ. CONRNG) GO TO 2000 
JA-I+1 (ALREADY DONE) 

DO 1900 J-JA, CONRNG 
NSUB-NSUB+1 
LS=°SUBMAP( NSUB ) 

MOV BIADDR, B2ADDR 
SOB NDF, BIADDR 
SUB# CONRNG, J 
BNE LOOP13; 

DEC BIADDR 

JMP CONTIS 


"LOAD B2 ADDRESS 
"REAJUST B1 ADDRESS 
"TEST IF ICONRNG 
"IF I IS NOT EQUAL, CONTINUE 

"IF I-CONRNG, GO TO 2000 


LOOP13: ADD# MPBASE , NSUB ; SETMA "GET LS-SUBMAP( NSUB ) 

DPY(-1)<DB; DB=40000; WRTMAN; "DPY(-1)-1.0 

INC NSUB 

DPY(-1)<DB; DB-1015; WRTEX "DPY(1)=»1.0 

DPX( -1)<MD 

FSUBR nPY(-l),MDPX(-l); MOV N27,N27 "LS-1 

FADD 

M 

" CALCUALTE THE OFFSET WITH THE NEW LS VALUE 

M 

FMUL DPY(3),FA "MPY ( LS-1 )»NDF*^NDF 

FMUL 
FMUL 

DPX(-1)<FM; 


STORE THE RESULT 
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MOV NDF.N 
FIX DPX(-1)> 

MOV NDF,M 


•LOAD INNER COUNTER 


•'LOAD OUTER COUNTER 


•• GET READY FOR THE 1800 LOOP CALCUALTIONS 

•t 

•• FORTRAN: DO 1800 ICOL-l.NDF 
•• DO 1800 IROW-l.NDF 

" 1800 S(IROW,ICOL,LS)-S(IROW,ICOL,LS)-B(IROW,I,K)*B(ICOL,J,K) 

H 

•• PERFORM THE ADDRESS SET UP FOR THE LOOPS 


FADD 

DPX(-1)<FA; 

MOV B2ADDR, B2ADDR; SETMA 
LDSPI SADDR; DB-DPX(-l) 

ADD SBASE, SADDR 
DEC SADDR 


-STORE THE S ADRESS 
"GET ADDRESS OF S 
"DEC FOR LOOP 


" PERFORM THE OUTER 1800 LOOP CALCULATION 
•• SINCE B(ICOL,J,K) DOES NOT CHANGE IN THE INNER 
•• LOOP. IT IS ACCESSED ONLY IN THE OUTER LOOP 


LOOP14: DPY(0)<MD 

INC BIADDR; SETMA 
NOP 

INC SADDR; SETMA 
FMUL DPY(0),MD 
FMUL; 

DEC N 


"STORE B2 
"GET B1 
"TIMING NOP 
"GET S ARRAY ELEMENT 
"DO B1*B2 

••DEC INNER COUNT 


PERFORM THE INNER MOST 1800 LOOP CLALCULATIONS 


LOOP15 : FMUL 

FSUBI FM,MD; 

INC BIADDR; SETMA 

FADD 

DPX( -1)<FA; 

INC SADDR; SETMA 

DEC N; 

FMUL DPY(0),MD 

END15 ; MI <DPX( -1 ) ; DECMA; 

BNE LOOP15; 

FMUL 

MOV NDF,N; 

FMUL 

SUB N, BIADDR; 

FSUBR FK,MD 

INC B2A'.m; SETMA; 

FADD 


"DO S-B1*B2 
"GET NEXT B1 ELEMENT 

"STORE S-B1*B2 

"GET NEXT S ELEMEMT 
••TEST INNER LOOP 
"DO (B1+1)*B2 
"STORE S-S-B1*B2 


••RESET INNER LOOP COUNT 

"RESET B1 ADDRESS 
"(S+1)-(B1+1)*B2 
"GET NEXT B2 ELEMENT 


DEC M; 


END14: 


BEQ GO; 


JMP LOOP14 


’•TEST OUTER 1800 LOOP 

DPX(-1)<FA "STORE ( S+1 )-( Bl+1 )*B2 

"SAVE RESULT 

MI<DPX(-1); MOV SADDR,SADDR; SETMA 


" TEST THE OUTER LOOPS FOR COMPLETION 

•t 

" FORTRAN: 1900 CONTINUE 
•• 2000 CONTINUE 

" 2100 CONTINUE 


GO: INC J 

SUB* J.CONRNG 
BEQ C0NT16 
JMP LOOP13 

C0NT16 : INC ICNT 

SUB# ICNT.CONRNG 
BEQ CONTI 7 
JMP LOOP9 

C0NT17: INC K 

SUB* K.NZERO 
BEQ CONTI 8 
JMP LOOP8 

" GET READY FOR THE 2600 LOOP CALCULATIONS 

” FORTRAN: 2500 DO 2600 L«l,CONRNG 

" K-SUBMAP( L ) 

" DO 2600 J-1,NDF 

" DO 2600 I«1,NDF 

" 2600 S(I,J,K)-0. 

RETURN 


"TEST 1900 LOOP 
"TEST IF J-CONRNG 
"IF YES, CONTINUE 
"IF NOT, BRANCH BACK 
"TEST 2000 LOOP 
"IS I-CONRNG? 

"IF YES, CONTINUE 
"IF NOT, BRANCH BACK 
"TEST OUTER 2100 LOOP 
"IS K=N2ERO 
"IF YES, CONTINUE 
"IF NOT, BRANCH BACK 


" FIND THE ADDRESS OFFSET OF S(1,1,K) BY FINDING 
" (K-1)*NDF»NDF 

t* 

COOT18: MOV CONRNG,CNT "LOAD COUNTER 

DEC MPBASE "LOOP SET-UP 

•• 

" DPY(3) HAS TIE VALUE NDF»NDF 

H 

OUTLOP: INC MPBASE; SETMA; 

FIX DPY( 3 ) 

FADD 

DPX(1 )<FA 

DPX(0)<MD 

FSUBR DPY( -1),MDPX(0); MOV N27,N27 
FADD; 


"GETK 

"INT( ^DF*NDF) 


"STORE K 

"K-1 
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LDSPI COUNT; DB-DPX( 1 ) "SET COUNT-NDF*NDF 

FMUL DPY(3),FA; "K-1*NDF*NDF 

DEC COUNT 

FMUL 

FMUL 

DPX(0)<FM "STORE K-1*NDF*NDF 

FIX DPX(O) "FIX RESULT 

FADD 

DPX(0)«FA 

LDSPI SADDR; DB-DPX( 0 ) "LOAD S ADDRESS OFFSET 



ADD 

SBASE, SADDR; SETMA; MI < ZERO 

"CLR FIRST ELEMENT 

INLOP: 

DEC 

COUNT 

"CHECK LOOP 


INC 

SADDR; SETMA; MI < ZERO; BNE INLOP 

"CLR S(I,J,K) 



DEC CNT 

"CHECK OUTER LOOP 



BEQ CONT19 

"IF DONE, GO TO END 

V* 


JMP OUTLOP 

"IF NOT, BRANCH BACK 

ONE: 

$FP 

1.0 


ZEROD: 

SFP 

. 00001 


C0NT19: 

RETURN 



$END 



APPENDIX G 


Listings of 


F0R8.F0R 

APF0R8.F0R 


APAFEX.APM 
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TYPE FOR8 

FORTRAN IV V02.S-2 Tue 04-NOV-80 10 s 24 i 09 PAGE 001 

COM 0!:>-14-80 RPI» 66666 66666 66666 77777 00000 33333 

C 

^ «*ir «*•«***«•*••**««*•«*« «»**«««**«*««*«« ***«•«***« 

c 

C This is Prograr F0R8.F0R which represents a portion of 

C Subroutine AFEX in Processor INV, Used to obtain 

C timing information for FORTRAN execution. 

C 

C Corresponc ing Programs arei 

C APF0R8.F0R - FORTRAN of FOR8.FOR with AP Calls. 

C APAFEX.APM - FP120 Assembler Program replacement 

C of FORTRAN portion. 

C APAFEX.ABJ - Object CiwSe of APAFEX.APM. 

C APAFEX.SAV - Linked version of APAFEX.ABJ. 

c ; 

c 

0001 DItZNSION AK(1000),S(6,6,10),K4(1000) 

0002 DIMENSION ITIM1( 2 ), ITIM2( 2 ) 

0003 INTEGER CONRNG 

0004 DO 10 11-1,1000 

0005 K4( II )-l 

0006 10 AK(II)-1. 

0007 DO 17 11-1,6 

0008 DO 17 JJ-1,6 

0009 DO 17 KK-1,10 

0010 17 S(II,JJ,KK)-FLOAT(II) 

0011 NSUBS-lOO 

0012 LCON^l 

0013 CONRNG-1 

0014 NDF-6 

0015 LK-1 

0016 LCONX-0 

0017 LKSUB=1 

0018 LSUB-1 

0019 WRITE(6,3001) 

0020 3001 F0RMAT( ' INPUT THE NUMBER OF TIMES TO BE EXECUTED* ) 

0021 READ( 5 , • ) N*TEST 
C 

c 

0022 CPD1-SECNDS( 0 . ) 

0023 CALL GTIM(ITIMl) 

0024 DO 270 L-1,NTEST f 

C 

0025 DO 1000 I SUB-1, NSUBS 

0026 IF(ISUB.EQ.l) GO TO 600 

0028 J-IFIX(AK(LK)) 

0029 400 IF( J.EQ.K4(LC0N)) GO TO 600 



FORTRAN IV 


V02 . 5-2 


Tue 04-Nov-BO 10 i 24 t 09 


PAGE 002 


0031 

0032 

0034 

0035 

0036 

0037 
0036 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 
0061 


LCON-LCON+1 

IF(LCON.IiT.LCONX) GO TO 400 
VrRITE(6,500) 

500 FORMAT(29HO«*« MFIIiE/KMAP INCONSISTENCY) 

STOP 

600 LSUB-LCON+CONRNG 
K-K4( LSUB ) 

DO 700 J-1,NDF 
DO 700 I-1,NDF 
S(I,J,K)-S(I,J,K)+AK(LKSUB) 

700 LKSUB-IiKSOB+1 
LCON-LCON+1 
1000 LK-LK+1 
270 CONTINUE 

CALL GTIM(ITIK2) 

CPU2-SECNDS(0. ) 
vmiTE(6,50) CPUl 
WRITE(6,50) CPD2 

CALL CVTTIM( IT1M1,IHR,IMI,ISE,ITI) 

VTRITE(6,70) IHR, IMI , ISE, ITI 
CALL CVTTIM( ITIM2 , IHR, IMI , ISE, ITI ) 

VmiTE(6,70) IHR, IMI, ISE, ITI 

50 FORMAT( ’ NUMBER OF SECONDS PAST MIDNIGHT - ',P15.7) 
70 FORMAT( ’ TIME - ' , 12, * : • , 12 , * : • , 12, ' i ’ , 12 ) 
CPU-CPU2-CPU1 
WRITE( 6 , 12 )CPU, NSUBS 
12 PORMAT( 3X, 'TIME-* ,F16.B,I10) 

D WRITE(6, 2001 )(S(1, 1,1), 1-1,360) 

2001 FORMAT(P17.6) 

STOP 

END 
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FORTRAN IV 


V02 . 5-2 


Tue 04-Nov-eO 10 t 24: 31 


PAGE 001 


COM 

C 

c 

c 

c 

c 

c 

c 

c 


05-14-80 


RPI* 66666 66666 66666 77777 00000 33333 


This is Program APFORB.FOR which contains the AP Calls 
as replacement for FORTRAN code in F0R8.F0R. Represents 
a portion of Subroutine AFEX in Processor INV. obtains 
timing information for AP execution. 



C 

Corresponding Programs are: 


c 

F0R8.F0R 


c 

APAFEX.APM 


c 

APAFEX.ABJ 


c 

APAFEX.SAV 


c 



c 

«««***«««•««*«•••« **«««**«« < 


c 


0001 


DIMENSION AK(1000),S(6,6,10) 

0002 


DIMENSION ITIM1(2),ITIM2(2),: 

0003 


INTEGER CONRNG 

0004 


DO 10 11-1,1000 

0005 


K4(II)-1 

0006 


10 AK(II)-1. 

0007 


DO 17 11-1,6 

0008 


DO 17 JJ-1,6 

0009 


DO 17 KK-1,10 

0010 


17 S(II,JJ,KK)-FLOAT(II) 

0011 


NSUBS-100 

0012 


LCON-l 

0013 


CONRNG-1 

0014 


NDF-6 

0015 


LK-1 

0016 


IiCONX-0 

0017 


LKSUB-1 

0018 


LSUB-1 

0019 


WR1TE( 6 , 3001 ) 

0020 

3001 F0RMAT( * INPUT THE NUMBER OF 

0021 


READ(5,*) NTEST 


C 



c 


0022 


CPU1-SECNDS( 0 . ) 

0023 


CALL GTIM(ITIMl) 

0024 


DO 270 L-1, NTEST 


c 


0025 


CALL APCLR 

0026 


CALL APPUT(K4(1), 12, 1000,1) 

0027 


CALL APPUT(S( 1,1,1), 10x2, 360 

0028 


CALL APIUT(AK(1),1372,1000,2 


«« 9 »«**«*«**«*««*«« 
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FORTRAN IV V02.5-2 Tue 04-Nov-BO 10 i 24 t 31 


PACT 002 


0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 


0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 


CALL APWD 

D CALL GTTM( 1TIM2 ) 

CALL APAF£X( 1012,1372,12, NSUB , UC, LCON , LCONX, NDF , 
1C0NRNG,LKSUB) 

CALL APVm 

D CALL GTIM(1TIH3) 

CAIX APGET(S(1,1,1),1012,360,2) 

CALL APWD 
270 CONTINUE 

CALL GTIM(ITIM4) 

CPU2-SECNDS( 0 . ) 

CALL CVTT1M(ITIM1,IHR,IMI,ISE,IT1> 

WRITE(6,70) 

WRZTE(6,75) IHR, IMI , ISE, ITI 
D CALL CVTTIM(ITIM2,IHR,IMI,ISE,ITI) 

D WRITE(6,80) 

D WRITE(6,75) IHR, IMI , ISE, ITI 

D CALLCVTTIM(ITIM3,IHR,IMI,ISE,ITI) 

D WRITE(6,90) 

D WRITE(6,75) IHR, IMI , ISE, ITI 

CALL CVTTIM(ITIM4,IHR,IMI,ISE,ITI) 

WRITE(6,100) 

WRITE(6,75) IHR, IMI, ISE, ITI 
WRITE(6,110) CPUl 
WRITE(6,110) CPU2 
CP0-CPO2-CPU1 
WRITE(6,120) CPU 

70 F0RMAT{ ’ STIME AT START OF DATA INPUT - • ) 

80 F0RMAT( '3TIME AT COMPLETION OF DATA INPUT AND EXEaniON', 
1 • START - * ) 

90 FORMAT( '5TIME AT END OF EXECUTION AND START OF DATA' 

1 • OUTPUT - ’ ) 

100 FORMAT( • STIME AT COMPLETION OF DATA OUTPUT - ' ) 

75 FORMAT( * + * , 12, ' : ' , 12, ' : ' , 12, ’ i ' , 12 ) 

110 FORMAT( ' NUMBER OF SECONDS PAST MIIWIGHT - ’,F15.5) 

120 F0RMAT(5X, 'ELAPSED TIME - ',F10.5,' SECONDS') 

C DO 1000 1SUB-1,NSUBS 

C IF( ISUB.EQ.l ) GO TO 600 

C J-IFIX(AK{LK)) 

C 400 IF( J.EQ.K4(LC0N)) GO TO 600 
C LCON-LCOM+1 

C IF( LCON. LT. LCONX) GO 'PO 400 

C WRI'TE(6,500) 

C 500 FORMAT(29HO*** MFILE/KMAP INCONSISTENCY) 

C STOP 

C 600 LSUB-LCON-KTONRNG 
C K-K4( LSUB ) 

C DO 700 J-1,NDF 
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FORTKAN IV V02 . S>2 Tue 04-NOV-80 10 s 24 i 31 PAGE 003 



C 

DO 700 I-1,NDF 


c 

5( I , J,K )-S( I , J,K HAK( UCSUB ) 


c 

700 LKSUB-UCSOB+l 


c 

LC0N-U;0N41 


c 

1000 LK-LK41 

0054 


WR1TE( 6, 2001 )( S( 1. 1, 1 ), 1-1, 360) 

0055 

2001 PORMAT(F17.6) 

0056 


STOP 

0057 


END 


$TITLE APAFEX 
SENTRY APAFEX.IO 

ft 

" THIS ROUTINE SIMULATES A SECTION OF CODE IN THE hFZX SUBROUTINE 
" THAT IS LOCATED IN THE MATRIX INVERSION ROUTINE INV 

H 

" AUTHORi X. PERSON 
" DATE. APRIL 79 

** Revised: L. J. Feeser and X. Metis 
” Date: May 1980 


" “USUAGE 

- FORTRAN: CALL APAFEX( SBASE, AKBASE,K4BASE,NSUB,LK,LC0N,LC0NX,NDP, . 
" CONRNG,UCSUB) 

" ALL PARAMETERS MUST BE INTEGER VALUES] 


" MAIN DATA MEMORY MAP 


(STARTING ADDRESS) 


12+LR4+ISIZE*NDF*NUF+LRK 


AX ARRAY 


12+LR4+ISIZE*NDF*NDF 


S ARRAY 


12tLR4 


X4 ARRAY 


12 


SI 


6 
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" ARRAY EXPLANATIONS 


" AK(LRK) REAL ARRAY. TRANSFERS!) ONLY AFTER 

” BEING READ FROM THE DATA BASE. * 

M 

" S(ISIZE*NDF*NDF) REAL ARRAY. NEVER TRANSFERED. 

•« 

" K4(LR4 ) INTEGER ARRAY. TRANSFERED (MILY AFTER 

" BEING READ FROM THE DATA BASE. 


" S-PAD PARAMATERS 


SBASE 

$EQU 

0 

"BASE ADDRESS OP S(1,J,K) 

AKBASE 

SEQO 

1 

"BASE ADDRESS OF AK(LK) ARRAY 

K4BASE 

$Egu 

2 r 

••base ADDRESS OF K4(LC0N) ARRAY . 

NSUB 

SEQU 

3 

"OUTER LOOP COUNT 

LK 

$EQU 

4 - 

"ADDRESS POINTER FOR AK ARRAY “ 

LCON 

$EQO 

5 

"ADDRESS POINTER FOR K4 ARRAY 

LCONX 

$EQU 

6 _ 

- "CONSTANT-L4+CONRNG 

NDF 

$EQO 

7 

"NUMBER OF DEGREES OP FREEDOM 

NDF2 

SEQO 

7 

"NDF*NDP 

CONRNG 

$EQU 

10 

"NUMBER OF NONZERO SUBMATRICES 

LKSOB 

SEQO 

11 

"ADDRESS POINTER FOR AK ARRAY 

TEMP 

SEQU 

11 

"SCRATCH REGISTER-K4(LCON) 

J 

$EQU 

12 

"SCRATCH REGISTER-IPIX( AK( LK ) ) 

SIADDR 

SEQO 

12 

"ADDRESS POINTER FOR S(I,J,K) 

S2ADDR 

SEQU 

13 

"ADDRESS POINTER FOR TARGET S(I,J,K) 

AKADDR 

SEQO 

14 

"ADDRESS POINTER FOR AK ARRAY 

CNT 

SEQU 

15 

"INNER LOOP COUNT 

N27 

SEQU 

15 

"CONSTANT-27 . 

ISUB 

SEQO 

16 

"OUTER LOOP COUNTER 

LSUB 

SEQO 

17 

"ADDRESS POINTER FOR K4 ARRAY 


" CALCULATE NDF*NDF AND STORE THE RESULT 

" THIS VALUE IS USED TO CALCULATE THE ADDRESS OF S(NDF,NDF,K) 
" FOR ANY GIVEN K. 


APAFEX: LDSPI N27; DB-27. 

MOV NDF,NDF; DPX(0)<SPFN 

FADD ZERO,MDPX(0); MOV N27,N27 

FADD 

DPX(0)<FA 

FMUL DPX(0),DPX(0) 

C - 


"LOAD CONSTANT-27. 
"DPX( 0 )-NDF 
"FLOAT NDF 
"POSH 

"SAVE FLOAT(NDF) 
"NDF*NDF 


FNUL; 


CLR ISUB 


PMUL; 

MOV AKBi^E,AKM>DR 

DPX(2)<FM; 

ADD LKSOB,AKADDR 

FIX DPX(2) 

FADD 

DPX(1)<FA 

LDSPI NDF2; DB=DPX(1) 

JMP LOOP3 


"PUSH 

"CLEAR OUTER COtRlT 
"POSH 

"LOAD ADDRESS POINTER 
"STORE NDP»NDF IN DPX(O) 
"POINT TO AK(LKSUB) 
"FIX(NDF*NDF> 

"PUSH 

"SAVE NDF*NDF 

"NDF2“NDF*ftDF 


" FORTRAN: DO 1000 ISUB-1,NSUB 

" IF (ISUB .EQ. 1) GO TO 600 

" J“IFIX(AK(UC)) 

" 400 IF (J .EQ. K4(LCON)) GO TO 600 


LOOPl: ADD* LK,AKBASE; SETMA 

NOP 

ADD* LC0N,K4BASE; SETMA 

FIX MD 
FADD 

DPX(1)<FA 

LDSPI J; DB=DPX(1) 

LDSPI TEMP; DB>MD 
SUB* J,TEMP 
BEQ LOOPS 

n 

" FORTRAN: LCON=LCON+l 

" IF (LCON .LT. LCONX) GO TO 400 


"GET AK(LK) 

"GET K4(LCON> 
"INT(AK(LK)) 

"PUSH 
"SAVE J 

"STORE J=AK(LK) 
"STORE K4(LCON) 
"J»K4(LCON) ? 
"BRANCH IF YES 


LOOP2 : INC LCON 

SUB* LCON, LCONX 
BGT CONT 
OMP ERROR 


"LCON .LT. LCONX ? 
"IF YES, CONTINUE 
"IF NOT, ERROR 


CONT: ADD# LC0N,K4BASE; SETMA "GET K4(LC0N+1) 

NOP 
NOP 

LDSPI TEMP; DB-MD 

SUB* J,TEMP "K4(LC0N+1)-J ? 

BEQ LOOPS "IF YES, BRANCH 


JMP LOOP2 


CONTINUE TEST 


FIND THE STARTING ADDRESS OF THE S ARRAY 
GIVEN THAT DPX( 2 )-NDF*NDF 


" FORTRAN 1 K-K4(LSUB) 


LOOP3; NOV K4BASE,LSUB 
ADD LCON.I^UB 
ADD CONRNG.LSUB; SETHA 
RPSF ONE; DPY(0)<DB 
IDSPI N27; DB-27 . 
DPX(0)<MD 

FSUBR DPY(0),MDPX(0); MOV 
FADD 

FMUL DPX(2),FA 

FMUL 

FNUL; 

INC JSUB 

DPX(0)<FM; 

INC LCON 

FIX DPX(O) 

FADD 

DPX(1)<KA 

LDSPI SIADDR; DB-DPX( 1 ) 


"FIND LSUB 

"GET K4(LSUB) 

"DPY(0)-1. 

"LOAD CONSTANT 
"SAVE K 

, N27 "FLOAT K-1 

"PUSH 

"DO K-1 * NDF**2 
"PUSH 
"PUSH 

"INC OUTER LOOP COUNT 
"SAVE (K-1>*NDF**2 
"INC LCON 

"INT( ( K-1 )*NDF*»2 ) 

"PUSH 

"SAVE ADDRESS POINTER 

"SAVE ADDRESS POINTER 


" PERFORM TOE ADDITION OF AK AND S 


" FORTRAN: DO 700 J-1,NDF 

" DO 700 I-l.NDF 

" S(I,J,K)-S(I,a.K)+AK(LKSUB) 

" 700 LKSUB-LRSUB41 



ADD SBASE.SIADOR; SETMA 
MOV SIADDR, S2ADDR 
NOV AKADDR,AKADDR; SETNA 
DPY(0)<ND; 

MOV NDF2,CNT 

DEC S2ADDR 

L00P4: FADD DPX(0),ND 

INC SIADDR; SETHA; 

FADD 

NOP 

INC AKADDR; SETMA 
DPX(0)<ND; 

DEC CNT 


"GET S(1,1,K) 
"LOAD TARGET ADDRESS 
"GET AK(LKSUB) 

"SAVE S(1,1,K) 

"LOAD COUNTER 
"LOOP "fT-UP 
"S(I, J,K)+AK(LKSUB) 
"GET NEXT S(I,J,K) 
"PUSH 

"GET NEXT AK(LKSUB) 
"SAVE S(I,J,K) 

"TEST COUNT 
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INC S2MDDR; SETNAl NI<FAj 
BNC UX)P4 


"SAVE RESULT 



SUB» NSUB.ISUB 

"TEST OUTER HOST LOOP 


BEQ CONTI 

• 


JKP LOOPl 

"START AGAIN 

ONEi 

SFP 1.0 


ERROR 1 

NOP 


CONTI « 

RETURN 

SEND 



APPENDIX H 


Listings of: 

F0L9.F0R 

APF0R9.F0R 


APMLTX.APM 


101 


.TYPE FOR9 

FORTRAN IV V02.5-2 Tue 04-Nov-BO 10 j 25 j 36 PACE 001 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 
0017 
0v018 

0019 

0020 
0021 


COM 05>14-80 RPI# 66666 66666 66666 77777 22222 

CCOMDECK MULTEX 
C 

Q «t 

C 

C This is Program FOR9.FOR which represents a portion of 
C Subroutine MULTEX in COEC Library. Used to c^tain 
C timing information for FORTRAN execution. 

C 

C Corresponding Programs are: 

C APF0R9.F0R - FORTRAN of FOR9.FOR with AP Calls. 

C APMLTX.APM - FP120 Assembler Program replacement 

C of FORTRAN x>ortion. Called from APSMLT. 

C APMLTX.ABJ ~ Object code of APML'nC.APM. 

C APMLTX.SAV - 

C 

C FOR7.FOR 

C APF0R7.F0R 

C APSMLT. APM 

C APSMLT. ABJ 

C APSMLT. SAV 

C 

C *««••«•••*«•«•*«•*««•««*•••*«•«*««•*»««•*«««**« 

C 

C SUBROUTINE MULTEX 

C S{ N.NSUBS, JLIST,A,VIN,VOUT) 

C 8/74 WD WHETSTONE 

DIMENSION A(6,6,50),V1W(6.50),V0UT(6,50), JLIST(50) 
DIMENSION ITIM1(2).ITIN2(2) 

NSUBS-50 

N-6 

C 

DO 11 11-1,50,5 
JLIST( II )-l 
JLIST(II+1 )-2 
JLIST( 11+2 )-3 
JLIST(II+3)-4 
JLIST(II+4)-5 
11 CONTINUE 

DO 12 11-1,6 
DO 12 JJ-1,6 
DO 12 KK-1,50 
IF(II.NE.l) GO TO 12 
VIN( JJ , KK )-FL0AT( J J ) 

12 A(II,JJ,KX)-FLOAT(II) 

WRITE( 6 , 3001 ) 

3001 PORMAT( • INPUT THE NUMBER OF TIMES TO BE EXECUTED* ) 
READ(5,*) NTEST 
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0022 

0023 

0024 

0025 

0026 

0027 

0028 
0029 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 


CPU1-SECNDS(0. ) 

CALL GTIN(XTIMI) 

DO 270 LL-1,NTEST 
I-JLIST( 1 ) 

DO 100 L-1,N 
DO 100 N-1,N 

100 VOUT( L, I )-VOUT( L, I )+A( L, M, 1 )*VIH< M, 1 ) 

IF( NSUBS.lt. 2) GO TO 300 
DO 200 K-2,NSUBS 
J-JLIST(K) 

DO 200 L-1,N 
DO 200 N-1,N 

VOUT( L, I )-V00T( L, I )+A( L, M, K ) *VIN{ M, J ) 

200 VOUT( L, J )-VOUT( L, J )+A( M, L,K )*VIN( M, I ) 

270 CONTINUE 

CALL GTIH(ITxM2) 

CPU2->SECNDS(0. ) 

WRITE(6,50) CPUl 
HRITE(6,50) CPU2 

CALL CVTTIM(ITIM1,IHR,IMI,ISE,ITI) 
WRITE(6,70) IHR,IHI,ISE,ITI 
CALL CVTTIM(ITIN2,IHR,INI,ISE,ITI) 
VnaTE(6,70) IHR,IMI,ISE,ITI 
50 FORMAT( ' NUMBER OF SECONDS PAST MIDNIGHT 
70 FORMAT( ' TIME - * ,12, ’ J ’ ,12, • i M2, ' j M2 
CPU-CPD2-CPU1 
WRITE(6,33)CPD 
33 FORMAT(SX,5HTIME-,F16.7) 

STOP 

300 STOP 'NSUBS.lt. 2' 

END 
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COM 05-14-80 RPIi 66666 66666 66666 77777 22222 33333 

CCOMDECX MULTEX 
C 

C ****•****•«***••****«.«*****••******«********•** 

c 

C This is Program APFOR9.FOR which contains the AP Calls 

C as replacements for FORTRAN code in FOR9.F0R. Represents 

C a portion of Subroutine MULTEX in CDEC LibraLry. Obtains 
C timing information for AP execution, 

c 

C Corresponding Programs aret 

C FOR9.FOR 

C APMLTX.APM 

C APMLTX.ABJ 

C APMLTX.SAV 

C 

c ************************************************ 

c 

C SUBROUTINE MULTEX 

C S( N,NSUBS,JLIST,A,VXN,VOUT) 

C 8/74 Vn> WHETSTONE 

0001 DIMENSION A( 6,6,50 ), VIN( 6, 50 ) , VOUT( 6, SO ), JLIST( 50 ) 

0002 DIMENSION ITIM1( 2 ), ITIM2( 2 ) , ITIM3( 2 ) , ITIM4( 2 ) 

0003 NSUBS-50 

0004 N-6 
C 

0005 DO 11 11-1,46,5 

0006 JLIST( II )-l 

0007 0LIST( II+l )-2 

0008 JLIST( I 1+2 )-3 

0009 JLIST( II+3 )-4 

0010 JLIST( II+4 )-5 

0011 11 CONTINUE 

0012 DO 12 11-1,6 

0013 DO 12 JJ-1,6 

0014 DO 12 XX-1,50 

0015 IF(II.NE.l) GO TO 12 

i 0017 VIN( JJ,KK)-FLOAT( JJ> 

j 0018 12 A(II,JJ,KK)-FL0AT(1I) 

I 0019 WRITE(6,3001) 

0020 3001 F0RMAT( ' INPUT THE NUMBER OP TIMES TO BE EXECUTED' ) 

0021 READ( 5 , * )NTEST 

0022 CPU1-SECNDS(0. ) 

0023 CALL GTIM(ITIMl) 

0024 DO 270 LL-l.NTEST 

0025 CALL APCLR 

0026 CALL APPUT(JLIST(1), 0,50,1) 

0027 CALL APPUT(A(1,1,1),50,1800,2) 
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0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 
004.\ 


0042 

0043 

0044 

0045 
004«> 

0047 

0048 

0049 

0050 

005 \ 

0052 

0053 

0054 

0055 


CALL APPUT(VIN(1,1),1850.300,2) 

CALL APWD 

D CALL GTXN(IT1N2) 

CALL APNLTX( 50, 2150. 1850, 0,N,NSUBS) 

CAU. APWR 

D CALL GTIN(ITXH3) 

CALL APGET(VOUT(1,1),2150,300,2) 

CALL APMD 
270 CONTINUE 

CALL GTXM(XTXN4) 

CPU2-SECNDS(0. ) 

CALL GTXN(XTXN4) 

CPU2-SECNDS(0. ) 

CALL CVTTXN( XTXMl , XHR, INI . XSE , XTX ) 

WRITE(6,70) 

WR1TE(6,75) XRR, INI, XSE, XTX 
D CALL CVTTXN(XTXN2,XHR,XICX,ISE,XTX) 

D WRI!TE(6,80) 

D WRITE(6,75) XHR, INI , XSE, XTX 

D CALLCVTTTM( XTXM3 , XHR, IMX , XSE , XTX ) 

D WRXT£(6,90) 

D WRITE(6,75) XHR, XNX , XSE. XTX 

CALL CVTTIN( ITIM4, XHR, XMI , XSE, XTX ) 

WRITE(6,100) 

NRJTE(6,75) XRR, INI , XSE, XTX 
WR1TE(6,110) CPUl 
WRITE(6,110) CPU2 
CP0-CPO2-CPU1 
WRITE(6,120) CPU 

70 FORI«T{ • STINE AT START OF DATA INPUT - * ) 

80 FORNAT( 'STINE AT CONPLETXON OF DATA INPUT AND EXECUTKNf', 
1 • START - • ) 

90 FORNAT( 'STINE AT END OF EXECUTION AND START OF DATA' 

1 ' OUTPUT - • ) 

100 PORNAT( ' STINE AT CONPLETXON OF DATA OUTPUT - ' ) 

75 FORNAT( , X2, ' i ' , X2, ' j ' , X2 , ' j ' , X2 ) 

110 FORNAT( ' NUMBER OF SECONDS PAST NIDNXGHT - ',F15.5) 

120 F0RMAT(5X, 'ELAPSED TINE - ',F10.5,' SEC(»tDS' ) 

C I-JLIST(l) 

C DO 100 L*1,N 

C DO 100 N-1,N 

C 100 VOOT(L,X)-VOOT(L,X)+A(L,N,1>*VIN(N,X) 

C IF(NSUBS.LT.2) 00 TO 300 

C DO 200 K-2,NSUBS 

C J-JLIST(K) 

C DO 200 L-1,M 

C DO 200 N-1,N 

C VOUT( L, I ;-VOUT( L, I )+A( L, N, K )*VXN( N, J ) 
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C 200 V00T(I,,J)-V00T<L,J)+A(M,L,K)*V1N(M,I> 

0056 STOP 

0057 300 STOP 'NSUBS.LT.2' 

0056 END 


9TITIX AMILTX 
GENTRY APHLTX.e 

SINUIATES SUBROUTINE MULTEX(N,NSUBS, JLXST,A,VXN,VOUT) 

AUTHOR! K. FERSCM 
DATE: FEB 79 

Raviatdi L. J. PMMr And K. Itttia 
Date I Nay 1980 


USAGZ 

FORTRAN! CALL APNLnC( ABASE, INBASE, OTBASE.JBASE,NDF,NSUB) 
ALL PARAMETERS MUST BE INTEOCR] 

^MAIN DATA HENORY MAP 


IN ARRAY 


OUT ARRAY 


A OR lA ARRAY 


LSDO ARRAY 


(STARTING ADDRESS) 


NV’iA+NV*LVOUT 


NV+LR 


NV 


0 


ARRAY EXPLANATICMIS 

IN(NV*LVXN) REAL ARRAY. 

OUT(NV»LVOOT). . . REAI. ARRAY. 
lA(LR) REAL ARRAY. 


LSDO(NV) 


INTEGER ARRAY 


S-PM) PARAMETERS 


m 


A 

SEQU 

0 

"ADDRESS OP A OR A(L,N,K) 

ABASE 

SEQU 

0 

"BASE ADDRESS OF ARRAY A 

INBASE SEQO 

1 

"BASE AZX)RES3 OF ARRAY VIN 

OTBASE 

SEQU 

2 

"BASE A1»RE8S OF ARRAY VOUT 

JBASE 

SEQU 

3 

"BASE AI»R£SS OF ARRAY JLZ8T 

NDF 

SE{^ 

4 

"NDF-N NUMBER OF DEG11EES OF 

NSUB 

SEQU 

5 

"NUMBER OP SUBMATRICES 

VOUT 

SE(^ 

6 

"ADDRESS OF VOUT 

C 

SE(^ 

f 

"ADDRESS OF C OR VCXn'( L, X ) 

VIN 

SE(^ 

7 

"ADDRESS OF VIN 

F 

SE{^ 

7 

"AIX>RESS OF F OR VOUT(L,J) 

NDPSQ 

SEQU 

10 

"NDF SQUARED 

M27 

SEQU 

11 

"CCmSTANT 

OUTCT 

SEQU 

12 

"LOOP COUNTER 

CNT 

SEQU 

12 

"LOOP COUNTS 

COUNT 

SEQO 

13 

"LOOP COUNTER 

CT2 

SEQU 

13 

"LOOP COUNTER 

D 

SEQU 

14 

"AI»)RESS OF D OR A(M,L,K) 

E 

SEQU 

15 

"ADDRESS OF E OR VZN(N,Z) 

B 

SEQU 

16 

"A1X)RESS OF B OR VIN(N,J) 

CNTER 

SEQU 

17 

"LOOP COUNTER 


"CALCULAn: THE STARTING EI£MENT OP THE VOUT AND VXN ARRAYS 
" THEREFORE FIND (I-1)*NDF 

N 

" FORTRAN! I-JLZST(l) 


APNLTXt LDSPZ N27) OB-27. 

DPY(l)<DBt DB-40000) NRTMAM 
DPY(1)<DB; DB-1015) WRTEX 
MOV JBASE.JBASEi SETMA 

MOV NDF.NDF; DPX(2)<SPFM 
FAIH) ZERO,NDPX(2); NOV N27.N27 

DPX(0)<MD 

FSUBR DPY(1).MDPX(0)) NOV N27.N27 
FADD) 

OPX(2)<FA 
FMITI. DPX(2),FA 

FNUL DPX(2),DPX(2) 

FMJL 

0PX(0)<FN; 

FNDL 

FIX DPX(O); 

m>A(0)<FM 
FIX DPX(0) 

DPX(0)<FA 

U>SPZ VIM; DB-DPX(0)} 


"LOAD CONSTANT 
"DPY( 1 )-l 

"GET Z-JLZST(l) 
"FLOAT NOF 

"STORE X 
"I-l 


"STORE NDF 
"(I-1)*NOT 
"NDF'NW 

"STORE I-l ‘NDF 

"INT( I-1*NDF) 
"STORE HDF»mF 

"STORE( I-l) "NDF 
"LOAD VIN 
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PAOD 

DPX(0)<PA "STORE NOP'KDP 

U)SPX NDPSQ; DB-DPX(O) "STORE NDP*NDF 

" SET UP FOR 100 LOOP 


" FOKTRAHi DO 100 L-1,N 

" DO 100 M-l,» 

" 100 VOOT(L,I)-VOUT(L,I)+A(L,M,l)«VlN(M,I) 

•• 

MOV KDF.OUTCr 
MOV VIH.VOUT 
DEC VOOT 
ADD INBASE, VIN 
ADD OTBASE.VOUT 
DEC ABASE 
KDD NDF,VIN 

N 

"DO V0UT(L,I)-V00T(L,I)+A(L,M,X)*V1N(M,J) 

ft 

LOOPli INC VOUT;SE01« 

NOP 

INC ABASE ;SETHA 

FAIX3 ZERO,ND 

SUB NDF.VIN; SETMA; FADD 

DPX(1)<MD "STORE A(L,M,I) 

MOV NDF, COUNT 

L00P2I FMUL DPX(l),hO 

AIX> NDF, ABASE; SETMA; 

FMUL 

FMUL 

INC VIN; SETHS; 

FA',M) 2M.FA 

DPX(1)<MD; 

FAIN); 

DEC COUNT 

END2i MI<FA; NOV VOUT,VDUT; SETMA; 

BNE L00P2 
DEC OUTCT 
BEQ CONTI 

SUB NDFSQ, ABASE 

ENDl I JNP LOOPl 

M 

"CHECK FOR NSUB LESS THAN TWO (IE..NSUB-1) 

« 

" FORTRAN: IF(NSUBS .LT. 2 ) 00 TO 300 


"LOAD COUNTER 
"LOAD VOUT 

"GET VIN ADDRESS 
"GET VOUT ADDRESS 


"GET VOUT(L.I) 

"GET A(L,H,I) 
"VOUT(L,I>+0 
"GET VIN(M,I) 

"LOAD COUNT 
"VIN«A 

"CTT A(L,M-H,I ) 


"G^T VIN(M+1, 1 ) 
"VOUT+VIN*A 
"STORE A(L,M+1,I) 


"STORE VOUT 
"TEST M-NDF 

"TEST L-NDF 
"REAJUST A ADDRESS 


corn t DEC§ NSUB 
WE C(Nrr2; 
DEC NDFSQ 


GET NDF**2-1 
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JMP END 

W 

" SET DP FOR THE 200 U)OPS 


" FORTRAN! DO 200 K-2,NSUBS 

" J-JLIST(K) 

** DO 200 L-1,N 

" DO 200 1^1, N 

- VODT( L, I )-VOOT( L, I )+A( L,M,K )*VTN(M, J ) 

" 200 VOOT(L,J)-VDUT(L,J)+A(M,L,K)*VIN(M,l) 

M 

C0NT2: SUB NDF, ABASE 

MOV ABASE, D 
SOB NDF, ABASE 
IiDSPI CarPER; DB-1. 

INC A 
INC C 
NOV VIN,E 
SUB NDF,E 
DEC E 

nav NDF,CT2 


"ADD OP A(N,L,K) 
"REAJUST A 
"LOAD CNTER-1 


"LOAD VIN(M,I) 

"REAJUST C 
"LOAD COUNTER 


"CALCULATE THE STARTING ELEMENTS OF VOUT(L,J) AND VIN(M,J) 
" THEREFORE (J-1)*NDP 


LOOPS: INC JBASE; SETMA 

INC CNTER 

NOP 

DPX(0)<MD; SUB NDF,C 

FSUBR DPY(1),MDPX(0>; MOV N27,N27 

FADD; 

MOV NDF,CNT 

FMUL DPX(2),FA 

FMOL 

FMUL 

DPX(0)<FM 
FIX DPX(O) 

FADD 

DPX(0)<FA 

LDSPI F; DB-DPX(O) 

MOV F,B 
ADD OTBASE,P 
ADD INBASE, B 
DEC B 


"GET J 

"INCREMENT K COUNT 

"STORE J, INC VDD 
" J-1 


"LOAD COUNT 
"J-1*NDF 


"STORE J-1*NDP 


"STORE INCREMENT IN F 
"STORE INCREMENT IN B 
"GET F ADDRESS 
"GET B ADDRESS 


"DO ... VOUT(L,I)-VOUT(L,I)+A(L,M,K)*VIN(M,J) 
" VOUT(L,J)-VOUT(L,J)+A(M,L,K)*VIN(M,I) 


L00P4) 


GET K 


END4i 

CX)NT3j 


END3i 


A13D NDP.ABUSEi SETKE 
NOP 

INC B) smov 

DPX(1)<ND 

NOV C.C) SETNE 

FNUL DPX(1),ND 
INC D| SETNX) 

FNUL 

DPY(2)<ND; 

niuL 

INC E) SETNX) 

PADD rM,DPY(2) 

DPX(0)<ND) 

FADD 


NOV C.CtSETNAt HltFA 
FNUL DPX(0),ND 
NOV F.F;SETNA) 


FNUr, 

FNUL 

NOP 

FAOD FM,ND 
FADDi 

DEC CNT 

NI<FA« NOV P,F> SETMA; 

BEQ CONT3 
JNP LOOP4 


INC F 

NOV NDF.CNT 
INC C 
SUB NDF,E 
SUB NDFSQ. ABASE 
DEC CT2 

SUB NDF,B 


BEQ CONT4> 


JNP LOOP4 


“GET B 
“STORE A 
"GET C 
"A*B 
• "GET D 

“STORE C 

"GET E 
"C4A*B 
"STORE D 

"STORK C»C+A*B 

"K*D 

"GET F 


"F+D«E 

"TEST INNER LOOP 
"STORE P-F+D*E 


"INC F ADDRESS 
"RELOAD COUNTER 

"REAJUST E ADDRESS 
"REAJUST A ADDRESS 

"TEST OUTER LOOP 
"REAJUST B ADDRESS 


“ TEST FOR K-NSUBS 


CONT4I NOV D, ABASE 
NOV NDF,CT2 

SUB NDF, ABASE "REAJUST A 

SUB# NSUB.CNTER "TEST X 

BEQ END) INC ABASE 
JNP LOOP3 


" FORTRAN I 300 RETURN 


ENDi 


RETURN 

3END 


APPENDIX I 


Listings of: 


FORO.FOR 

APFORO.FOR 


APTRN6.APM 
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TYPE TORO 
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0001 

0002 

0003 

0004 

0005 


0006 

0007 

0008 
0009 


0010 

0011 

0012 

0013 

0014 

0015 

0016 
0017 
0019 

0019 

0020 
0021 
0022 
0023 


CCOHDECK TRAN6 
C 

C ************************************************ 

c 

C This l8 Program FORO.POR, which represanta a portion of 
C Subroutine TRAN6 in the SPAR Library. Ueed to obtain timing 

C information for FORTRAN execution. 

C 

C Corresponding Programs arei 

C APFORO.FOR - FORTRAN of PORO with AP Calls 

C APTRN6.APN - FP120 Assembler Program replacement of 

C FORTRAN portions. 

C APTRN6.ABJ - Object code of APTRN6.APM 

C APTRN6.SAV - Linked version of APTRN6.AZ9I 

C 

C *««***••*••«****•••**•«***«*•**«•*****••**•**•«* 

c 

C SUBROUTINE TRAN6(MAP,H,NDF,NN0DES, ITRANS,T) 

C 4/73 WD WHETSTONE 

COMNON/CONSTR/ JT, JDF,JDDF, INEX( 6 ) , MEXIN{ 6 ) , KSYM( 3 ) 
COMMON/TEMPS/ 

$ HKL(6,6),GKLTL(6,6) 

COMMON/LOCAL/S( 6,6,10) 

DIMENSION MAP(10),H(6,6,1),ITRANS(4),T(3,3,1) 

DIMENSION ITIM1( 2 ) , ITIM2( 2 ) 

C H(NDF,NDF,1) WHERE THE TOTAL DIM-NDF*NDP*KSIZE 
C SINSERT SYSCOM>A$KEYS 

CALL ASSI(a»( 4, • RKltSFILE.DAT ',13, ’RDO’ ) 

CALL ASSIGW(2, 'RKl: TFILE.DAT • ,0, 'RDO' ) 

CALL ASSIGN(3, 'RKl sHFILE.DAT' ,0, 'RDO' ) 

44 F0RMAT(F16.7) 

C R£AD(7,44)Z,S 

C READ(8,44)Z,T 

C READ(9,44)Z,B 

WRITE(6,3001) 

3001 F0RMAT( ' INPUT THE NUMBER IF TIMES TO BE EXECUTED* ) 

READ(5,*) NTEST 

READ(4,44) (((S(I,J,K),I-1,6),J-1,6),K-1,10) 

READ<2,44) ( ( T( I , J, 1 ), I-l, 3 ), J-1, 3 ) 

READ<3,44) ( ( H( I, J, 1 ), I-l, 6 ), J-1, 6 ) 

CPU1-SECNDS(0. ) 

CALL GTIM(ITIMl) 

DO 3000 ITEST-1, NTEST 
NDF-3 
NNODES-3 
DO 4 1-1,4 
4 ITRANS(I)-! 

KX-1 
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0024 DO 7 X-1.10 

0025 NAP( I )-KX 

0026 7 KX— XX 

0027 DO 9 X-1,6 

0028 9 XNEX(X)-1 

0029 N-0 

0030 50 PORKAT( ' NUMBER OF SECONDS PAST MIDNIGHT - '.FIS. 7) 

0031 DO 2000 L-l.NNODES 

0032 IJ.-ITRANS(L) 

0033 DO 2000 X-l.L 

0034 XK-ITRANS(K) 

0035 N-N-l-1 
C 

0036 DO 100 J-1,6 

0037 DO 100 1-1.6 

0038 GKLTL(I,J)- .0 

0039 100 HXL(I.J)- .0 

C FORM HKL- TK( TRANSPOSE ) *GKL *TL 

C FIRST. GXL*TI., 

C 

0040 DO 1100 J-1.3 

0041 DO 1100 1-1.3 

0042 DO 1100 M-1.3 

0043 TUUr-T(M. J.LI,) 

0044 GKLTL( I. J)- GKLTL( I. J) 4S( I. M.N) *TU(J 

0045 GXLTL( I.J+3)- GKLTI.( I.J+3) +S( I.M+3.N) *TLMJ 

0046 GXLTL(I+3. J)- GKLTI.(I+3. J) +S(I+3. M.N) *TLKJ 

C047 1100 GKLTL(I+3,J+3)- GKLTL( 1+3, J+3 ) +S( I+3,M+3,N) *TLHJ 

C 

C TK( TRANSPOSE )*(GKL*TL) 

0048 DO 1200 J-1,3 

0049 DO 1200 1-1.3 

0050 DO 1200 N-1,3 

0051 TXMI-T(M,I,XX) 

0052 HXL( I, J)- HXL( I, J) +TXMI*GXLTL( M, J) 

0053 HXL( I.J+3)- HXI.( I.J+3) +TXMI*GXLTL( M.J+3) 

0054 HKL(I+3, J)- HXL(I+3, J) +TXMI *GXLTI,( M+3 , J) 

0055 1200 HXL(I+3, J+3)- BK1.( 1+3, J+3 ) +TXMI*GXLTL(M+3, J+3 ) 

C 

C TRANSPOSE, IF REQ. 

0056 1F(MAP(N).GT.0) GO TO 1400 

0058 DO 1300 J-2,6 

0059 JM-J-1 

0060 DO 1300 1-1, JM 

0061 EIJ- RKL(I.J) 

0062 HXL(I,J)-HKL( J.I) 

0063 1300 HKL(J,I)-EIJ 

0064 1400 LOC-IABS(MAP(N)) 
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0065 

0C67 

0068 

0069 

0070 

0071 

0072 

0073 

0074 
0C75 

0076 

0077 

0078 

0079 

0080 
0081 
0082 

0083 

0084 

0085 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

0093 


1P(NDF.LT.6) GO TO 1600 
DO 1500 J-1,6 
DO 1500 1-1,6 

1500 H(I,J,LOC)- H(I,J,IiOCHHlCL(I,J) 

GO TO 2000 
1600 DO 1700 I-1,NDF 
NROW-INEX( I ) 

DO 1700 J-1,NDF 
NCOL-INEX( J) 

1700 H( I , J, LOG )-H( I , J, IOC )+Hiai( NROW, NCOL ) 

2000 CONTINUE 
3000 CONTINUE 

CALL GTIM(ITIM2) 

CPU2- SECNDS(0. ) 

WRITE(6,50) CPUl 
WRITE(6,50) CPU2 

CALL CVTTIM( ITIia , IRR, INI , ISE , ITI ) 

WRITE(6,70) IHR,IMI,ISE,ITI 
70 PORMAT( • TIME - ’,12, • t ' ,12, * j M2, ' t * ,12 ) 

CALL CVTTIM( ITIM2 , IHR, IMI , ISE , ITI ) 

WRITE(6,70) IHR, IMI, ISE, ITI 
CPU - CPU2 - CPUl 
WRITE(6,99)CPU 

99 FORMAT( ’ EXECUTION TIME -’,F12.7,' SECONDS') 

MRITE(6,88) ( H( 1, 1, 1 ), I-l, 36 ), ( GKLTL( 1, 1 ) , I-l, 36 ), 
X(HKL(I,1), 1-1,36) 

88 FORMAT(3E16.7) 

STOP 

END 


ns 
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GCONDECK TRAN6 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

0001 

0002 

0003 

0004 

0005 

0006 

C B(NDP,NDP,1) WHERE THE TOTAL DIN-NDF«NDF*K5IZE 



C 

$ INSERT SYSCON>A«KEYS 

0007 


CALL APCLR 

0008 


CAUi APPUT( 0,15000, 1.1) 

0009 


CALL APPUT( 0,15500, 1,1) 

0010 


DO 20 L00P»1,200 

0011 


IADRO-( LOOP-1) *9 

0012 


CALL APPUT(XADRO,1500CH-LOOP,1,1) 

0013 


IADRO-( LOOP-1) *6 

0014 


CALL APPUT(IADRD,155O0tL0OP,l,l) 

0015 

20 

CONTINUE 

0016 


CALL APWD 

0017 


DO 18 LOOP^l,6O0 

0018 


CLR(LOOP)-0.0 

0019 

18 

CONTINUE 

0020 


CALL APPUT(CLR,0,600,2) 

0021 


CALL ASSICai(l, • RKliSFILE.DAT* ,0, *RD0' ) 

0022 


CALL ASSlCai( 2, • RKliTFILE.DAT ’,0, RDO' ) 

0023 


CALL ASSIGN(3, • RKliHFILE.DAT ’,0, ’RDO’ ) 

0024 


44 PORMAT(F16.7) 


C 

READ(7,44)Z,S 


C 

READ(8,44)Z,T 


nila is Program APF01H).F0R Which contains tha AP Calls 
as rsplacsments for FORTRAN Cods in FORO.FOR. Rsprsssnts 
a portion of subroutine TRAN6 in the SPAR Library. 

Obtains timing information for AP execution. 

Corresponding Programs are> 

FORO.FOR 

APTRN6.Ant 

AFTRN6.ABJ 

APTRN6.SAV 


SUBROUTINE TRAN6(NAP,B,NDF.NN00ES, ITRANS,T) 

4/73 WD WHETSTONE 

CONNCm/CONSTV •rr,JI»P,JDDF,INEX(6),KEXIN(6),KSYN(3) 
CONNON/TQfPS/ 

$ HKL(6,6),aa.TL(6,6) 

C0NN0N/L0CAL/S( 6 , 6 , 10 ) 

DIMENSION NAP(10),H(6,6,1).XTRANS(4),T(3,3,1),C'LR(600) 
DIMENSION ITINl( 2 ) , ITIM2( 2 ) , ITIN3( 2 ) , ITIM4( 2 ) . ITIie( 2 ) 
DIMENSION XTIN6(2),ITIN7(2) 
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FORTRAN IV V02.5-2 Tue 04-NOV-80 10 t 27 i 20 PACE 002 

C READ(9,44)Z,R 

0025 WR1TE(6,3001) 

0026 3001 PORNAT( ' INPUT THE NUMBER IP TINES TO BE EXECUTED* ) 

0027 READ(5,*) NTEST 

0028 READ(1,44) ( ( ( S( I, J,K), 1-1,6 ), J-1,6 ),K-1, 10 ) 

0029 READ(2,44) ( (T( X, J,1 ), I-l, 3 ), J-1, 3 ) 

0030 READ(3,44) ( (H( I, J, 1 ), 1-1,6 ), J-1,6 ) 

0031 CPU1-SECNDS( 0 . ) 

0032 CALL GTZM(1TIM7) 

0033 DO 3000 ITEST-1 , NTEST 

D CALL GTZM(ITZNl) 

0034 NDF-3 

0035 NNODES-3 

0036 DO 4 1-1,4 

0037 4 ITRANS( I )-l 

0038 XK-1 

0039 DO 7 Z-1,10 

0040 NAP( I )-KX 

0041 7 KX— KK 

0042 DO 9 Z-1,6 

0043 9 INEX( I )-l 

0044 N-0 

0045 50 FORMAT( ' NUMBER OF SECONDS PAST MIDNIGHT - *,F15.7) 

C 

0046 CALL APPUT(INEX, 0,6,1) 

0047 CALL APPUT(ITRANS(1),6,4,1) 

0048 CALL APPUT(S(1,1,1),82,360,2) 

0049 CALL APPUT(H(1,1,1),442,36,2) 

0050 CALL APPUT(T(1,1,1),488,9,2) 

0051 CALL APPUT(NAP(1), 478, 10,1) 

0052 CALL APND 

D CALL GTIM(ITIN2) 

0053 CALL APTRN6(48B,NDF,NN0DES,478) 

0054 CALL APHR 

D CALL GTIM(ITZN3) 

0055 CALL APGET(H(1,1,1),442,36,2) 

0056 CALL APWD 

D CALL GTIM(IT1M4) 

C 

c ******************************************** 

c 

C FORTRAN Code replacement - Begin 
C 

C DO 2000 L-1,NN0DES 

C LL-ITRANS(L) 

C DO 2000 K-1,L 

C KX-ITRANS(K) 

C N-N+1 
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FORTRAN XV V02.5-2 Tue 04-NOV-80 10 t 27 t 20 


PACZ 003 


CC 

C DO 100 J-1,6 

C DO 100 1-1,6 

C GKLTL(I,J)- .0 

C 100 HKL(I,J)- .0 

cx: Pom BXL- nc( transpose) *001 *tl 

cx: FIRST, aCL*TI.. 

CC 

C DO 1100 J-1,3 

C DO 1100 1-1,3 

C DO 1100 N-1,3 

C TUtJ-T(N,J,U.) 

CNC GKT,TL( I, J)- GKLTLC I, J) -t-S( I, N,N) •TLHJ 

C GKLTL( I,J+3)- GKLTI.( I,J+3) +S< I,lf<-3,N) miMJ 

C Gia.TI,(I+3, J)- Gia.TL(I+3, J) +S(I+3, N,N) *TUIJ 

C 1100 GXLTL(I43,J-l-3)- GKXTL( 143, J+3 ) 4S( X43,M43,N) »TU(J 

CC 

C TK( TRANSPOSE )*(aCL*TL) 

C DO 1200 J-1,3 

C DO 1200 1-1,3 

C DO 1200 N-1,3 

C TKMI-T(M,I,iaC) 

C HKL( I, J)- ma,( I, J) 4TKMI*Gia.TL( M, J) 

C BKL( I,J43)- HKL( I,J43) 4TKMI*GKLTL( M,J43) 

C HKL(I43, J)- Hia,(I43, J) 4TKMI *GKLTL( H43 , J) 

C 1200 HKL(I43. J43)- HKL( 143, J43 ) 41XMI*Gaa,TL(H43, J43 ) 

CC 

CC TRANSPOSE, IF REQ. 

C IF(MAP(N).GT.O) GO TO 1400 

C DO 1300 J-2,6 

C JM-J-1 

C DO 1300 I-1,JN 

C EIJ- BKL(I,J) 

C 00.(1, J)-HKL( J,I) 

C 1300 BK1,(J,Z)-EIJ 
C 1400 LOC-IABS(NAP(N)) 

C IF(NDF.LT.e) GO TO 1600 

C DO 1500 J-1,6 

C DO 1500 1-1,6 

C 1500 H(I,J,LOC)- H(I,J,L0C)4HKL(I,J) 

C GO TO 2000 

C 1600 DO 1700 I-1,NDF 
C NROIf-INEX( I ) 

C DO 1700 J-1,HDF 

C NCOIi-INEX(J) 

C 1700 H(I, J,L0C)-H(1, J,L0C)4HKL(NR0W,NC0L) 

C 2000 CONTINUE 
C 
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FORTRAN TV V02.S-2 Tu« 04-Nov-BO 10 i 27 i 20 PAGE 004 

C FORTRAN Codt xtplacCMnt - End 

C 

c ••*•*«•*******••****•**••.*••*******•**•*• 

c 

D CALL GTIN(ITIHS) 

D NR2TE(6,4000) ITEST 

D CALL CVTTZN(ZTIM1,XHR,XN2,XSE,ZTI) 

D WRXTE(£,70) XRR, XMX , XSE, XTX 

D CALL CVTrXM(XTZM2,XHR,XMX,XSE,XTX) 

D WRXTE(6,70) XHR, XNZ, XSE, XTX 

D CALL CVTTXN(XTIM3,XHR,XNI,XSE,XTX) 

D WRXTE(6,70) XRR,I1I1,XSE,XTX 

D CALL CVTTXM(XTXN4,XBR,XNX,XSE,XTX) 

D WRXTE(6,70) XBR,XNX,XSE,XTX 

D CALL CVTTXN( XT1M5 , XRR, XNX , XSE , XTX ) 

D NRITE(6,70) XHR,XMX,XSE,XTX 

0057 3000 CONTXNUE 

0058 CALL GTXM;XTXM6) 

0059 CPU2- SECNDS(0. ) 

0060 WRXTE( 6,50) CPUl 

0061 WRXTE(6,S0) CPU2 

0062 CALL CVTTXM(XTIM7,XHR,XNX,XSE,XTX) 

0063 WRITE(6,70) XHR, XNX, XSE, XTX 

0064 4000 P0RMAT( ' CYCLE - ',X5) 

0065 70 PORMAT( • TIME - M2, • I *,X2, ' J *,X2, • « *,X2) 

0066 CALL CVTTXM(XTIM6,XBR,IMX,X5E,XTI) 

0067 WRXTE(6,70) XBR, XMl , XSE, XTX 

0068 CPU - CP02 - CPUl 

0069 WRITE(6,99)CP0 

0070 99 FORMATC ' EXECUTION TIME -•,P12.7,’ SECONDS') 

0071 WRITE(6,88) ( H( X , 1, 1 ), I-l, 36 ) 

0072 88 P0RMAT(E16.7) 

0073 STOP 

0074 END 
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"MODIFIED 29-MAY-80 TO USE M.D. INSTEAD OF TMRAM FOR ADDRESS 
"TABLE 


*TITLE TRAN6 
♦ENTRY APTRN6,4 


"THIS ROUTINE PERFORMS THE TRANS SUBROUTINE LOCATED IN THE 
" SPAR LIBRARY. IT IS CALLED EXTENSIVELY IN THE K AND M PROCESSORS. 

M 

" AUTHOR: K. FESON 
" DATE: MAY 1979 


— USUAGE 


" FORTRAN: CALL APTRN6<TBASE,fCF,NN0DES,MAPBAS) 


" tPAGE 

M 


•I 


M 


MAIN DATA MEMORY MAP 


♦ t 

* * 

* T ARRAY OR B3<LE> 

t t 

* * 

t t 

t t 

t K4<LREC5) OR MAP 

t t 

* * 

t t 

t H<NDF*NDF»KSI2E> t 

* t 

* » 

♦ SLOC<360> t 

♦ t 

♦ t 

t GKLTL <30 * 

t * 

t t 

t HKL( 3S > » 

» » 


(STARTING ADDRESS) 

♦ <-- LR (ADDRESS RELATIVE TO BASE) 
442+NDF»NDF*KS I 2E+LREC5 
t <— L- (ADDRESS RELATIVE TO BASt) 
442+NDFCNDF»KSI2E 

442 

82 


46 


» » 

t 1TRAHS<4> t 

»: t 

% % 

% n€x<6> * 

t t 


le 


€ 


e 


ARRAY DESCRIPTIONS 


INEX(6) INTEGER ARRAY. TRANSFERED ONLY ONCE 

ITRANSC* ,» INTEGER ARRAY. TRANSFERED ONLY ONCE. 

HKL<36) REAL ARRAY. HE\ER TRANSMITTED. 

GKLTL<36> REAL ARRAY. NEVER TRANSMITTED. 

SLOC<360) REAL ARRAY. TRANSFERED ONCE FOR EVERY CALL TO 

THE ROUTINE. 


H<N0F,NDF,KSI2E). . REAL ARRAY. TRANSFERED AND RETURNED FOR EVERY 

CALL TO THE ROUTINE. 

K4(^LREC5> OR MAP , INTEGER ARRAY. TRANSFERED ONLY AFTER BEING READ 

FROM THE DATA BASE. 

B3(LE) OR T<3-3,.) REAL ARRAY. TRANSFERED ONLY AFtER BEING READ 

FROM THE DATA BASE. 


TABLE MEMORY USAGE 

A TABLE OF INDIRECT ADDRESS OFFSETS MUST BE LOADEO AT 
TMRAM LOCATION 4096 (DECIMAL). 

THE TABLE IS USED TO CALCULATE AN ADDRESS OFFSET FOR 
AN ARRAY OF THE FORM X(3,3-LOC) THE VALUE 'LOC* IS 
ADDED TO 4096 TO FlfO THE ADDRESS OFFSET IN THE 
TACl-E FOP EXAMPLE IF L0C=3, THE TABLE LOCATION 
4096+3 OR 4099 IS READ AND A VALUE 18 IS LOADED 
AS THE ADDRESS OFFSET FOR ARRAY Y(3,3,3>. THIS 


REPRESEHTS THE ADDRESS OF XU,1.3> WHEN ADDED TO THE 
BASE ADDRESS C«^ X< 1,1,1) IN THE MAIN DATA MEMORY 


II 


II 

4096 

— 

0 

N 

4097 

— 

0 

H 

4098 

— 

O 

M 

4099 

— 

18 

M 

4100 

— 

2'’ 

M 

4101 

— 

36 

M 

4102 

— 

45 

II 

4103 

— 

54 


4104 

— 

63 

H 

4105 

— 

72 

M 

4106 

— 

81 


4107 

— 

90 

N 

4108 

— 

99 

N 

4109 

— 

108 

N 

4110 

— 

109 

M 

4111 

— 

117 

N 

4112 

— 

126 


N 


M 

M 

N 

" TABLE ADDRESS OFFSET ADDRESS 

M 


" < THE TABLE SHOIA-D COlfrAIN SEVERAL HWCRED ENTRIES) 


" TABLE fCMORY USAGE 


" ANOTHER TABLE OF INDIRECT ADDRESS MUST BE LOADED AT 
“ TMRAM LOCATION 4196 (DECIMAL). 

" THE TABLE IS USED TO LOCATE AN ADDRESS OFFSET FOR AN 
" ARRAY OF THE FORM X(6>NC0L) THE VALl€ CT ‘MCO.' IS 
" ADDED TO 4196 TO FIND THE TABLE ICMORY LOCATION OF 
" THE ADDRESS OFFSET 


4196 — 0 

4197 — 6 

4198 — 12 

4199 — 18 

4200 — 24 

4201 — 30 

4202 — 36 

4203 — 42 
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« 4204 — 42 

« 4205 — 54 

'■ 4202 — 60 

« 4207 “ 66 

It 

ft 

tt ... 

M 


TABLE ADDRESS ADDRESS OFFSET 

(THE TABLE SHOULD CONTAIN SEVERAL HUftDRED ErmiES) 


” S-PAD PARAMETERS 


TBASE »EQU 0 "BASE ADDRESS OF T ARRAY 
NOF #EQU 1 -DEGREES OF FREEDOM PER JOINT 

IMDOES »EQU 2 "NUMBER OF NODES IN ELEMENT 

T2INC *EQU 3 “OFFSET ADDRESS TO LOCATE T< 1 , I , KK > 

THTSIX 4E0U 3 "CONSTANT OF 36 

TUNTYl »EQU 3 "CONSTANT OF 21 

EIGHTN »EQU 3 "CONSTANT OF 18 

MAPBAS *EQU 3 "BASE ADDRESS OF MAP ARRAY (PASSED AS ARGIW 

NROUl $EQU 3 "LOOP COUNTER 

GKLTLB »E^ 4 "BASE ADDRESS OF CKLTL ARRAY 

BASEGK $EQU 4 "BASE AOWESS OF CKLTL ARRAY 

HBASE tEQU 4 "BASE ADDRESS OF H ARRAY 

HKINC »EOU 4 "ADDRESS POINTER OF H ARRAY 

THREE tEQU 5 "CONSTANT OF 3 

HKADDR tEQU 5 "ADDRESS POINTER OF HKL ARRAY 

N tEQU 6 "INDEX FOR MAP ARRAY 

XNT tEQU 6 "LOOP COUNTER 

J tEQU 6 "LOOP COL»tTER 

GKLTLA tEQU 7 “ADDRESS POINTER FOR CKLTL ARRAY 

CKLTL 1 tEQU 7 "ADDRESS POINTER FOR CKLTL(1,J> ARRAY 

HKLl tEQU 7 "ADDRESS POINTER FOR HKL(I,J> 

HADDRl tEQU 7 "ADDRESS POINTER FOR H ARRAY 

SIADDR tEQU 10 "ADDRESS POINTER FOR S(I,M,N) ARRAY 

GKLTL2 tEQU 10 "ADDRESS POINTER FOR GKLTL(2-J) 

HKL2 tEQL 10 "ADDRESS POINTER FOR HKL(J,l) ARRAY 

HADDR2 tEQU 10 "ADDRESS POINTER FOR H(I,J.LOC) 

TMA tEQU 11 "TABI.E ADDRESS POINTER 

S2ADDR tEQU 11 "ADDRESS POINTER FOR S(I,M+1,N> 

CKLTL3 tEQU 11 "ADDRESS POINTER FOR CKLTL(3,J) 

HKTMP tEQU 11 "ADORESS POINTER FOR HKL ARRAY 

I TRAN tECHJ 12 "BASE ADDRESS OF I TRAN ARRAY 
LL tEQU 12 "INDICY FOR T ARRAY * ITRAN(L) 

KK tEQU 12 "INDICY FOR T ARRAY » 2TRAN(K) 


m 
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S3ADDR 

$EQU 

1£ 

"ADDRESS POINTER FOR S<I,M+S,N> 

HKL 

t£QU 

12 

"ADDRESS POINTER FOR HKL ARRAY 

HADDR 

tEQU 

1£ 

"ADDRESS POINTER FOR H ARRAY 

K 

$EQU 

13 

"INNER LOOP COUNTER FOR 2000 LOOP 

ICNT 

$EQU 

13 

"LOOP COUNTER 

JM 

tEQU 

13 

"LOOP COUNTER 

CNT 

$EQU 

13 

"LOOP COUNTER 

ZERO 

fEQU 

13 

"CONSTANT OF 0 

I 

$EQU 

13 

"LOOP COUNTER 

L 

$EQU 

14 

"OUTER LOOP COUNTER FOR £000 LOOP 

DPA 

$EQU 

14 

"DATA PAD REGISTERS BASE 

LOC 

^EQU 

14 

"TEMPORARY STORAGE OF MAP<N) 

INEX 

$EQU 

14 

"ADDRESS POINTER FOR INEX ARRAY 

SIX 

$EQU 

15 

“CONSTANT OF 6 

TUNC 

$EQU 

16 

"OFFSET ADDRESS TO LOCATE T<1,1,KK) BASE 

SADDR 

$EQU 

16 

"TEMPORARY STORAGE OF S BASE 

N£7 

tEQU 

16 

"CONSTANT OF £7 

NCOL 

$EQU 

16 

"LOOP COUNTER 

MAP 

$EQU 

16 

"ADEM?ESS OFFSET FOR MAP ARRAY 

SBASE 

$EQU 

17 

"BASE ADDRESS OF S ARRAY 

BASEHK 

$EQU 

17 

"BASE ADDRESS OF HKL ARRAY 

TADRl 

$EQU ; 

15000. 

"ADDRESS OF ADDRESS TABLE 1 

TADR£ 

fEQU : 

15500. 

"ADDRESS OF ADDRESS TABLE 2 


FORTRAN: H=0 

DO £000 L=l..NNODES 
LL=ITRAHS<L) 


APTRN&: LDDPA; DB=1£. 

CLR N 

DEC N; DPX<£XSPFN 
CLR Li DPX<0XSPFN 
MOV MAPBAS.MAPBAS; DPY<0XSPFN 
LDSPI SBASE; DB=46. 

MOV SBASE, SBASE; DPXOXSPFN 


"GET HIGH DATA PADS 
"SET N=0 
"INIT N TO -1 
“INIT L 

"STORE BASE OF MAP ARRAY 
"LOAD SBASE - 36 
"SAVE IN HIGH DATA PAD 


LOOPl: LDSP7 L; DB=DPX<0) 

LDSPI ITRAN; DB=6. 
ADO# L,1:RAN; SETMA 
LDSPI TMA; OB=TADRl 
CLR » ; DPXCIXSPFN 
LDSPI LL; DB=MD 
ADD# LL,TMA; SETMA 
INC L, DPX<0XSPFN 
NOP 


"RESTORE L 

LOAD BASE OF I TRAN ARRAY 
"GET LL=ITRAN<L) 

"GET BASE OF TABLE 

"RESET INNER LOOP COUNT 
"SAVE LL 

FIND ADDRESS OFFSET IN TABLE 
"SAVE L+1 IN HIGH DATA PAD 
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LDSPI T£lNC.iDB-MD “SAVE ADDRESS OFFSET 

MOV T£lNC,TaiNCi DPY<2KSPFN “SAVE IN HIGH DATA PAD 


M 


II 


FORTRAN' DO £000 K«1,L 
KK*ITPANS<K> 


'• >. N=N+1 IS PERFORMED LATER IN THE CODE > 

14 


II 


LOOP£^ LDSPI I TRAN; DB»6 
LDSPI K.; DB«OPXa> 

ADD# KMTRAN; SETMA 
LDSPI TMA,i DB=TADR1 
NOP 

LDSPI KK; DB=MD 
ADD# KK,THA,i SETMA 
INC Ki DPXaXSPFN 
NOP 

LDSPI TUNC,; DB=MD 

MOV T 1 1 NC . T 1 1 NC DP YU X SPFN 


“LOAD BASE OF I TRAN ARRAY 
"RESTORE K 
“GET KK=ITRAN<K> 

“GET BASE OF TABLE 
“WAIT FOR M,D, 

“SAVE KK 

"FIND ADDRESS OFFSET IN TABLE 
"SAVE K+1 IN HIGH DATA PAD 
"WAIT FOR M,D. 

"STORE ADDRESS OFFSET 
“SAVE IN HIGH DATA PAD 


" GET TLMJ=TCM, J,LL> 

II 

" TBASE+T£INC GIVES THE ADDRESS OF T<1,1,LL> 
" FOUND IN THE PREVIOUS SECTION 

II 

" THE NINE VALUES OF T FOR M=1.3 AND J=l,3 
“ ARE STORED IN THE DATA PAD X REGISTERS 


LDSPI T£INC.; DB=DPY(£> 

LDSPI THREE; DB»3, 

ADD TBASE,T£INC; SETMA 
CLR DPA, SETDPA 
INC T£INC; SETMA 
DPXCOXDB. DB=MD 
INC T2INC. SETMA 
DPXC 1 XDB. OB*MD 
INC T£INC; SETMA 
DPX^aXDB. DB=MD; 

ADD THREE, DPA; SETDPA 
INC T£INC; SETMA 
DPXCOXDB; DB»MD 
INC T2INC; SETMA 
DPXaXDB, DB=MD 


"RESTORE T2INC 
“LOAD CONSTANT 
"GET T<1,1,LL> 

"CLEAR DATA PAD BASE AO 
"GET T>;£,t,LL) 

"SAVE Tk 1,1,LL> 

"GET Tk3,l,LL) 

"SAVE T(£,1,LL> 

"GET T<1,£,LL) 

"SAVE Tk3,l,LL> 

"GET NEXT SET OF DATA PADS 
"GET T<£,£,LL) 

"SAVE Ta,£,LL) 

"GET T(£,3.LL> 

“SAVE T<£,£,LL> 
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INC T£1NC,: SETMA 
DPX<£KDB; DB=MD: 

ADD THREE, DPA.i SETDPA 
INC T£INC,i SETMA 
DPX< 0 HMD 
INC T2INC; SETMA 
DPX< 1 XMD 
NOP 

DPX(£Kmq 


" PERFORM THE 11 00 LOOP CALCULATIONS IN A SERIES OF STEPS 

II 

'• BECAUSE THE INDIVIDUAL CALCULATIONS ARE INDEPENDENT, EACH GIVEN 
" CALCULATION IS PERFORMED IN A SEPERATE IDENTICAL ROUTINE 

H 

’■ THE INNER MOST LOOP ACTUALLY PERFORMS THE FORTRAN M LOOP INTERNALLY 
■' WITH 'HE CALCULATIONS 

II 

THE FORTRAN I LOOP IS THE INNER LOOP FOR THE ROUTINE 

It 

'■ THE J LOOP IS OUTSIDE BOTH AS NORMAL 


"GET T<l,3,LL> 

"SAVE T<3,£,LL) 

"GET NEXT SET OF DATA PADS 
"GET T(S,3,LL) 

"SAVE Ta,3,LL> 

"GET T<3,3,LL) 

"SAVE T<£’,3,LL) 

"SAVE T<3,3,LL> 


" FORTRAN^ DO 160 J=l,6 
" DO 100 1=1,6 

" GKLTL<I, J)=.0 

" 100 HKL(1,J>= .0 

" DO 1100 J=l,3 

" DO 1100 1=1,3 

" DO 1100 M=l,3 

" TLMJ=TCM,J,LL) 

" GKLTL<I,J> = GKLTL<I,J> + S< I,M,N)<-TLMJ 

" GKLTL<I,J+3) = GKLTL<I,J+3) + S< I , M+3 , N ):»:TLM J 

" GKLTL(I+3,J> = GKLTL<I+3,J) + S< I+3,M,N>*TLMJ 

" 1100 GKLTL<I+3,J+3> = GKLTL< 1+3, J+3 ) + S<I+3,M+3,N) 

II 

" CALCULATE GKLTL<I,J> = GKLTL<I,J) + Sa,M,N)*TLMJ 


LDSPI SIX; DB=6. 

CLR DPA; SETDPA 
MOV THREE, JCNT 
LDSPI GKLTLB; DB=46. 
MOV GKLTLB, GKLTLA 
DEC GKLTLA 
LDDPAi DB=1£. 


"LOAD CONSTANT 

"CLEAR DATA PAD BASE 
"LOAD J LOOP COUNT 
"LOAD BASE ADR OF GKLTL 
"LOAD BASE ADDRESS 
"LOOP SET-UP 
"GET HIGH DATA PADS 
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L00P3 ^ 


L00P4 = 


LOSPl SBASE; DB=DPX<3> 

LDSPI THTSIX; 0B=36. 

ADD THTSIX, SBASE; DPX<3XSPFN 
LDDPA; DB=0, 

MOV SBASE, SADDR 
MOV SADDR,S1A0DR; SETMA 
MOV S1ADDR,SEADDP 
ADD SIX,S£ADDR; SETMA 
MOV S2AD0R,S3ADDR. 

FMUL DPX<0>.MD 
ADD SIX,S3ADDR; SETMA. 

FMUL 

FMUL DPX<1),MD 
MOV THREE,ICNT; 

FMUL; DPY<0XFM 
FMUL DPX<£),MD 
INC SIADDR; SETMA; 

FADD FM,DPY<0); FMUL 
FADD; FMUL 
INC S£ADDR; SETMA; 

FADD FM,FA 
FMUL DPX<0),MD; 

FADD 

INC S3ADDR; SETMA; FMUL 
FMUL DPXa>,MO; 

DEC ICNT 

FMUL; DPY<0XFM; 

INC GKLTLA; SETMA; MKFA; 
BNE L00P4 


•'RESTORE SBASE 
"LOAD CONSTANT 
"GET Sa,M,N+l > BASE 
"GET LOW DATA PADS 
"LOAD S<I,M,N> BASE 
"LOAD Sa,M,N) 

"LOAD S(I,M+1,N> 

" S1»TLMJ(1,J,LL) 

"LOAD Sa,M+£,N> 

"PUSH 

" TLM J( £ , J , LL )*S< I , M+ 1 , N ) 
"LOAD INNER COUNT 
"STORE SI PRODUCT 

" S< I , M+£ , N )*TLM J< 3 , J , LL ) 
"GET NEXT SI 

"S£ PRODUCT+ SI PRODUCT 

"GET NEXT S£ 

"S£+S1+S3 (PRODUCTS) 

"Sl + UTLMJ 

"GET NEXT S3 
"S£+UTLMJ 
"CHECK LOOP 
"SAVE Sl+1 PRODUCT 
"STORE RESULT 


THE INNER I AND M LOOPS ARE COMPLETE 
REA JUST AND CHECK THE OUTER J LOOP 


ADD THREE, DF A; SETDPA 
DEC JCNT 

BEQ CONTI; ADD THREE, GKLTLA 
JMP L00P3 


"GET NEXT DATA PAD UNITS 
"CHECK OUTER LOOP 

"PFAJUST GKLTLA ADDRESS 


CALCULATE GKLTL(I+3,J) = GKLTL(I+3,J) + S( 1+3, M, N >4:TLMJ 
TLMJ IS STILL IN THE DATA PADS 


CONTI ^ 

CLR DPA; SETDPA 

"GET FIRST SET OF 1 


MOV THREE, JCNT 

"LOAD COUNTER 


MOV GKLTLB.. GKLTLA 

"LOAD BASE ADDRESS 


ADD THREE, GKLTLA 

"GET GKLTL(I+3,J) ADDRESS 

• 

DEC GKLTLA 

"LOOP SET UP 
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MOV SBASE,SADDP 

"LOAD S BASE 


ADD THREE, SADDP 

"GET S<I+3,M,N) ADDRESS 

L00P5 

MOV SADDR.SIADDR; SETMA 
MOV SIADDR.SSADDR 

"LOAD SI ADDRESS 


ADD SIX,S£ADDR^ SETMA 
MOV S£ADDR,S3ADDR; 

"LOAD S£ ADDRESS 


FMUL DPX<0),MD 

"DO SUTMLJ 


ADD SIX,S3ADDRi SETMA: 
FMUL 

"GET S3 ADDRESS 


FMUL DPXa>,MD 

"DO SE4TLMJ 


MOV THREE.ICNT. 

"LOAD INNER COUNT 


FMUL, DPY<.eXFM 

"STORE Sl*TLMJ 

LOOPb ^ 

FMUL DPX<2>,MD 

"DO S3ATLNJ 


INC SIADDR, SETMA: 

"GET NEXT SI ELEMENT 


FADD FM,DPY<,0); FMUL 
FADD; FMUL 

"S£ PRODUCT + SI PRODUCT 


INC SSADDR; SETMA; 

FADD FM,FA 

"GET NEXT S£ ELEMENT 


FMUL DPX<0),MD; 

FADD 

"SI + imMJ 


INC S3ADDR; SETMA; FMUL 

"GET NEXT S3 ELEMENT 


FMUL DPX<l),MD; 

"DO S£+1*TLMJ 


DEC ICNT 

"CHECK LOOP 


FMUL; DPY<0XFM; 

"SAVE Sl+1 PRODUCT 


INC GKLTLA; SETMA :MI< FA; 

"STORE RESULT 

M 

BNE LOOPS 


" RECHECK THE OUTER J LOOP 

It 



ADD THREE, DPA; SETDPA 

"GET NEXT DATA PADS 


DEC JCNT 

"CHECK OUTER LOOP 


b’EQ CONTE; ADD THREE, GKLTLA 

"RE A JUST GKLTLA ADDRESS 

It 

JMP LOOPS 



CALCULATE GKLTL<I,J+3) = GKLTL<I.J+3) + S(I , M+3, N >ATLMJ 


" TLMJ IS STILL IN DATA PADS 

H 


M 

CONT£^ CLR DPA, SETDPA 

"GET FIRST SET OF DATA 

MOV THREE, JCNT 

"LOAD OUTER COUNT 

MOV GKLTLB , GKLTLA 

"LOAD BASE ADDRESS 

LDSPI EIGHTN; DB=18. 

"LOAD CONSTANT 

ADD EIGHTN, GKLTLA 

"GET ADDRESS OF GKLTLU,J+3) 

DEC GKLTLA 

"LOOP SET UP 

MOV SBASE,SADDR 

"LOAD S BASE ADDRESS 

ADD EIGHTN, SADDR 

"GET S<I,M+3,N) AS BASE 

L00P7^ MOV SADDR, SIADDR; SETMA 

MOV S1ADDR,S£A0DR 

"LOAD SI WITH BASE 
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LOOPS ^ 


It 


ADD SIX,S£ADDR; SETMA 
MOV S2ADDR-S3A0DR: 

FMUL DPX<0),ND 
ADD SIX,S3ADDRi SETMA; 

FMUL 

FMUL DPX<1>,MD 
MOV THPEE^ICNT; 

FMUL; DPY<0KFM 
FMUL DPX<2> .MD 
INC SIADDP; SETMA; 

FADD FM,DPY<0 .V; FMUL 
FADD; FMUL 
iNC S£ADDR; SETMA; 

FADD FM..FA 
FMUL DPX<0),MD; 

FADD 

INC S3ADDR; SETMA; FMLH. 

FMUL DPX<1),MD; 

DEC ICNT 

FMUL; DPYCeXFM; 

INC GKLTLA; SETMA; MI< FA; 
BNE LOOPS 


“GET S£ 

"DO SIATLMJ 
"GET S3 

"DO S£»TLMJ 
"LOAD INNER COUNT 
"SAVE SIATLMJ 
"DO S3*TLMJ 
"GET NEXT SI ELEMENT 
"SI PRODUCT +S£ PRODUCT 

"GET NEXT S£ ELEMENT 
"S1+S2+S3 PRODUCTS 
"Sl+UTLMJ 

"GET NEXT S3 
"DO S£+1*TLMJ 
"TEST LOOP 
"SAVE Sl+l PRODUCT 
"STORE RESULT 


" CHECK THE OUTER J LOOP 


ADD THREE, DPA; SETDPA 
DEC JCNT 

BEQ CONT3; ADD THREE, GKLTLA 
JMP LOOP? 


"GET NEXT DATA PAD SET 
"CHECK LOOP 

"RE A JUST ADDRESS 


" CALCULATE GKLTL(I+3, J+3) = GKLTLCI+3, J+3 ) + SU+3,M+3,N)*TLMJ 

t| 

" TLMJ IS STILL IN DATA PADS 


C0NT3 ^ 


LOOPS ^ 


CLR DPA, SETDPA 
MOV THREE, JCNT 
MOV GKLTLB, GKLTLA 
LDSPI TWNTYl; DB=£* 

ADD TUNTYl, GKLTLA 
DEC GKLTLA 
MOV SBASE,SADDR 
ADD TUNTY1,SAD0R 
MOV SADDR.SIADDR; SETMA 
MOV S1ADDR,S£ADDR 
ADD SIX,S£ADDR; SETMA 
MOV S£ADDR,S3ADOR; 

FMUL DPX<0>,MD 
ADD SIX,S3ADDR; SETMA; 


"RESET DATA PAD ADDRESS 
"LOAD OUTER COUNT 
"LOAD BASE 
"LOAD CONSTANT 
"GET GKLTL<I+3,J+3) ADDRESS 
"LOOP SETUP 
"LOAD BASE FOR S 
"GET S <I+3,M+3,N> ADDRESS 
"GET SI 

"GET S£ 

"SUTLMJ 
"GET S3 
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FMUL 

FMUL DPX(1>,MD 
MOV THREE.ICHT, 

FMUL; DPY<0XFM 
LOOPIO^ FMUL DPX(£>,MD 

INC SIADDR; SETMA; 

FADD FM,DPY<0); FMUL 
FADD; FMUL 
INC S2ADDR; SETMA. 

FADD FM,FA 
FMUL DPXCe^.MD; 

FADD 

INC S3ADDR; SETMA; FMUL 
FMUL DPXa)>MD; 

DEC ICNT 

FMUL; DPYCeXFM; 

INC GKLTLA; SETMA; MKFA; 
BNE LOOP 10 


"SS^TLMJ 

"LOAD INNER COUNT 
"SAVE SUTLMJ 
"DO S3»TLMJ 
"GET NEXT SI ELEMENT 

"SI PRODUCT + S2 PRODUCT 

•GET NEXT S£ ELEMENT 
"ADD SI + S£ +S3 PRODUCTS 
"DO Sl + UTLMJ 

"GET NEXT S3 ELEMENT 
"DO S2+UTLMJ 
"CHECK INNER COUNT 
"SAVE Sl+1 PRODUCT 
"STORE RESULT 


" CHECK THE OUTER J LOOP 

II 

ADD THREE. DPA; SETDPA "GET NEXT DATA PADS 

DEC JCNT "CHECK LOOP 

BEQ CONT4; ADD THREE • GKLTLA 
JMP LOOPS 

II 

" THIS ROUTINE PERFORMS THE 12O0 LOOP CALCULATIONS 


" GET TLMJ=T(M.J.KK> 

II 

" TBASE+TIINC GIVES THE ADDRESS OF T<1.1.KK> 

" FOUND IN THE PREVIOUS SECTION 

II 

" THE NINE VALUES OF T FOR M=1.3 AND J=1.3 
" ARE STORED IN THE DATA PAD X REGISTERS 

II 

" EACH CALCULATION IS PERFORMED INDEPENDENTLY IN A SEPERATE 
" AND INDENT I CAL ROUTINE JUST LIKE THE 1100 LOOP. 


FORTRAN^ DO 1200 J=-1..3 
DO 1200 1=1.3 
DO 1200 M=1.3 
TKMI=T(M.I.KK) 

HKLU.J) = HKL<I.J) +TKMI*GKLTL<M.J) 

HKL<I.J+3> = HKL<I.J+3) +TKMT.t,;,WLTL(M. J»3> 

HKL< I +3 . J ) = HKL< I +3 . J ) + TKM I »GKLTL< M-> 3 . J > 

1200 HKL<I+3.J+3) = HKL(I+3.J+3) + TKMI* GKLTU'M+3. J+3) 
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It 


C0NT4^ LDDPA; DB=1£ 

LDSPI TUNC; DB=DPY<1 > 

ADD TBASE,T1INC; SETMA 
CLR DPA; SETDPA 
INC TIIHC; SETMA 
DPX<0XDB; DB=MD 
INC TiINC; SETMA 
DPX-; 1 KDB; DB=MD 
INC TUNC. SETMA 
DPX<£ KDB; DB=MD; 

ADD THREE. UP A; SETDPA 
INC TUNC; SETMA 
DPX<0KDB. DB=MD 
INC TUNC; SETMA 
DPXaXDB; DB=MD 
INC TUNC; SETMA 
DPX<£XDB; DB=MD; 

ADD THREE. DPA; SETDPA 
INC TUNC; SETMA 
DPX<0XMD 
INC TUNC; SETMA 
DPX< I XMD 
NOP 

DPX<£XMD 


"GET HIGH DATA PADS 
"RESTORE TIINC 
"GET T<1.1.KK) 

"CLEAR DATA PAD BASE AD 
"GET T<£,1.KK> 

"SAVE T(l.l.KK) 

"GET TO.l.KK) 

"SAVE T<£.1.KK) 

"GET Ta.£.KK) 

"SAVE T<3.1.KK> 

"GET NEXT SET OF DATA PADS 
"GET T(£.£.KK> 

"SAVE T<1,£.KK) 

"GET T<£.3.KK> 

"SAVE T<£.2.KK) 

"GET T<1.3.KK) 

"SAVE T<3.£.KK> 

"GET NEXT SET OF DATA PADS 
"GET T<£.3,KK'> 

"SAVE T<1.3.KK) 

"GET T(3.3.KK) 

"SAVE T<£.3.KK) 

"SAVE T(3.3,KK> 


" CALCULATE HKL(I.J) = HKL(I.J> + TKMUGKLTLC M. J > 


LDSPI SIX; DB=6. 

LDSPI THREE; DB=3. 

LDSPI BASEGK; DB=46 , 

MOV BASEGK. GKLTLl 
LDSPI BASEHK; DB=10 
MOV BASEHK. HKL 
DEC HKL 

MOV THREE. JCNT 

II 

" BEGIN THE J LOOP CALCULATIONS 

II 

L00P11-- CLR DPA; SETDPA 

MOV GKLTLl. GKLTL2; SETMA 
INC GKLTL2 

MOV GKLTL2.GKLTL3; SETMA 
FMUL DPX<0).MD 
INC GKLTL3; SETMA; FMUL 
FMUL DPX<n.MD; 


"LOAD CONSTANT OF 6 
"LOAD CONSTANT OF 3 
"LOAD GKLTL BASE 
"LOAD GKLTL<M.J) ADDRESS 
"LOAD BASE OF HKL 
"LOAD BASE ADDRESS OF HK 
"LOOP SET UP 
"LOAD OUTER COUNTER 


"GET FIRST SET OF DATA 
"GET GKLTLCM.J) 

"GET GKLTLCE.J) ADDRESS 
"GET GKLTL(£.J) 

"TKMK l.J.LL )»GKLTL1 
"GET GKLTL<3.J) 

"TKMI(£. J.LL )*GKLTL£ 


MOV THREE,ICNT 
DPY<0KFMi FMLIL 
FMUL DPX<£>,MD 
FADD FM,nPY<0); FMUL 
LOP HA MOV GKLTL1,GKLTLI; SETMAi 

FADDi FMUL 
FADD FM,FA^ 

ADD THREE. DPA; 
MOV GkLTL£.GKLTL£; SETMA; 

FADD 

FMUL DPXao.MD 

MOV GKLTL3,GKLTL3i SETMA; 

FMUL 

FMUL DPXU).MD 
DPYCeXFM; FMUL 
FMUL DPX<2>,MD; 

DEC ICNT 

FADD FM>DPY<0); FMUL; 

INC HKL; SETMA; MKFA; 

BNE LOP 11 A 


"LOAD INNER COUNT 
"SAVE GKLTLl PRODUCT 
"TKMK3. J.LL )»GKLTL3 
"GKLTLl PRODUCT ♦ GKLTL2 PROD 
"GET CKLTL<1.J> 

"PUSH 

"GKLTL- 1+2+3 PRODUCTS 
SETDPA "GET NEXT DATA PADS 

"GET GKLTL<2.J) 

"PUSH 

"GKLTL< 1 . J >*TKMK 1 . J. LL > ' 

"GET GKLTL<3.J) 

"PUSH 

"GKLTL< 2, J )*TKMI< 2, J, LL ) 

"SAVE GKLTLl PRODUCT 
"GKLTL3*TKMI<3,J.LL> 

"CHECK INNER LOOP 
"GKLTL 1+2 PRMUCT 
"SAVE RESULT 


" CHECK THE OUTER LOOP AND REAJUST THE ADDRESSES 


ADD SIX. GKLTLl 


"GET GKLTL<M.J+1) 


ADD THREE. HKL 
DEC JCNT 
BEQ C0NT4A 
JMP LOOPll 


"GET HKL<I.J+1> 
"CHECK OUTER LOOP 
"IF DONE. CONTINUE 
"IF NOT. JUMP BACK 


" PERFORM THE 1200 LOOP CALCULATIONS FOR THE SECOND EQUATION 

" HKL<I+3.J) = HKL(I+3.J) + TKMI*GKLTL<M+3. J) 

11 


C0NT4A: MOV BASEHK.HKL 

MOV BASEGK. GKLTLl 
ADD THREE. HKL 
ADD THREE. GKLTLl 
MOV THREE. JCNT 
DEC HKL 

•I 

" BEGIN THE J LOOP CALCULATIONS 

n 

LOPllB^ CLR DPA; SETDPA 

MOV GKLTL 1.GKLTL2; SETMA 
INC GKLTL2 


"LOAD BASE ADDRESS 
"LOAD BASE ADDRESS 
"GET HKL<I+3.J> ADDRESS 
"GET GKLTL(M+3.J) ADDRESS 
"LOAD OUTER COUNTER 
LOOP SET UP 


"GET FIRST SET OF DATA 
"GET GKLTL(M.J) 

"GET GKLTL<5.J) ADDRESS 
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MOV CKLTL2-CKLTL3i SETMA 
FMUL DPX<0),MD 
INC GKLTL3; SETMA; PMUL 
FMUL DPX(1)/MD; 

MOV THREE, ICNT 
DPY<0XFM; FMUL 
FMUL DPX<2),MD 
FADD FM,DPY<0); FMUL 
LOPllC^ MOV GKLTLl/GKLTLl; SETMA; 

FADD; FMUL 
FADD FM,FA; 

ADD THREE, DPA; 
MOV GKLTL2,GKLTL2; SETMA; 

FADD 

FMUL DPX<0),MD 

MOV GKLTL3,GKLTL3; SETMA; 

FMUL 

FMUL DPX<1),MD 
DPY<0XFM; FMUL 
FMUL DPX<2),MD; 

DEC ICNT 

FADD FM,DPY<0>; FMUL; 

INC HKL; SETMA; MKFA; 

BNE LOPllC 


"GET GKLTL<5,J) 

"TKMKl, J,LL )»GKLTL1 
"GET GKLTL(6,J) 

"TKMK2,J,LL )^GKLTL£ 

"LOAD INNER COUNT 
"SAVE GKLTL! PRODUCT 
"TKMK3,U,LL ^tGKLTL3 
" GKLTL 1 PRODUCT + GKLTL2 PROD 
"GET GKLTL<4,J> 

"PUSH 

"GKLTL 1+2+3 PRODUCTS 
SETDPA "GET NEXT DATA PADS 

"GET GKLTL<5,J> 

"PUSH 

"GKLTL< 4, J )!»:TKMI< 1 , J, LL ) 

"GET GKLTL<6,J) 

"PUSH 

"GKLTL< 5, J )*TKMI< 2, J, LL > 

"SAVE GKLTL 1 PRODUCT 
"GKLTL3*TKMK3,J,LL) 

"CHECK INNER LOOP 
"GKLTL 1+2 PRODUCT 
"SAVE RESULT 


" CHECK THE OUTER LOOP AND REAJUST THE ADDRESSES 

M 


ADD SIX, GKLTL 1 
ADD THREE, HKL 
DEC JCNT 
BEQ C0NT4B 
JMP LOP 1 IB 


"GET GKLTL<M,J+1) 
"GET HKL<I,J+1> 
"CHECK OUTER LOOP 
"IF DONE, CONTINUE 
"IF NOT, JUMP BACK 


" PERFORM THE 1200 LOOP CALCULATIONS FOR THE THIRD EQUATION 

II 

" HKL<I,J+3) = HKL(I,J+3> + TKMIAGKLTLCM, J+3) 

•I 


C0NT4B^ LDSPI EIGHTN; DB=18. 
MOV BASEHK,HKL 
MOV BASEGK, GKLTL 1 
ADD EIGHTN, HKL 
ADD EIGHTN, GKLTL 1 
DEC HKL 
MOV THREE, JCNT 


"LOAD CONSTANT 
"LOAD BASE ADDRESS 
"LOAD BASE ADDRESS 
"GET HKL<I,J+3) ADDRESS 
"GET GKLTL<M,J+3) ADDRESS 
"LOOP SET UP 
"LOAD OUTER COUNTER 


BEGIN THE J LOOP CALCULATIONS 
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LOPllO CLR DPA> SETDPA 

MOV CKLTLl,CKLTL2i SETMA 
INC GKLTL2 

MOV GKLTL2,GKLTL3; SETMA 

FMUL DPX(0),MD 
INC GKLTL3; SETMA; FMUL 
FMUL DPX(1),MD; 

MOV THPEE.ICNT 
DPY<OKFM. FMUL 
FMUL OPX(£>,MD 
FADD FM..DPY<0>; FMUL 
LOPllE^ MOV GKLTLl.GKLTLl; SETMA; 

FADD; FMUL 
FADD FM.FA; 

ADD THREE, DPA; 
MOV GKLTL2,GKLTL2; SETMA; 

FADD 

FMUL DPX(0),MD 

MOV GKLTL3,GKLTL3; SETMA; 

FMUL 

FMUL DPX<1),MD 
OPY<0XFM; FMUL 
FMUL DPX<2),MD; 

DEC ICNT 

FADD FM,DPY(0); FMUL; 

INC HKL; SETMA; MKFA: 
PNE LOPllE 


"GET FIRST SET OF DATA 
"GET GKLTL<M,J+3) 

"GET GKLTL<2,J+3) ADDRESS 
"GET GKLTU2,J+3) 

"TKMK1,J,LL )*GKLTL1 
"GET GKLTL(3,J+3> 

"TKMK2, J,LL >tGKLTL2 
"LOAD INNER COUNT 
"SAVE GKLTLl PRODUCT 
"TKMK3, J,LL )CGKLTL3 
"GKLTLl PRODUCT + GKLTL2 PROD 
"GET GKLTL<l,J+3) 

"PUSH 

"GKLTL 1+2+3 PRODUCTS 
SETDPA "GET NEXT DATA PADS 

"GET CKLTL<2,J+3> 

"PUSH 

"GKLTL< 1 , J+3 >»TKMI< 1 , J/ LL > 

"GET GKLTL<3,J+3) 

"PUSH 

" GKLTLC £ , J+3 )»TKM I < 2 , J , LL ) 
"SAVE GKLTLl PRODUCT 
"GKLTL3!»:TKMI<3,J.LL) 

"CHECK INNER LOOP 
"GKLTL 1+2 PRODUCT 
"SAVE RESULT 


CHECK THE OUTER LOOP AND REAJUST THE ADDRESSES 


ADD SIX, GKLTLl 
ADD THREE, HKL 
DEC JCNT 
BEG C0NT4C 
JMP LOP I ID 

PERFORM THE 120O LOOP CALCULATI 


"GET GKLTL<M,J+1) 
"GET HKL<I,J+1) 
"CHECK OUTER LOOP 
"IF DONE, CONTINUE 
"IF NOT, JUMP BACK 

FOR THE FOURTH EQUATION 


HKL<I+3,J+3) = HKL<I+3,J+3) + TKMI*GKLTL< M+3, J+3 ) 


C0NT4C^ LDSPI TUNTYl, DB*£1 . 
MOV BASEHK,HKL 
BASEGK, GKLTLl 
ADU TUNTYl, HKL 
ADD TUNTYl, GKLTLl 
DEC HKL 


"LOAD CONSTANT 
"LOAD BASE ADDRESS 
"LOAD BASE OF GKLTLl 
"GET HKL<I+3,J+3) ADDRESS 
"GET GKLTL<M+3,J+3) ADDRESS 
"LOOP SET UP 
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MOV THREE -XMT 


"LOAD OUTER COUNTER 


" BEGIN THE J LOOP CALCULATIONS 

LOP I IF ^ CLP DP A; SETDPA 

MOV GKLTL1.GKLTL2; SETMA 
INC GKLTL2 

MOV CKLTL£,GKLTL3; SETMA 
FMUL DPX<0)4C' 

INC GKLTL3; SETMA, FMUL 
FMUL DPX<1>,MD. 

MOV THPEE,ICNT 
DPYCOKFM; FMUL 
FMUL DPX(£),MD 
FADD FM,DPY<0); FMUL 
LOOPl£^ MOV GKLTLl/GKLTLl; SETMA; 

FADD: FMUL 
FADD FM,FA; 

ADD THREE, DPA; 
MOV CKLTL£,GKLTL£; SETMA; 

FADD 

FMUL DPX<0),MD 

MOV GKLTL3,GKLTL3; SETMA; 

FMUL 

FMUL DPX<1>,MD 
DPY<0XFM; FMUL 
FMUL DPX<£XMD; 

DEC ICNT 

FADD FM,DPY<0); FMUL; 

INC HKL; SETMA; MKFA, 
BNE LOOPl£ 


"GET FIRST SET OF DATA 
"GET GKLTL<M+3,J+3> 

"GET GKLTU5,J+3) ADDRESS 
"GET GKLTL<5,J+3> 

"TKMIv 1, J,LL)»GI'LTL1 
"GET GKLTL<:S, J+3> 

"TKMI<£, J,LL>*GKLTL£ 

"LOAD INNER COUNT 
"SAVE GKLTLl PRODUCT 
"TKMK3,J,LL )»GKLTL3 
"GKLTLl PRODUCT + GKLTL2 PROD 
"GET GKLTL<6,J+3+3) 

"PUSH 

"GKLTL l+£+3 PRODUCTS 
SETDPA "GET NEXT DATA PADS 

"GET GKLTL<5,J+3) 

"PUSH 

"GKLTL< 6, J+3+3 )»TKMI< 1 , J, LL ) 
"GET GKLTL(6,J+3) 

"PUSH 

" GKLTL< 5 , J+3 )*TKM I < £ , J , LL ) 
"SAVE GKLTLl PRODUCT 
"GKLTL3».TKMI<3,J,LL) 

"CHECK INNER LOOP 
"GKLTL l+£ PRODUCT 
"SAVE RESULT 


" CHECK THE OUTER LOOP AND REAJUST THE ADDRESSES 

H 


ADD SIX, GKLTLl 
ADD THREE, HKL 
DEC JCNT 
BEQ C0NT5 
JMP LOPllF 


"GET GKLTL<M,J+1) 
"GET HKL<I,J+1) 
"CHECK OUTER LOOP 
"IF DONE, CONTINUE 
"IF NOT, JUMP BACK 


" FORTRAN: IF <MAP(N) CT, 0) GOTO 1400 

H 


C0NT5 LDDPA; DB*1£ 

LDSPI MAPBAS, DB=DPY<:0> 
LDSPI N; DB=DPX<£) 

INC N 


"GET HIGH DATA PADS 
“RESTORE BASE OF MAP 
"RESTORE N 
“N*N+1 


1 
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ADO» N,MAPBAS, SETMA 
MOV N.N; 0PXC2XSPFN 
LODPA, DB=e 
LOSPI MAP; DB«MO 
CLP ZERO 
SUB MAP, ZERO 
BGT COHTc. 

JMP C0NT7 


GET MAP<H> 

“SAVE N 

“GET NORMAL DATA PADS 
“SAVE MAP<H) 

"SET 2ERO»0 
“TEST MAP(N) .CT.0 
"CONTINUE- IF FALSE 
“OTHERWISE. GOTO 1400 


" PERFORM THE 1300 LOOP CALCULATIONS 

" FORTRAN : DO 1300 J=£.6 

“ JM»J-1 

" DO 1300 1*1. JM 

“ ElJ*Hf<La. 

■ HKL<I.J)*HKL<J,I) 

" 1300 HKL<J.1>*E1J 


C0NT6= MOV BASEHK.HKTMP 
LDSPI SIX; DB*6. 

LDSPI J; DB-1. 

N 

" BEGIN THE OUTER LOOP 

M 

L00P13^ MOV J.JM 

MOV BASEHK.HKL2 
ADD SIX.WCTMP 
MOV HKW.HKLl 
DEC HKLl 
ADD J.HKL2 
SUB SIX.HKL2 
INC J 

II 

" PERFORM THE INNER LOW> 

ft 

L00P14^ NOP 

ADD SIX.HF^LE; SETMA 
NOP 

IKS HKLl, SETMA 
NOP 

MOV HKLl. HKLl; SETMA, MKMO 
DEC JM 

MOV HKL2.HKL2; SETMA; MKFtt); 
BNE LOOP 14 

M 

" TEST OUTER LOOP 


“LOAD BASE 
"LOAD CONSTANT 
"LOAD OUTER COUNTER - 1 


“LOAD THE INNER COUNTER 
“LOAD BASE 

"GET HKL<I,J) ADDRESS 
"LOAD HKLU.J) 

"LOOP SET UP 
"GET HKL(J,I) ADDRESS 
"LOOP SET UP 
“GET J 


“GET HKL<J,I) ( 

j I 
1 

“GET HKLU.J) 'j 

"STORE HKL<I,J)»HKL<J.I) ti 

i 

“HKL<J.I>*HKL(I.J> j 

“LOOP DONE? 1 

I 


4i 


SUBi SIX,J 
BED COMT-’ 
JMF LOOP 13 


"TEST j«e 
IF EQUAL. END 


FOPTPAH 1400 L0C-1ABS<MAP<N>> 


C0MT7 LDSPI SIX. DB*6 
LDSPI HE7; DB«27 
LDDPA; DB*l£. 

LDSPI N, DB=DPX<2) 

LDSPI MAPBAS; DB«DPY(0> 
ADD# N, MAPS AS; SETMA 
LDDPA; DB*0. 

CLR ZERO 

LDSPI LOC; DB»MD 
MOV LOC.LOC 
BGE C0NT8 
SUB LOC.ZERO 
MOV ZERO.LOC 


"LOAD CONSTANT 
"LOAD CONSTANT 
"GET HIGH DATA PADS 
"RESTORE N 

"RESTORE BASE OF MAP ARRAY 
"GET MAP<>N> 

"GET NORMAL DATA PADS 
"GET e 

"LOC»MAP<N> 

"IS LOC NEGATIVE 
"IF POSITIVE. CONTINUE 
"GET -LOC 
"STORE -LOC 


" CALCULATE THE BASE ADDRESS OF Hd.J.LOC) 


CONTS^ DEC LOC 

MOV LOC.LOC; DPX<0XSPFN 

FADD ZERO.MDPX<0); MOV N27.N2? 

FADD; 

RPSF THYSIX 

DPY<0KDB 
FMUL DPYCO.FA 
FMUL 
FMUL 

DPXCOXFM 
FIX DPX<0> 

FADD 

C^X<0XFA. 

LDSPI 1«ASE; 08*442. 

LDSPI HADDR; DB*DP‘X<0) 


"GET uOC -1 
"FLOAT LOC-1 


"DPY<0)*36. 

"36CL0C-1 


"GET INT(36»L0C-1 ) 

"STORE RESULT 
"LOAD H BASE 

"GET H<1. l.LOC; ADDRESS 


" FORTRAN^ IF<NDF LT 6) GOTO 1600 


" THE ADDRESS OF H( 1.1. LOC) HAS BEEN CALCULATED FOR USE 
" IN EITHER THE 1500 OR 1700 LOOPS 


SUBi SIX.IOF 
BGE CC»<T9 
JC CO»m0 


TEST FW? hOF LT 6 
IF TPL€. CONTINUE 
OTHERSUIICE, GOTO 1600 
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M 


" BEGIN THE 15O0 LOOP FORTRAN 


FORTRAN: DO 1500 J=1.6 
" DO 1500 1=1/6 

■■ 1500 H<I/J/LOC) = HC 1/ J/LOC) + HKL( I/J) 

'■ GOTO £000 


C0NT9: ADD HBASE/HADDR; SETMA 

MOV HADDR/HADOR£ 

MOV BASEHK/HKLl; SETMA 
DPX<0KMD; 

DEC HADDR2 
LDSPI CNT; 0B=36. 


"GET HU/l/LOC) 

"LOAD DESTINATION ADDRESS 
"GET HKL<1/1> 

"STORE H<l/l/LOC) 

"LOOP SET UP 


" PERFORM THE 1500 LOOP CALCULATIONS 


L00P15-- FADD DPX<0),MD 

INC HADDR; SETMA; 

FADD 

INC HKLl; SETMA 
DPX<0XMD/ 

DEC CNT 

INC HADDR£; SETMA; MKFA; 

BNE LOOP15 


"H( I /J/LOC >+HKL<I/J) 
"GET NEXT H< I/J/ LOO 
"PUSH 

"GET NEXT H<I/J) 

"STORE H<I,J/LOO 
"TEST LOOP 
"SAVE RESULT 
"FINISHED? 


JMP CONTI 1 "GOTO £000 

II 

" PERFROM THE 1700 LOOP CALCULATIONS 

II 

" FORTRAN: 1600 DO 170© 1=1 /NDF 
" NROW=INEX< I ) 

" DO 1700 J=1/NDF 

" NCOL=INEX<J) 

" 1 700 H< I / J / LOC >=H< I / J / LOC )+HKL< NROW / NCOL ) 


CONT10: LDSPI SIX; DB=6. 

CLR I HEX 
CLR I 

ADD HBASE/ HADDR 
DEC HADDR 

LDSPI TMA; DB=TADR2 


"LOAD CONSTANT OF 6. 

"GET I HEX BASE ADDRESS 
"CLEAR OUTER COUNTER 
"GET H<1/ 1/LOO ADDRESS 
"GET H(1/1/LOO ADDRESS -1 
"SET ROM BASE ADDRESS 


OUTER LOOP 
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LOOPie^ ADD# I,lNEXi SETMA 
MOV BASEHK . HKADDP 
DEC HKADDF 
LDSPl NROWi DB=MD 
ADD HROW* HKADDP 
INC HADDR 
MOV HADDR. HADDR 1 


“GET INEX<I) 

"GET HKL<1. l> ADDRESS 
"LOOP SET UP 
"STORE NROUI=iNEX<I) 

"GET HKLCNROUI. 1> ADDRESS 
"GET Ha+l.J.LOC> ADDRESS 
"LOAD H ADDRESS 


" INNER LOOP SET UP 


CLR J 

ADD# J. INEXi SETMA 
NOP 

MOV HADDR I.HADDR2; SETMA 
LDSPl NCOLj DB=MD 
ADD# NCOL.TMA; SETMA 
DFX<0XMD. 

SUB SIX.HADDR2 

NOP 

LDSPl HKINC; DB=MD 
ADD# HKINC.HKADDR; SETMA 
INC 1 
NOP 


"CLEAR INNER COUNT 
"GET INEXCJ) 

"LOAD TARGET ADDRESS 
"NCOL=INEX<J) 

"GET INCREMENT FROM ROM 
"STORE H(NCOL.l) 

"LOOP SET UP 
WAIT FOR M D. 

"STORE TABLE INCREMENT 
"GET HKLCNROUI.NCOL) 
"INC OUTER LOOP COUNT 


" INNER LOOP CALCULATIONS 


LOOP17^ INC J.i 

FADD DPX<0 > .MD 
ADD# J.. INEX.; SETMA; 

FADD 

NOP 

ADD SIX.HADDRl; SETMA 
LDSPl NCOL; DB=MD 
ADD# NCOL.TMA; SETMA 
NOP 

DPXv'OKMD 
LDSPl HKINC; DB=MD 
ADD# HKINC. HKADDRi SETMA 
SUB# NDF.J 

ADD SIX.HADDR2; SETMA; MKFA; 
BNE LOOP 17 


"INCREMENT INNER LOOP COUNT 
"H< I . J. LOC >+HKL< NROW. NCOL > 
"GET INEX<J> 

"PUSH 

"GET Ha.J+l.LOC) 

"STORE NCOL=INEX<J) 

"GET INCREMENT IN ROM TABLE 

"SAVE HU. J+l.LOC) 

"SAVE INCREMENT 
"GET HKL< NROW. NCOL) 

"TEST INNER LOOP 
"STORE RESULT 


H 


SUB# NDF.I 
BEQ CONTU 
JMP LOOP I 6 


TEST OUTER LOOP 
IF DONE, CONTINUE 
IF NOT. GBRANCH BACK 
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CONTll LDDPA. DB=l£ 

LDSPI L. DB=DPX<8> 
LDSPI Ki DB=DPX<1) 
SUB<< K .L 
BEQ C0HT12 
JMP LOOP£ 

COini£ SUB»« L-NNODES 
BEQ CONTI 3 
JMP LOOP I 
TH^SIX^ itFP 36 
CONTI 3^ RETURN 
»ENO 


'•GET HIGH DATA PADS 
RESTORE L 
RESTORE K 

IS INNER K LOOP DONE 
IF YES. CONTINUE 
OTHERWISE, BRANCH BACK 
TEST OUTER MOST £000 LOOP 
"IF FINISHED, END 
OTHERWISE, BRANCH BACK 
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APPENDIX J 


Listings of: 


FUSEL 

LUT 


(XQT TAB 

START eo, 5$ ROTATIONS ABOUT Y EXCUIUED 

TITLE" FUSELAGE MODEL, PSPARl 

TEXT 

" MEMBRANE-ROD-BEAM FUSELAGE MODEL 
"NONREPEATABLE PART 

JLOC5 FUSELAGE DIA. BOO. CM. , LENGTH-BOO. CM. 

FORMAT- 2 SCYLINDRI CAL COORDINATES 
1 400. 0. 0. 400. 337. S 0. 16 1 S 

16 400. 0. 600. 400. 337.5 600. 

NREF 

FCRMAT-2 

1 -2 O. 0. lOOOOOOO. 

2 1 0. 0. lOOOOOOO. 

MATC 

1 .7246 0.3 .0026 22.-6$ AL-ALLOY, METRIC UNITS 
E23 SECTKHi PROPERTIES $ROD ELEMENTS 
1 4.168$AREA OF TOE RODS 
SHELL SECTI(»i PROPERTIES 
1 O.ISSKIN THICKNESS 
E21 SECTION PROPERTIESSBEAN ELEMENTS 
DSY 1 16804. 0. 1262.7 0. 108. 144. 0. 6.0764 O. 0. 

0. O. 0. -8.7776 17. 3.2222 17. 3.2222 -17. -8.7776 -17. 
CONSTRAINT CASE 1 

ZERO 1,2,3;1,16$ CANTILEVER THE FUSELAGE 
[XQT ELD 

E23$ROD ELEMENTS 

NSECT-IS 

NREF-2 

1 17 1 4 3 1 $ 

4 20$ 

5 21$ 

6 22$ 

52 68$ 

53 69$ 

54 70$ 

7 23 1 4 10 1$ 

E41$ MEMBRANE PANELS 
NSECT-IS 

1 17 18 2 2 16 1 $ 

49 65 66 50 2 16 1$ 

17 33 34 16 1 1 2 2 16$ 

23 39 40 24 1 1 9 2 16$ 

32 48 33 17$ 

48 64 49 33$ 

E21 

NSECT-1$ 

NREF-1 

1 2 2 16 2 16 $ 

49 SO 2 16 2 16$ 


33 34$ 

34 35$ 

39 40$ 

40 41$ 

41 42$ 

42 43$ 

43 44$ 

44 45$ 

45 46$ 

46 47$ 

47 48$ 

48 33$ 

[XQT E 
[XQT EKS 
[XQT TQPO 
[XQT K 
[XQT ZNV 
[XQT AOS 

SYSVEC;APPLIED FORCES 1 
Z-l;J-65; -lOOOO. 

SYSVEC) APPLIED FORCES 2 
I-l;J-69>77> -20000. 20000. 
SYSVEC; UNIT VEC 
I-l; J-1,80; 1.0 
DEFINE WT-DEM DIAG 0 O 
DEFINE UN-DNIT VEC 
OBJF AUS 1 l-XTY(ON,»fr) 

[XQT DCU 

PRINT 1 OBJF AUS 1 1 
[XQT SSOL 
[XQT GSF 
[XQT PSF 
[XQT VPRT 

PRINT APPL FORC 1 1 
PRI T STAT DISP 1 1 
[XQT SSOL 
RESET SET-2 
[XQT OSF 
RESET SET-2 
[XQT PSF 
RESET SET-2 
[XQT VPRT 

PRINT APPL FORC 2 1 
PRINT STAT DISP 2 1 
[XQT EXIT 
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(TQ TAB . GENERATE BASIC TABIXS DEFINING STRUCTURE 

START 372$ 

TITLE"SATURN V LAUNCHER IMBILICAL TONER (LOT) 

TEXTS 

"SATURN V LAUNCHER UMBILICAL TONER (LOT) 

M 

"OAT ESIGNEOS-i'O EXEMPLIFY THE USE OP "INC" 

"AND "IK>D" COMMANDS IN THE "ELD" PROCESSOR. 

MATERIAL CONSTANTSS 

$ MATERIAL PROPERTY IS DEFINED NITH A 

$ NEIGHT DENSITY 

1 3.-f7 .29962668 .28$ 

CONSTRAINT CASE 1$ 

$ JOINTS 1-4 ARE COMPLETELY CONSTRAINED 

ZERO 12345 6| 14 1$ 

BEAM ORlENTATICm SPECIFICATIONSS 


1 1 : 

i 1 

o.$ 






2 1 : 

L 1 

o.$ 






E21 SECTI(»I PROPERTIESS 





GIVN : 

L 9012 . 1 

0 

250.4 

0 

44.16 9.47$ 

GIVN 

2 

2096.4 

0. 

76.5 

0. 

22.4 

8.4 $ 

GIVN 

3 

63.35 

0. 

63.35 

0. 

7.3 

126.7 $ 

GIVN 

4 

4461. 

0. 

135.1 

0. 

31.77 

4.66$ 

GIVN 

5 

1814.5 

0. 

63.8 

0. 

20. 

1.7 $ 

GIVN 

6 

9012.1 

0. 

250.4 

0. 

44.16 

9.47$ 

GIVN 

7 

2096 . 4 

0. 

76.5 

0. 

22.4 

8.4 $ 

GIVN 

8 

28.14 

0. 

28.14 

0. 

5.58 

56 . 28$ 

GIVN 

9 

2824.6 

0. 

95.7 

0. 

24.71 

2.57$ 

GrVN 

10 

1140.7 

0. 

44. 

0. 

16.18 

1.14$ 

GIVN 

11 

1140.7 

0. 

44. 

0. 

16.18 

1.14$ 

GIVN 

12 

21.7 

0. 

2.89 

0. 

3.53 

.08$ 

GIVN 

13 

105.3 

0. 

2.79 

0. 

4.86 

.1 $ 

GIVN 

14 

446.3 

0. 

22.1 

0. 

10.59 

.5 $ 

GIVN 

15 

15.16 

0. 

15.16 

0. 

4.3 

30.32$ 

GIVN 

16 

21200. 

0. 

21200. 

0. 

276. 

31100. $ 

GIVN 

17 

475.7 

0. 

475.7 

0. 

26.04 

951.4 $ 

GIVN 

18 

3988.6 

0. 

116.9 

0. 

29.11 

3.47$ 

GIVN 

19 

16667. 

0. 

16667 . 

0. 

228. 

24554. $ 

GIVN 

20 

475.7 

0. 

475.7 

0. 

26.04 

951.4 $ 

GIVN 

21 

3988.6 

0. 

116.9 

0. 

29.11 

3.47$ 

GIVN 

22 

14284. 

0. 

14284. 

0. 

205. 

21538. $ 

GIVN 

23 

475.7 

0. 

475.7 

0. 

26.04 

951.4 $ 

GIVN 

24 

3266.7 

0. 

115.1 

0. 

27.65 

3.8 $ 

GIVN 

25 

12609. 

0. 

12609 . 

0. 

182. 

18683. $ 

GIVN 

26 

361.5 

0. 

361.5 

O. 

19.24 

723. $ 

GIVN 

27 

3266.7 

0. 

115.1 

0. 

27.65 

3.8 $ 

GIVN 

28 

10773. 

0. 

10773. 

0. 

160. 

16000. $ 

GIVN 

29 

361.5 

0. 

361.5 

0. 

19.24 

723. $ 

GIVN 

30 

3266.7 

0. 

115.1 

0. 

27.65 

3.8 $ 

UiVN 

31 

9068. 

0. 

9068. 

0. 

138. 

13476. $ 


GIVN 

32 

361. S 

0. 

361.5 

0. 

19.24 

723. $ 

GTVN 

33 

3266.7 

0. 

115.1 

0. 

27.65 

3.8 $ 

GIVN 

34 

7459. 

0. 

7459. 

0. 

117. 

11233. $ 

GIVN 

35 

361.5 

0. 

361.5 

0. 

19.24 

723. $ 

GIVN 

36 

2364.3 

0. 

86.3 

0. 

24.7 

3.57$ 

GIVN 

37 

5969. 

0. 

5969. 

0. 

96. 

8917. $ 

GI>m 

38 

361.5 

0. 

361.5 

0. 

19.24 

723. $ 

GIVN 

39 

2364.3 

0. 

88.3 

0. 

24.7 

3.57$ 

GIVN 

40 

4585. 

0. 

4585. 

0. 

76. 

6859. $ 

GIVN 

41 

248.5 

0. 

248.5 

0. 

12.88 

497. $ 

GIVN 

42 

2096 . 4 

0. 

76.5 

0. 

22.4 

8.4 $ 

GIVN 

43 

2402.4 

0. 

930.1 

0. 

56.73 

34.4S$ 

GIVN 

44 

2«8.5 

0. 

248.5 

0. 

12.88 

497. $ 

GIVN 

45 

2096.4 

0. 

76.5 

0. 

22.4 

8.4 $ 

GIVN 

46 

2402.4 

0. 

930.1 

0. 

56.73 

34.45$ 

GIVN 

47 

248.5 

0. 

248.5 

0. 

12.88 

497. $ 

GIVN 

48 

1814.5 

0. 

63.8 

0. 

20. 

1.7 $ 

GIVN 

49 

1266.5 

0. 

454.9 

0. 

32.65 

7.63$ 

GIVN 

50 

192.3 

0. 

192.3 

0. 

9.84 

384.6 $ 

GIVN 

51 

1814.5 

0. 

63.8 

0. 

20. 

1.7 $ 

GIVN 

52 

1266.5 

0. 

454.9 

0. 

32.65 

7.63$ 

GIVN 

53 

192.3 

0. 

192.3 

0. 

9.84 

384.6 $ 

GIVN 

54 

1814.5 

0. 

63.8 

0. 

20. 

1.7 $ 

GIVN 

55 

641.5 

0. 

107.3 

0. 

17.94 

2.01$ 

GIVN 

56 

192.3 

0. 

192.3 

0. 

9.84 

384.6 $ 

GIVN 

57 

1814.5 

0. 

63.8 

0. 

20. 

1.7 $ 

GIVN 

58 

641.5 

0. 

107.3 

0. 

17.94 

2.01$ 

GIVN 

59 

192.3 

0. 

192.3 

0. 

9.84 

384.6 $ 

GIVN 

60 

1814.5 

0. 

63.8 

0. 

20. 

1.7 $ 

GIVN 

61 

2987.3 

0. 

203.5 

0. 

29.43 

4.52$ 

GIVN 

62 

1140.7 

0. 

44. 

0. 

16.18 

1.14$ 

GIVN 

63 

1326.8 

0. 

53.1 

0. 

18.23 

1.71$ 

GIVN 

64 

446.3 

0. 

22.1 

0. 

10.59 

.5 $ 

GIVN 

65 

15.16 

0. 

15.16 

0. 

4.3 

30.32$ 

GIVN 

66 

5886.9 

0. 

170.3 

0. 

34.71 

4.9 $ 

GIVN 

67 

2987.3 

0. 

203.5 

0. 

29.43 

4.52$ 

GIVN 

68 

2987.3 

0. 

203.5 

0. 

29.43 

4.52$ 

GIVN 

69 

2987.3 

0. 

203.5 

0. 

29.43 

4.52$ 

GIVN 

70 

2987.3 

0. 

203.5 

0. 

29.43 

4.52$ 

GIVN 

71 

15.16 

0. 

15.16 

0. 

4.3 

30.32$ 

GIVN 

72 

26478. 

0. 

26478. 

0. 

326. 

38584. $ 

GIVN 

73 

562. 

0. 

562. 

0. 

18.41 

1124. $ 

GIVN 

74 

732. 

c. 

732. 

0. 

24.35 

1464. $ 

GIVN 

75 

1157. 

0. 

1157. 

0. 

40.19 

2314. $ 

GIVN 

76 

26300. 

0. 

26300. 

0. 

294. 

38700. $ 

GIVN 

77 

562. 

0. 

562. 

0. 

18.41 

1124. $ 

GIVN 

78 

1556. 

0. 

1556. 

0. 

56.6 

3112. $ 


JOINT LOCATIONS^ 

1 360. 907. 

2 -360. 907. 


-2880. 270. 

-2880. -270. 


589. -2160. 2 48$ 
589. -2160. 2 48$ 


3 

-360. 

-427, -2880. 

-270. 

-109. 

-2160. 2 

48$ 

4 

360. 

-427. -2880. 

270. 

-109. 

-2160. 2 

48$ 

5 

316. 

348. -2520. 

270. 

348. 

-2160. 2 

24$ 

6 

-316. 

348. -2520. 

-270. 

348. 

-2160. 2 

24$ 

7 

-316. 

132. -2520. 

-270. 

132. 

-2160. 2 

24$ 

8 

316. 

132. -2520. 

270. 

132. 

-2160. 2 

24$ 

9 

316. 

480. -2520. 

270. 

480. 

-2160. 2 

24$ 

10 

-316. 

480. -2520. 

-270. 

480. 

-2160. 2 

24$ 

11 

-316. 

0. -2520. 

-270. 

0. 

-2160. 2 

24$ 

12 

316. 

0. -2520. 

270. 

0. 

-2160. 2 

24$ 

13 

108. 

480. -2520. 

108. 

480. 

-2160. 2 

24$ 

14 

-108. 

480. -2520. 

-108. 

480. 

-2160. 2 

24$ 

15 

108. 

348. -2520. 

108. 

348. 

-2160. 2 

24$ 

16 

-108. 

348. -2520. 

-108. 

348. 

-2160. 2 

24$ 

17 

108. 

132. -2520. 

108. 

132, 

-2160. 2 

24$ 

18 

-108. 

132. -2520. 

-108. 

132. 

-2160. 2 

24$ 

19 

108. 

0. -2520. 

108. 

0. 

-2160. 2 

24$ 

20 

-108. 

0. -2520. 

-108. 

0. 

-2160. 2 

24$ 

21 

316. 

240. -2520. 

270. 

240. 

-2160. 2 

24$ 

22 

0. 

751. -2520. 

0. 

589. 

-2160. 2 

24$ 

23 

-316. 

240. -2520. 

-270. 

240. 

-2160. 2 

24$ 

24 

0. 

-271. -2520. 

0. 

-109. 

-2160. 2 

24$ 

25 

316. 

751. -2520.8 





26 

-316. 

751. -2520.8 





27 

-316. 

-271. -2520.8 





28 

316. 

-271. -2520.8 





S3 

240. 

348. -1920. 

240. 

348. 

1440. 15 

20$ 

54 

108. 

480. -1920. 

108. 

480. 

1440. 15 

20$ 

55 

-108. 

480. -1920. 

-108. 

480. 

1440. 15 

20$ 

56 

-240. 

348. -1920. 

-240. 

348. 

1440. 15 

20$ 

57 

-240. 

132. -1920. 

-240. 

132. 

1440. IS 

20$ 

58 

-108. 

0. -1920. 

-108. 

0. 

1440. 15 

20$ 

59 

108. 

0. -1920. 

108. 

0. 

1440. 15 

20$ 

60 

240. 

132. -1920. 

240. 

132. 

1440. 15 

20$ 

61 

108. 

348. -1920. 

108. 

348. 

1*40. 15 

20$ 

62 

-108. 

132. -1920. 

-108. 

132. 

1440. 15 

20$ 

63 

-108, 

348. -1920. 

-108. 

348. 

1440. 15 

20$ 

64 

108. 

132. -1920. 

108. 

132. 

1440. 15 

20$ 

65 

240. 

240. -1920. 

240. 

240. 

1680. 16 

20$ 

66 

0. 

480. -1920. 

0. 

480. 

1680. 16 

20$ 

67 

-240. 

240. -1920. 

-240. 

240. 

1680. 16 

20$ 

68 

0. 

0. -1920. 

0. 

0. 

1680. 16 

20$ 

69 

240. 

480. -1920. 

240. 

480. 

1680. 16 

20$ 

70 

-240. 

480. -1920. 

-240. 

480. 

1680. 16 

20$ 

71 

-240. 

0. -1920. 

-240. 

0. 

1660. 16 

20$ 

72 

240. 

0. -1920. 

240. 

0. 

1680. 16 

20$ 

353 

240. 

290. 1680.8 





354 

50. 

480. 1680.8 





355 

-50. 

480. 1680.8 





356 

-240. 

290. 1680.8 
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357 

-240. 

190. 

1680.3 

358 

-50. 

0. 

1680.3 

359 

50. 

0. 

1680.3 

360 

240. 

190. 

1680.3 

361 

50. 

290. 

1680.3 

362 

-50. 

190. 

1680.3 

363 

-50. 

290. 

1680.3 

364 

50. 

190. 

1680.3 


RMASSS 

$ 


REPEAT 
49 100 
REPEAT 
89 43 

REPEAT 
109 24 

REPEAT 
129 47 

REPEAT 
149 30 

REPEAT 
169 64 

REPEAT 
189 56 

REPEAT 
209 84 

REPEAT 
229 72 

REPEAT 
249 110 
REPEAT 
269 116 
REPEAT 
289 121 
REPEAT 
309 106 
REPEAT 
329 151 
REPEAT 
349 231 
REPEAT 
369 308 
[XQT ELD 
E21 3 
GROUP 1 

NSECT»1 


4 13 
0253 
4 13 
8803 
4 13 
7883 
4 13 
4803 
4 13 
6323 
4 13 
7103 
4 13 
9383 
4 13 
4603 
4 13 
4723 
4 13 
5853 
4 13 
0483 
4 13 
7503 
4 13 
5983 
4 13 
3783 
4 13 
4783 
4 13 
2823 


"FLOOR 2 
9 25 
5 9 

5 21 
8 21 

6 23 


1 14 1 

1 14 1 

1 12 2 


RIGID LUMPED MASSES 


. READ ELEMENT DEFINITICWS 


3 

3 

3 

3 

3 


I 


NSECT-2 


NSECT-3 ; 


NSECT-4 


NSECT-5 ; 


22 

26 

1 

1 

2 

2 


22 

25 

1 

1 

2 

2 


13 

25 

1 

1 

2 

1 


13 

22 

1 

1 

2 

6 


14 

22 






20 
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